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FOREWORD 

This  report.  Volume  II,.  is  one  of  a  series  of  three  volumes,  reporting 
work  accomplished  under  Contract  AF19(268)-239  initiated  in  December  1961. 

The  purpose  of  this  program  was  to  investigate  by  theory  and  experiment 
the  pertinent  factors  influencing  the  creation,  propagation  and  detection  of 
a  type  of  magnetohydrodynamic  wave,  the  Alfven  wave,  that  could  exist  in  a 
partially  ionized,  low  density  medium. 

The  theoretical  aspects  of  magnetohydrodynamic  waves,  as  might  be 
generated  under  ionospheric  conditions  and  excited  in  a  simulated  environment, 
are  investigated  in  this  report.  Experimental  design  and  subsequent  test 
results  of  wave  propagation  in  three  McDonnell  plasma  facilities,  are  reported 
in  Volume  III.  Volume  I  of  this  series  provides  the  reader  with  a  brief 
summary  of  Volumes  II  and  III. 

This  program  required  the  cooperation  and  assistance  of  many  persons. 

The  chief  contributors,  in  addition  to  C.  D.  Joerger,  Project  Leader  were: 

J.  L.  Hickerson  and  E.  W.  Hobbs  for  plasma  wave  theory;  T.  R.  McPherron, 

E.  S.  Thompson,  G.  L.  Elder,  and  J.  L.  Walker  for  experimental  design  and 
execution;  and  H.  J.  Fivel  and  J.  N.  Holsen  for  thermodynamic  analyses. 
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1.  SUMMARY 

This  report  is  submitted,  in  fulfillment  of  Contract  AF19(268)-239  for  "the 
investigation  of  magnetohydrodynamic  waves.  Under  this  contract  McDonnell  has 
successfully  generated  Alfven  waves  in  an  arc  discharge  tube  and  further,  has 
evidence  that  this  phenomena  was  also  detected  in  the  ionized  flow  of  a 
Hypervelocity  Impulse  Tunnel  (HIT).  The  goal  of  this  program  was  to  determine 
by  means  of  theory  and  experiment  the  pertinent  factors  influencing  the  creation, 
propagation,  and  detection  of  the  Alfven  wave  that  can  exist  in  a  partially 
ionized,  low  density  atmosphere.  This  volume  discusses  the  theoretical  studies 
conducted  during  this  program  to  describe  plasma  wave  propagation  in  the  iono¬ 
sphere  and  in  a  laboratory  generated  plasma  environment. 

The  equations  that  describe  the  wave  motion  in  a  plasma  have  been  treated 
by  a  number  of  authors;  however,  solutions  are  available  only  for  a  limited 
number  of  cases,  treating  conditions  in  which  coupling  between  the  ion  and  the 
neutral  gases  is  either  very  small  or  very  large  for  the  ionosphere,  and  whereas 
for  laboratory  plasmas  the  coupling  is  usually  intermediate.  The  inability  to 
generate  a  plasma  of  even  modest  size  has  prevented  other  experimenters  from 
trying  to  scale  laboratory  plasmas  to  ionospheric  phenomena.  The  effect  of 
neutral  particle  damping,  which  is  a  major  contributor  to  Alfven  wave  attenua¬ 
tion  in  the  ionosphere  has  not  been  previously  studied  because  of  the  severe 
attenuation  of  the  wave  in  the  experimental  facilities  (the  mere  presence  of  a 
large  percentage  of  non- ionized  gas  rapidly  attenuates  the  wave). 

Previous  major  experimental  investigations  have  treated  the  thermonuclear 
and  not  the  ionospheric  problem.  Important  requirements  for  ionospheric 
simulation  are  availability  of  a  facility,  large  in  size,  generation  of  a 
realistic  percentage  of  gas  ionization,  and  a  low  ion  neutral  collision  frequency. 

L 


MCI 


1 


REPORT  NO.  A219 
30  NOVEMBER  1963 


IV1 H  D  Wave  Investigation 

It  -  theory 

It  was  hypothesized  that  these  parameter  requirements  could  he  met  in  the  flow 
of  the  Hypervelocity  Impulse  Tunnel  (HIT)  under  high  arc  chamber  temperatures. 
This  type  of  facility  had  previously  been  used  only  for  aerodynamic  testing 
in  non-ionized  gas  flow. 

The  first  part  of  this  volume  discusses  the  theory  of  plasma  waves  based 
upon  a  three  fluid  description  (neutral,  ion,  and  electron  gases).  A  general 
dispersion  equation  was  derived  and  then  specialized,  to  enumerate  the  properties 
of  transverse  and  longitudinal  waves  propagating  along  and  across  magnetic  field 
lines .  Additionally,  special  solutions  of  the  equation  were  used  to  describe 
propagation  in  a  plasma  waveguide  and  to  enumerate  the  existence  criteria  for 
Alfven  waves  in  the  ionosphere,  HIT,  and  plasma  waveguide. 

The  complete  dispersion  relationship  for  the  Alfven  wave  was  programmed 
on  an  IBM  709U  computer.  Parametric  data  was  obtained  for  the  attenuation 
distance  and  phase  velocity  of  the  wave  as  a  function  of  disturbance  frequency, 
ion  concentration  and  ion  percentage,  magnetic  field  strength,  temperature,  and 
molecular  weight.  These  solutions  were;  found  to  agree  with  the  asymptotic 
solutions  obtained  by  other  authors.  For  a  low  frequency  disturbance  the  Alfven 
wave  has  long  attenuation  distances  and  propagates  with  a  phase  velocity  deter¬ 
mined  by  the  total  density  of  the  gas.  As  the  coupling  between  the  ions  and 
neutrals  breaks  down,  the  phase  velocity  increases  to  a  maximum  determined  by 
the  ion  density.  The  computed  attenuation  distances  did  not  indicate  the 
increased  damping  predicted  for  the  intermediate  coupling  case;  instead  it 
indicated  that  the  phase  velocity  of  the  wave  increased  so  as  to  compensate  for 
the  shorter  time  of  wave  existence.  The  range  of  frequencies  for  which  the 
dispersion  equation  is  valid,  includes  the  ion  cyclotron  frequency. 

The  second  part  of  this  volume  relates  the  characteristics  of  the  iono¬ 
sphere  to  the  characteristics  of  the  ionized  flow  existing  in  the  initial 
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phases  of  the  HIT  flow,  and  describes  the  properties  of  two  small  plasma  facil¬ 
ities,  an  electromagnetic  shock  tube  and  an  arc  discharge  tube. 

The  Hypervelocity  Impulse  Tunnel  can  be  described  as  a  large  nozzle,  whose 
flow  is  generated  by  arc  heating  a  precharged  volume  of  gas ,  nitrogen  in.  this  case . 
Two  ionized  flow  regions  are  present:  (l)  the  shock  excited  region  created  by 
the  initial  shock  generated  by  the  bursting  of  a  diaphragm  during  the  arc 
heating  of  the  gas,  and  (2)  the  initial  phase  of  blow-down  flow  which  follows 
the  initial  shock. 

The  shock  excited  flow  is  generated  within  a  fraction  of  a  millisecond  of 
the  initiation  of  the  arc  discharge.  As  the  shock  travels  through  the  divergent 
nozzle  it  loses  strength.  Three  analyses  have  been  performed  to  predict  the 
magnitude  of  this  shock  attenuation  and  the  available  test  time.  While  the 
duration  of  ionized  flow  between  the  shock  and  the  initial  gases  expanding  from 
the  arc  chamber  theoretically  lasts  a  fraction  of  a  millisecond,  experi¬ 
ments  have  shown  that  test  times  up  to  2  milliseconds  actually  exist.  The  high 
temperature  low  density  gas  behind  this  shock  is  highly  ionized. 

Theblow-down  period,  with  about  four  milliseconds  of  ionized  flow,  provided 

/ 

the  conditions  in  which  the  Alfven  wave  experiments  were  performed.  At  the 
12,000°K  arc  chamber  temperature  used,  nitrogen  is  highly  dissociated  and  highly 
ionized.  The  ionization  at  the  throat  is  principally  IT*"  ions.  During  the  rapid 
expansion  of  gas  through  the  nozzle  into  the  test  section,  the  recombination  of 
these  ions  through  three-body  interactions  is  not  realized,  and  thus  the  N+  con¬ 
centration  remains  intact.  On  the  other  hand,  the  N2+  ion  rapidly  disappears 
through  dissociative  recombination  and  does  not  exist  in  the  test  section  in  the 
charged  state.  Theoretically,  ion  concentrations  of  greater  than  1011  ions/cm^ 
should  exist  in  the  test  section.  As  the  arc  chamber  temperature  drops,  j , 

the  number  of  N+  ions  created  rapidly  decrease^  and  ionization  in  the  test 


AfCOOWWVLL 


3 


M  M  D  Wave  Investigation  30  NOVEMBER  1963 

1 1  -  theory 

section  ceases.  Two  analyses  were  performed,  to  verify  this  effect,  a  simplified 
frozen  flow  analysis  and  a  non-equilibrium  flow  analysis  using  the  IBM  7094 
computer. 

The  two  smaller  plasma  sources  used  for  experimentation  and  instrumentation 
development  were  an  electromagnetic  shock  tube  and  an  arc  discharge  tube.  The 
electromagnetic  shock  tube  creates  a  moving  plasma  by  discharging  a  short  high 
current  pulse  across  a  coaxial  electrode.  A  piston  of  fully  ionized  hot  gas 
is  propelled  down  the  tube,  compressing  and  heating  the  gas  in  front  while  itself 
losing  strength.  The  processes  in  the  discharge  and  the  acceleration  of  the 
gases  are  not  well  understood;  however,  the  high  ionization  can  be  readily 
measured. 

Similar  to  the  electromagnetic  shock  tube  is  the  second  facility  employed 
in  these  experiments,  an  arc  discharge  tube.  An  arc  is  struck  along  the  axis 
of  an  evacuated  tube,  and  the  experiments  are  performed  while  the  arc  is  still 
burning.  The  resulting  plasma  is  more  highly  ionized  than  that  in  the  shock  tube. 

Although  the  analysis  of  arc  heated  gas  effects  presented  many  theoretical 
difficulties  the  experiments  performed  have  provided  valuable  empirical  data. 


Mcoowwm. 


4 


IV!  H  D  Wave  Investigation  30  NOVEMBER  *963 

II  -  theory 

2.  WAVES  IN  PLASMAS 

2.1  INTRODUCTION .  -  Prior  to  starting  any  basic  experimental  program  a 
theoretical  analysis  is  required  to  provide  the  proper  background  for  designing 
the  experiment.  The  most  logical  starting  point  from  which  to  launch  this 
analysis  is  by  a  thorough  literature  survey  in  order  to  achieve  an  indication 
of  the  existence  or  non-existence  of  the  desired  phenomenon.  For  this  program 
a  realistic  theory  was  sought,  capable  of  predicting  and  describing  general 
plasma  oscillations  in  a  partially  ionized  gas.  A  review  of  pertinent  litera¬ 
ture  revealed  the  theoretical  approach  of  Lehnert  (Reference  j-l)  to  be  the  most 
suitable,  since  it  could  be  adapted  to  experiments  in  practically  any  facility. 

In  his  development  Lehnert  considered  a  three -component  gas  consisting  of 
electrons,  ions  and  neutral  particles.  By' means  of  the  basic  force  equations 
(derivable  from  the  Boltzmann  equation),  the  continuity  equations,  and  Maxwell's 
equations,  lie  obtained  a  set  of  linearized  equations  which  he  combined,  after 
specialization,  to  give  a  wave  equation  for  transverse  oscillations  of  the 
magnetic  flux  density.  Later  papers  by  other  authors,  notably  Tanenbaum  and 
Mintzer  (Reference  2-2),  developed  a  somewhat  more  general  dispersion  relation¬ 
ship,  describing  the  various  types  of  wave  motions  that  could  exist  in  partially 

and  fully  ionized  fluids  permeated  by  constant  magnetic  fields.  The  direct  use 

0 

and  adaptation  of  these  works  furnished  a  theoretical  background  for  the 
experimental  program. 

In  the  following  development  a  basic  dispersion  relationship,  valid  for  a 
variety  of  plasma  disturbances,  and  a  reasonably  general  wave  equation  for 
transverse  oscillations  are  derived  using  Lehnert' s  approach.  A  discussion  of 
interesting  specializations  of  the  dispersion  equation  with  emphasis  on  Alfven 
wavea  in  unbounded  media  and  wave  guides  is  carried  out  in  section  2.4.  In 
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section  2.4.4  the  dispersion  relationship  is  reformulated  for  parametric  analysis 
by  a  digital  computer.  Sections  2.5,  2.6  and  2*7  are  concerned  with  the  applica¬ 
tion  of  tie  theory  to  two  of  McDonnell's  experimental  facilities  and  the  F  region 
of  the  ionosphere.  The  propagation  velocity  and  attenuation  of  Alfve'n  waves 
for  large  ranges  of  parameter  values  are  derived  from  the  computer  plots  and 
then  compared  with  qualitative  curves  for  the  limiting  cases  investigated  by 
Tannenbaum  (Reference  2-3 ) -  In  each  case,  existence  criteria,  as  determined  by 
theoretical  considerations  and  the  computer  results,  are  summarized.  Frequency 
range,  phase  velocity  and  attenuation  properties  of  waves  capable  of  existing 
in  the  three  media  are  examined.  The  possibility  of  wave  channeling  in  Hie 
ionosphere  is  also  investigated. 

2.2  THREE  FLUID  DESCRIPTION.  -  The  basic  plasma  force  equations,  continuity 
equations,  perfect  gas  equations  of  state  and  Maxwell's  equations  can  be  utilized 
to  derive  expressions  for  momentum  transfer  between  the  gas  constituents  and  a 
generalized  Ohm's  law.  From  these  equations  a  basic  wave  equation  can  be  implied, 
and  its  corresponding  dispersion  relationship  obtained.  Following  the  derivation 
of  Lehnert  (Reference  2-1),  a  three-component  ideal  gas  of  Ne  electrons,  ions 
and  Nn  neutrals  per  unit  volume  is  considered.  From  the  basic  plasma  equations, 
the  rate  of  momentum  transfer  between  the  particles  and  between  the  fields  and 
particles  is  expressed  by  the  following  set  of  relations: 

Neme6eve  -  -eNe(E  +  ve  xlT)  -  NeMej  v#i  (ve  -’vj) 

-  Ne  Men  ven  (ve  -  vn)  -  V  Pe  -  pe  V  <t> 

for  the  electron  gas, 

Nj  mj  Vj  ■  e  Nj  (E  +  Vj  X  B)  -  NjMjei/ie  (vj  -  ve) 

-  Nj Mjn vjn  (vj  —  vn)  -  V Pj  —  pj 


for  the  ion  gas,  and: 

Nn  mn  ^n  vn  “  -  NnMnj  un-t  (vn  —  vj)  —  NnMnevne(vn  -  v#) 

-VP 

for  the  neutral  gas. 


(2-1) 

(2-2) 

(2-3) 
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In  these  equations,  Ns  is  the  s-particle  concentration  where  s  is  either  e,  i 
or  n,  and  ms,  its  mass.  The  positive  charges  are  taken  to  he  singly  ionized 
neutral  particles  with  net  charge  e.  The  macroscopic  average  for  s-particle 
velocity  is  denoted  hy  vs,  and  the  notation  Msr  represents  the  reduced  mass 
for  s  and  r-particle  interactions.  The  corresponding  "collision"  frequency 


is  v  .  The  scalar  pressure  and  gravitational  potential  are  denoted  by  Ps  and 
0  ,  while  E  and  B  represent  the  macroscopic  electric  field  and  magnetic  flux 
density  due  to  both  the  plasma  constituents  and  external  sources.  The  operator 
Ds  is  the  hydrodynamic  derivative  for  the  s-particle  gas.  Conservation  of  mass 


and  charge  is  implied  by  the  use  of  the  continuity  equations : 

d  Ns  _ 

~ —  +  V-  (Ns  vs)  *  0;  s  =.  e,  i,  n 


(2-4) 


Spatial  variations  in  the  pressure  will  be  assumed  to  occur  rapidly  enough  for 


the  adiabatic  condition  for  ideal  gases  to  apply.  Then: 


Vps-C  VPS;  s-e,  i,  n 


(2-5) 


where : 


p  -  Nsms 


<r«yskTs/ms 


with  k  being  the  Boltzmann  constant,  7S  the  ratio  fo  specific  heats  and  Ts  the 
temperature  of  the  s-particle  gas.  Under  the  assumption  of  electrical  neutrality 


for  the  whole  gas: 
N.  -  Ni  -  N 


(2-6) 


For  the  manipulation  of  equations  (2-1)  through  (2-3/  it  is  convenient  to 


define  the  following  relations: 


Pn  “  Nn  mn 

p  -  Nj  mj  +  Neme 

m  «  mj  -f 


P  +  P  n 


(2-7) 
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i  -  Ne  (vj  -  ve) 


(2-8) 


ve  »  v  — -  i 

Nem 


v  -  N  (mi  vj  +  meVe)/p 


(2-9) 


P-Pi  +  Pe 


(2-10) 


The  symbol  j  denotes  the  current  density  due  to  the  moving  charges,  irrespective 
of  the  cause  of  their  motion. 

The  three -component  gas  can  now  be  reduced  to  a  two-component  gas,  whose 
constituents  are  referred  to  as  the  plasma  and  neutral  gas.  Substitution  of  the 
relations  (2-7)  through  (2-10)  into  (2-1)  and  (2-2)  yields: 


A  /*•  mj  _A  /-*.  -A  Uij  *  ♦  1  /  mj  \— 

—  NMen  i-'en  (v  —  v,)  —  ^  P*  —  V 


VPe_p  V0 


(2-11) 


A  L.  I"«  -A  L-  —  -A  me  1  /  me  \ 

N’iTtarT,"7“l’'“'  + — M|n  v'i 

-  -  VPj  _  Pi 


/\  /s 


Adding  (2-11)  and  (2-12)  and  expanding  D_  and  gives: 

i  e  1 

A-*-  memj  L.  \  j  —•  — ■  -¥■-*■  . 

pDv+  - 1  j  •  Vj  - -  j  x  B  +  0  j  -  pa(v  -  vm)  -  VP  -  p 


where  the  definitions: 


A  d  -*• 
D-—  +(v.V) 

(7  1 


(2-12) 


(2-13) 


(2-14) 


a“  —  (Men^n +  Min^in) 


(2-15) 


,L 
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expressions  (2-15)  and  (2-l8)  yield: 

-  m*  1  .  . 

ama“  —  ven  +  Kin  (2-23) 

mi  / 

This  result  substantiates  the  interpretation  of  a  ,  since  in  order  to  maintain 
a  non-equilibrium  condition,  momentum  is  not  conserved  unless  the  disturbing 
energy  source  is  included.  Only  then  is  "aV  a , 

The  quantity  (3  takes  on  the  form  of  a  scalar  magnetic  field  in  the  center 
of  mass  frame  of  an  elementary  volume  of  the  gas.  With  the  gas  constituents  at 
thermal  equilibrium,  (2-21)  and  (2-22 )  imply  that: 

_  me 

/3  =  /3=  (2-24) 

e 

Since  (2-ll)  and  (2-12)  are  the  force  equations  for  the  charged  particles'  of 
an  elementary  volume  of  the  cps,  a  generalization  of  Ohm's  law  can  be  obtained 

,  .  mi  .  .  me 

by  multiplying  (2-11)  by  —  and  (2-12)  by  —  and  then  subtracting  the 

ep  ep 

resulting  expressions;  thus: 


^(7) 


mjme  — — ►  mjme(mi-  me)  — ,  .  , 

+  (i -V)v  + - ^ -  (i  -V)(-)  = 


e2  p 


mi  -  me  _ 
E  +  v  x  B - (j 


;T>-y(±) 


p 


+  /S(v  ~  vn)  +  -  (mj  VPe  -  m0VPj) 

ep 


where : 


i  M:, 


T  /me\  /mi\  m  [M  ^  /^in\  V'n  /m'\  ^  /^®n\  "1 

\~J  l'ie  +  V/  ei  +  7  [\7T/  \7~~ )~2  +  \T  /  "enJ 


For  the  condition cf  thermal  equilibrium: 

m 


ime  me 

~  l/en  +  "ei+  2mj  Vin, 


(2-25) 


(2-26) 


(2-27) 


In  a  three  constituent  gas  the  usual  resistivity  term  must  include  the  effects  of 

y 

momentum  transfer  between  the  plasma  and  neutral  gas,  therefore  —  is  the  proper 

P 

first  order  expression  for  the  gas  resistivity.  It  can  be  argued  that,  for  a 
negligible  difference  in  the  accelerations  of  the  ions  and  electrons  and  small 
gradients  and  rotations  of  the  ion  and  electron  velocities,  the  terms  involving 


j  and  the  operator  (7*7)  are  of  no  relative  consequence  in  (2-25). 
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With  these  terms  discarded,  equations  (2-13)  and  (2-25)  can  "be  written  as: 

p  D  v  =  j  x  B  +  /3  j  -  pa  (v  -  vn)  -  VP  -  p  V<£ 


— mj  —  m, 
E  +  v  x  B  m  — 


y  -»■  mi  me 

-  I  x  B  +  —  i  -  /3  (v  -  vn) - (VPe - V  Pj) 

e  p  p  e  p  mj 


(2-28) 

(2-29) 


If  the  particle  masses  are  included  in  (2-11),  the  result  is: 

dPs 

—  +  V  •  p,  vs  =  0;  s-e  n  (2-30) 

The  set  of  equations  (2-17),  (2-28),  (2-29)  and  (2-30)  are  sufficient  to 
describe  the  evolution  of  particle  velocities  and  densities  in  space-time;  and 
hence,  account  for  the  net  current  density.  Maxwell's  equations  then  furnish  the 
connection  between  the  time  varying  currents  and  the  general  space-time  varying 


magnetic  field  produced  in  the  plasma.  The  required  relationships  are: 


-r  <?B 

Vx  E  = - 

d  t 


(2-31) 


and: 


_  T  —  d  E 
V  x  B  =  p0  j  +  pQ  f  — 
at 


(2-32) 

where  /<0  is  the  permeability  of  the  gas,,  essentially  that  of  free  space,  and 
is  the  permittivity,  assumed  here  to  be  a  scalar  quantity.  Except  for  large 
phase  velocities  of  the  disturbances  in  the  plasma,  the  displacement  current  in 
(2-32)  makes  a  negligible  contribution  to  VX  B. 

As  is  customary  in  the  investigation  of  complex  phenomena,  a  first  order 
analysis  of  the  above  set  of  equations  will  noy  be  performed.  Each  quantity  is 
mathematically  decomposed  into  a  constant  component,  which  includes  any  initial 
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value,  and  a  perturbation,  thus: 

B~=  B'o  S' 

E  =  Eo  +  E' 

i  =  i0  +  i  ' 

V  =  Vo  +  v' 

Tn=  Vo+T^ 

Ps  =  Pso  +  p's'  s  =  e'  •'  n 
P  =  Po  +  p' 

where  all  the  zero  subscripted  quantities  are  taken  to.  be  constant,  in  space  and 
time.  Incorporation  of  (2-33)  into  +'ie  equations  expressing  momentum  transport 
(2-28)  and  (2-17),  the  generalized  Ohm's  law  (2-29),  the  equations  of  continuity 
(2-30),  and  the  field  and  current  relationships  from  Maxwell's  equations  (2-31) 
and  (2-32),  produces  a  set  of  equations,  where  all  second  order  terms  of  the 


type 

(jV),  (p'v'),(p'V-v' 
equations  are ; 


4-f) 


and  p'  are  neglected.  Thus  the  linearized 


,D0v'=  j  x  B0  +  /Si  -  p0a(v  -  vn)  -  ^  V p'  +  Cf  Vp.'^ 


(2-34) 


Pno  =  -  ^7+  P0  «  7  S  )  -  ^  Vp  '  (2-35  ) 


+  ni|  —  me  -►  y  -►  -►  -►  mj  n  m®  9 

E  +  v'x  B0  =  -  j x  B0  +  —  i-/3(v'-v')  r -  [C ZVP' - C/vPi'](2-36) 

ep0  P0  ®Po  e  mj  1 


8  P's 


db' 

V  x  E  - - 

dt 


—  —  <9E' 

7  *b'-M  *1*0*57 


(2-37) 

(2-38) 

(2-39) 
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In  (2-34)  and  (2-35)  the  psV<f>  terms  have  also  been  neglected,  since  in 
almost  any  experiment  or  terrestrial  atmospheric  disturbance  they  represent  a 
constant  but  insignificant  part  of  the  total  gas  energy.  The  zero  subscript 
on  the  hydrodynamic  derivatives  means  that  the  variable  velocity  is  replaced  by 
the  constant  velocity  V0. 

Since  all  of  the  quantities  expressed  in  (2-33)  are  macroscopic  averages  and 
the  primed  quantities  are  capable  of  Fourier  analysis,  for  any  physical  case,  si 
normal  mode  analysis  is  justified.  This  procedure  requires  all  the  primed 
variables  in  (2-33)  to  vary  in  proportion  to  _■(<“♦  +  «  •  r)  „  a*.  +->.=> 


,  where  k  is  the  wave 


vector  and  r  is  the  directed  distance  from  the  origin  of  coordinates  to  a  point 
on  the  wave  front.  To  be  physically  applicable,  the  gas  in  which  such  distur¬ 
bances  take  place  must  essentially  be  homogeneous  and  infinite  in  extent.  Upon 
making  the  above  normal  mode  substitution  for  the  dependent  variables,  the  new 
linearized  set  becomes: 

ip0  (co-i-T-  V0)v'=  j  x  B0  +  /3  j  -  p0a  (v '  -v^)  -  i  (C^  p'  +  C^p-)  k  (2-40) 

'Pno^  V°)vn  =~  +  Po“fy/-  vn’)  ~  '  Pn  *  (2-4l) 


E  +  v'xBo--! - -  j  x  B0  +  2— 

eP  o  P  o 

i  m.(c2n'  m e  r2  V 

-  I  - pe - C  Pi  Jk 

ep  \  e  mj  '7 


(2-42) 


Ps  1  ^  _ 

— “  =  -  -ft  (v'  • «  );  s  =  e,  i,  n 
P  so  U  K 


(2-43) 


F'-u^xiT') 


(2-44) 


7=  —  I  -(-T 

W 


(2-45) 
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vith: 


0=^ 

K 


(2-46) 


as  the  complex  phase  velocity  in  the  direction  of  T,  or  simply,  the  phase 
velocity  operator.  The  permittivity  e  has  been  taken  to  be  a  scalar  'whose 
value  is  that  of  the  permittivity  of  free  space,  f0  .  Then  in  the  MKS  system 
the  vacuum  speed  of  light  becomes: 


In  these  equations  is  the  unit  vector  in  the  direction  of  7;  hence:  Jt 

&  ‘K=  ft 

Upon  introduction  of  the  parameters: 


(i 

Pr.  a  \ 


(2-47) 


Po  T 

G  = - F 

Pno  a 


(2-48) 


equation  (2-4l)  may  be  solved  for  v'  and  this  result  used  in  (2-40)  to  relate 
v'to  j  and  the  p,  .  Thus: 

<SJ»> 


In  the  same  manner  equation  (2-42)  becomes: 


(2-50) 


-*•  iFflv'  mi  —  /fi/a  —  y  \  ♦ 

Box'''--; — —  “  B0  x  j  +  I  -  —  0 - )I  +  E 

1  +  iF  e  pQ  \P00  +  'F)  P0 j  U 

+  i'‘[4:J>7Tf  (2-5 

1  '-N 

The  parameter  -  is  interpreted  as  a  coupling  factor  between  the  plasma  and 

F  -0 

Po 

neutral  gases,  since  -  a  is  the  expected  number  of  plasma -neutral  collisions 

^  Pno 

per  unit  time,  and  -  is  the  period  over  which  a  disturbance  ■takes  place.  It 

OJ 
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follows  that  for  -  »  1,  the  neutrals  readily  participate  in  the  motion  of 
F 

the  ions;  while  for  -  «  1,  the  plasma  and  neutral  gas  retain  their  individuality. 
F 

This  interpretation  of  i  is  reasonable  for  nearly  stationary  gases,  where 
|Vo|«|0|  ;  hut  its  meaning  for  the  general  case  is  somewhat  obscure  since 

the  real  part  of  U  is  not,  in  general,  the  phase  velocity  of  the  disturbance. 


Using  the  definitions: 


'.-4- £3 


a/ a  — 

P*  =  P  ~  -T— ^  P 

I  +  iF 


*•-  -£-fM 

PoPo  \  1  +  1 F  / 


(2-51) 


equation  (2-49)  and  (2-50)  can  be  written  as: 


/V0  JV* _ -  /  2  2  2  (2-52) 

+"u~  i  -  Box  i-  iK^Ce  P;  +  C.  p{+ai,Cnp^iK 

and: 

B0  X  v'-  /3*v'-  -  B0X  i  -  g0A*j  +  E  (2-53) 

e  po 

.  f m'  (j2  ,  me  c2  \  Pa*  2  ,1  - 

[ep0  V*  mi  /  “Po  "  J 

The  density  perturbations  are  related  to  the  current  density  and  plasma  velocity 
by  means  of  equation  (2-43),  with  the  substitutions  for  v^  and  f rom  (2-9)>  and 

for  v^'from  (2-41).  The  resulting  expression  for  P„  is: 


A  /  UClpno  \ 

U(l  +  iF)  ll+-^ - 1 

'  U  ( 1  +  i  F)  a  pa  ’ 


~  /3  J  ^ 

v  ~  ~r~=  1  '  ‘k 

L  p°  J 


(2-54) 


The  starred  quantities  of  (2-52)  have  definite  physical  interpretations  for 
nearly  stationary  plasmas  at  thermal  equilibrium.  For  strong  plasma-neutral 
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coupling,  F  «  1,  p *  becomes  the  total  unified  gas  density,  and  /9*  ,  the 
effective  scalar  magnetic  field,  assumes  the  form: 

/3*  -  iF/3  (2-55) 

which  is  smaller  than  /9  by  just  F,  thus  indicating  the  reduction  in  forces 
that  are  perpendicular  to  j  X  B0  on  the  gas  elements .  Accordingly,  the 
resistivity  of  the  gas,  A*  ,  is  reduced  in  the  direction  of  the  elementary 


currents  as  exhibited  by: 


y  /34 

pQx*m - — 

P0  aPo 


(2-56) 


For  small  percentages  of  ionization: 


uin'  "ei  <<ven 


(2-57) 


and  therefore: 


which  is  essentially: 


At  the  other  extreme,  that  of  weak  coupling,  F»1  implies  that  is 


(2-58) 


(2-59) 


essentially  p0  ,  the  plasma  density,  and  /8*  is  just  /6  .  The  resistivity  term 
is  the  usual  scalar  resistivity  for  "the  predominant  electron  contribution,  v/p0  f 


provided  F  » 


Q£) 


,  which  is  generally  oflhe  order  of  unity.  In 


this  case,  <x  ,  the  scalar  conductivity,  it  is  properly  expressed  as: 

kj  .2 


(2-60) 


The  neutral  gas  acoustic  coupling  is  determined  by: 


and  hence,  for  weak  coupling,  a*  tends  to  zero.  This  implies  that  for  weak 
coupling,  acoustic  distr-’bances  of  the  plasma  do  not  involve  the  neutrals, 
except  perhaps,  as  a  dissipation  mechanism.  The  role  of  in  the  theory  is 
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quite  important  since  the  magnetostriciive  pressure,  which  is  responsible  for 
transverse  plasma  oscillations,  is  effective  when  |B0|  is  significantly  larger 
than  |/3*|  .  This  discussion  further  illustrates  the  meaning  of  F,  in  that  for 
F «  1  the  gas  motion  is  unified,  while  for  F  »  1,  the  component  gases  act  nearly 
independent  of  one  another. 

To  facilitate  manipulation  of  the  combination  of  equations  (2-33)>  (2-44), 
(2-45),  (2-52),  (2-53)  and  (2-54)  the  dyadic  formalism  is  employed,  and  the 
nonlan  form  of  the  operators  and  Tg  x  is  utilized,  where: 


‘Bo"  |Bo|  ' 

After  using  (2-44)  and  (2-45)  to  reduce  E  in  (2-53),  the  momentum  equation 
(2-52)  and  the  generalized  Ohm’s  law  (2-53),  can  be  written  as: 

(a  1 1  +  a2"£TK )  -T'  +  (a3 1  +  a4XtK  +  °5*^b0  x )  *T»  (°6 1  +  a7^Bo  x )  •  i«  ( 2-6l ) 


(eiT+  c2  *«  + c3tb0  x : >  • ' v'  +  (c4  i 1  +  c5  1k, v +  c6rBo  x )  -  r 
-  (cp~+  c8ttk + c4rBo  x)  -  T0  +  E0 

where  I  is  the  unit  dyadic.  The  coefficients  aj  and  are  as  follows: 


(2-62) 


’l"* 

jx  (r2  _2  Ji 

j2-y^  ^+ciPio+’»cn 


°-j  ■  -  p* 


-ix  "*i  /  2  Oh  r2  \  P  2  1 

u  iC  ,C"'WJ 


(2-63) 


ac"  B0 


k  f3  2 

B6"^*+i  TTiy^n^o 


°7*  °5 
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cr°3 


c2-a4 


e3-°5 


e4“^o 


i  K 

C.m  7T 

5  U 


K  + 


A 

U 


FI 


(2-64) 


c6-“ 


•Po 


o 


where: 


c7“-c4 


C8  — 
aa  p 


c9  “  "  C6 


9CnPno 


V 


(T  +  iF)  11  +  i  r 


-2 

*CnPno 


aPoLl(l  +  iF )} 


Equations  (2-6l)  and  (2-62)  can  be  more  compactly  expressed  by  the  metric  equation 


A11  a12 


A21  A22 


(2-65) 


O 


with  the  A^  as  the  dyadic  coefficients  and  R]_  and  Rg  as  the  constant  vectors  on 
the  right  side  of  (2-6l)  and  and  (2-62),  respectively.  It  can  be  shown  that  the 
inverse  of  the  square  matrix  exists,  an  therefore,  (2-65)  may  be  solved  for  y/ 

j 

and  J?  0 
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2.3  THE  BASIC  DISPERSION  RELATIONSHIP.  -  Prom  the  homogenous  form  of 


equation  (2-65),  that  is. 


Rj  m  R2  ■  0 


(2-66) 


a  general  dispersion  relationship  can  he  found.  Equation  (2-66)  implies  that 


no  constant  electric  field  or  current  is  maintained  in  the  fluid  by  external 


energy  sources;  hence s 


Eo-io-° 


(2-67) 


With  this  restr,J  ,j/on,  and  the  definition  of  T  (equations  2-33  with  2-9) , 


(2-65)  can  be  written  as: 


A11  a12 
A21  A22 


A11  A12 
A21  A22 


(2-68) 


Since  the  column  matrices  are  independent,  the  left  hand  terms  of  (2-68)  must 
vanish  separately,  in  order  for  the  metric  equation  to  be  valid.  Tanenbaum  and 
Mintzer  (Reference  2-1)  have  analyzed  only  the  first  term.  However,  from  the 
standpoint  of  developing  a  dispersion  equation,  either  term  may  be  used.  If 
(2-65)  Is  expanded  and  solved  for  v',  the  following  result  is  obtained,  when 
the  inverse  of  A22  exists : 

(A1 1  “  A12  ’  A22 ^ "  A21  ^  * *'m  ©  (2-69) 

Once  a  base-system  of  vectors  Is  chosen  (e.g.  the  unit  vectors  of  a  rectangular 
Cartesian  coordinate  system),  equation  (2-69)  is  found  to  be  satisfied  only  if 
certain  relationships  exist  between  the  aj  and  cj^  of  (2-63)  and  (2-64).  These 
relationships,  taken  together,  specify  a  dispersion  equation  for  wave  propagation 
in  an  infinite  gas  of  any  degree  of  ionization.  This  equation  will  involve  u  and 
k ,  either  or  both  of  which  may  be  complex. 


The  dispersion  equation  may  be  regarded  as  a  polynominal  in  either  u  or  «, 


whichever  is  appropriate  to  the  physical  situation.  If  a  spatial  disturbance 
of  wavelength: 
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A-  (2-70  f 

K 

is  maintained  throughout  the  measuring  of  propagation  characteristics  such  as 
attenuation  and  propagation  velocity,  k  is  real  and  specified,  and  therefore,  the 
polynomial  should  "be  solved  for  complex  u  .  The  roots  will  give  the  eigen- 
frequencies  and  associated  attenuation  time,  defined  as  the  time  for  the  wave 
amplitude  to  reach  l/e  of  its  initial  value.  Since  the  time  dependent  variables 


are  taken  as  proportional  to  e'“  ,  the  attenuation  of  time  r  becomes 


(2-71) 


If  instead  of  wavelength  the  angular  frequency  is  specified,  the  dispersion 
polynomial  should  be  solved  for  complex  k.  The  roots  will  then  yield  the  proper 
wavelength  for  a  given  w  ;  and  since  the  space  variables  go  as  •  '*  * " ,  the 


attenuation  distance  can  be  obtained  from: 


|lm  (k)| 


(2-72) 


Should  the  wave  be  excited  by  an  energy  pulse  localized  in  space,  either  w  or  k 
may  be  specified  as  complex  quantities  and  the  dispersion  relationship  solved  for 
the  other.  Attenuation  in  both  time  and  space  will  occur.  ( 

The  propagation  velocity  will  be  the  phase  ■velocity  for  the  first  two  cases 


just  discussed.  It  Is  always  expressed  as: 

R  .(*>) 

Uph  “  ReU) 

For  small  spatial  damping,  that  is  Im(*)«Re(K): 


(2-73) 


Uph-Re(U)  (2-74) 

In  the  situation  where  a  localized  pulse  initiates  the  disturbance ,  the  propagation 
velocity  of  the  signal  is  the  group  velocity  Ugj.  defined  by: 

•-SS 

Each  frequency  component  or  wavelength  will  have  its  own  phase  velocity,  as 
given  by  (2-73). 
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2.4  SPECIALIZATION  OF  BASIC  DISPERSION  EQUATION.  -  It  is  obvious  that  the 
dispersion  equation  will  contain  more  information  than  necessary  to  construct  a 
reasonably  complete  picture  of  the  types  of  waves  and  their  propagation 
characteristics.  It  is  also  apparent  that  such  an  expression  would  be  retract¬ 
able  even  for  a  qualitative  description  of  the  phenomena.  Therefore,  it  is 
desirable  to  examine  specific  cases  which,  when  taken  as  a  whole, yield  a  good 
qualitative  description  of  permissible  disturbances.  The  approach  taken  by  various 
investigators  is  to  discuss  longitudinal  and  transverse  wave  propagation,  under 
thermal  equilibrium  of  a  stationary  gas,  perpendicular  to  and  parallel  to  a  : 

permeating  constant  magnetic  field,  whose  flux  density  is  BQ.  Such  was  the 
approach  selected  in  the  investigations.  The  conditions  for  which  the  remaining 
discussions  of  section  2  are  valid  are  enumerated  in  the  following  statements: 

Thermal  Equilibrium:  Te  -  Tj  -  Tn;  a  -~a,  /9  -  7? 

F  Pnow  a) 

Stationary  Gas:  V0  -  0;  —  - - ,  G  «  — 

a  Po  a  “  _ 

A  rectangular  coordinate  system  is  used,  and  the  tz  axis  is  taken  parallel  to 


With  the  above  substitutions  the  factors  in  (2-63)  and  (2-64)  become: 

°i  “  •  <yp* 

°2“  "D"  (2c2p0  +  T'CnP"0) 

°3m-P* 
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and: 


where : 
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c  1  "  °3 


c2“°4 


c3  -  °5 


=4  “Me 


A 

u 


‘5‘0 


c6‘  Q 


^c2p°t(^)2,c"p~ 


(1  +  iF) 


1  +  iF  U  7T- 


Cz- 


y  k  T 


Q=  N. 


'Po 


(2-77) 


(2-78) 


These  expressions  will  now  he  used  to  find  the  explicit  dispersion  relationships 
for  the  four  cases  mentioned  above. 

2.4.1  Wave  Propagation  Perpendicular  to  Magnetic  Field.  -  The  relation 
between  the  propagation  yi'ctor  and  magnetic  field  vector  in  this  case  is  such 
that  Ik  is  perpendicular  to  lb0  .  For  convenience  Ik  will  be  taken  parallel  to 
the  x-axis.  Equation  (2-69)  now  reduces  to: 


W11  w12  0  T 

W21  w22  ° 

0  0  w33! 


(2-79) 
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where: 

W11  “  <°1  +  “2^  +  "T  *e3<c3c4  -  2c’e6>  -  [2e2(=3c6)  +  (c1  +  c2)2  e4+  e^cjll 


w12-  -  —  ((ci  +  c2^(c3e4  +  cjc^)  +  c^c3e6  -  c,  (c4  +  Cj)]) 


w21  -  -  —  {(c1  +  c2)(c3c4  -  cjcg)  +  cj|[c3c6  +  cj  (c4+  cj)]} 


(2-80) 


"22“  “1  +  T  ,c3(c3c4  *  clc6)  *  cl[c3c6+  C1  (c4+  C5W1 


o 


w33  “  °1  “  — 


and: 


c2-  c4  <c4  +  c5)+  c6 

Equation  (2-79)  implies  the  existence  of  a  transverse  wave,  linearly  polarized  in 
the  direction  of  the  magnetic  field.  Hence: 

w33  “  0  (2-81) 
is  required,  since,  v^may  have  any  value  that  is  consistent  with  the  linearization 
used  in  deriving  (2-65).  For  moderately  strong  plasma-neutral  coupling  (2-81) 


can  he  expressed  as  follows,  when  the  appropriate  substitutions  are  aade: 


a 

2  2 
k* c* 


1  -■ 


& 


2  2 
K  C 


K) 


(2-fie) 


where: 


•■t  t-m-) 
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and 


K)fe)’(9 


This  iq  essentially,  a  result  obtained  in  magneto-ionic  theory  for  electromagnetic 
waves  which  are  altered  by  the  charged  particle  response  to  the  radiation  fields. 
Such  disturbances  can  propagate  without  attenuation  when  v  =  0  and  w>n  ;  but 
the  attenuation  is  a  fairly  complicated  function  of  u  if  v  is  allowed  arbitrary 
values.  For  v  tending  to  infinity,  the  phase  velocity  operator  («■>/«)  is  real, 
and  the  wave  propagation  is  non-dispersive  and  unattenuated.  The  phase  velocity 

o 

is  therefore  equal  to  the  speed  of  light  in  a  vacuum. 

The  system  of  equations  in  the  vx'and  vy' components  of  (2-79)  imply  a  coupling 
of  longitudinal  and  transverse  disturbances.  Physically,  this  is  due  to  the 
charged  particles  which  sire  moving  in  the  direction  of  propagation  being  subjected 
to  a  Lorentz  force  from  the  presence  of  B^.  The  dispersion  equation  is: 

wllw22  ~  W21  w12“  0  (2-83) 

When  expanded,  this  expression  is  a  complicated  fourth  order  equation  in  k.  It 
has  been  amply  discussed  in  Reference  2-1  for  very  high  and  very  low  wave 
frequencies,  showing  that  at  low  frequencies  coupled  acoustic  and  magneto¬ 
hydrodynamic  waves  involving  the  entire  gas  occur.  As  the  frequency  is  increased, 
several  critical  regions  are  passed,  with  the  oscillations  eventually  going  from 
plasma  acoustic  waves  to  ion  acoustic  waves.  A  qualitative  description  of  the 
complicated  wave  behavior  for  propagation  across  the  magnetic  field  is  presented 
in  Figure  2-1,  which  is  reproduced  from  that  article. 

2.4.2  Longitudinal  Wave  Propagation  Parallel  to  Magnetic  Field.  -  In  this 
instance  (2-69)  can  be  expanded  to: 


WU  w12  0 

< 

V 

-w12  W11  0 

- 

0 

0  0  w33_ 

0 
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with: 

WM  “  °1  “  -T  [cl(clc4  +  C3C6)  ‘  c3(c3c4  “  CVC6)] 


w12  “  "T  tc3(clc4+  C3C6)  +  cl(c3e4  *  C1C6W 

Cl 


(2-85) 


and: 


w33 


<M  +  °2  “ 


(c'i  +  c2)2 
c4+c5 


c 


2 

o 


+  c 


2 

6 


Again,  an  appropriate  solution  exists  if: 

W33  +  w^2)  -  0  (2-86) 

Since  the  vz  component  for  a  longitudinal  wave  cannot  he  identically  zero, 
the  requirement  for  wave  propagation  ’becomes: 

w33  “  0  (2-87) 

As  pointed  out  by  Reference  2-1,  examination  of  the  dispersion  equation  (2-87) 
reveals  that  at  low  frequencies  the  constituent  gases  unify,  and  an  acoustic 
wave  with  velocity  characteristic  of  the  entire  gas  is  propagated.  When 
<u  >  Max  [vin,  vni]  the  coupling  factor  (l/F)  becomes  quite  small  for  low  degrees  of 
ionization.  Hence,  two  acoustic  waves,  one  associated  with  the  plasma  and 
another  with  the  neutral  gas,  are  disseminated.  As  «  becomes  even  larger,  such 

1 

that  <u~<uci,  the  plasma  breaks  up.  At  this  point  an  ion  acoustic  wave  and  an 
electron  acoustic  wave  propagate  with  their  separate  phase  velocities  going  to  the 
appropriate  acoustic  velocities.  In  the  range  o>ci  <0  <&>ce  ,  the  wave  velocity 
mathematically  may  exceed  c  by  a  considerable  amount.  This  may  be  interpreted 
as  a  "stiffening"  of  the  electron  gas,  or  a  velocity  resonance  condition.  When 
the  phase  velocity  of  the  electron  disturbance  is  greater  than  c  or  the 
wavelengths  become  so  short  that  the  continuum  assumptions  in  the 
basic  equations  are  invalid,  the  wave  cannot  be  regarded  as  a  physical  entity. 
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Such  conditions  exist  in  part  of  the  range  cocl<co  <<uce  .  '' 

2.4.3  Transverse  Wave  Propagation  Parallel  to  Magnetic  Field.  -  For  the 
over-all  objectives  of  the  program  this  case  is  the  most  relevant.  As  will  be 
shown  in  the  ensuing  discussions,  a  singular ily  important  type  of  wave,  called 
an  Alfven  wave,  emerges.  This  wave  usually  prevails  for  a  wide  range  of 
disturbance  frequencies.  At  the  higher  frequencies  a>  >  <uci ,  the  transverse  waves 
of  the  magneto-ionic  theory  come  into  existence. 

The  proper  dispersion  equation  is  again  derived  from  equation  (2-69).  Now, 
however,  v'  must  be  identically  zero,  while  vx' and  vy' cannot  both  vanish.  Thus, 
from  (2-69):  * 


W11  w  12 1  vx 


12  W11  vy 


(2-88) 


Therefore: 


is  required.  This  implies  that: 


2  2  n 

wll  +  w12-° 


Wn»i  IWi 


If  this  result  is  then  used  in  (2-88)  it  is  found  that: 


(2-89) 


(2-90) 


(2-91) 


or  that  the  transverse  waves  are  circularly  polarized.  When  viewed  along  the 
normal  to  the  plane  of  rotation,  the  plus  sign  is  seen  to  properly  represent  a 
left  circularly  polarized  wave  (LCP),  and  the  minus  sign,  a  right  circularly 


polarized  wave  (RCP) .  Expanding  (2-90)  yields: 

2  2 

°1  (c4  +  c6^  “  C(cic4  +  c3c6)  t  i  (c3c4  -  c1c6)](c1  ±  ic3)  =  0 

After  substituting  the  values  from  (2-76)  and  (2-77)  and  manipulating  the 
resulting  expression,  while  neglecting  (u/c)^  as  compared  to  unity,  the 


(2-92) 
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dispersion  relationship  becomes: 


AO  CJ 

(U2  ±  — 


«&-«?(•■  3^-' 


(2-93) 


If  the  gas  is  not  stationary  but  moving  with  velocity  Vo  along  <z,  equation 
(2-93)  is  modified  to 

(a2‘^)('4M'44'“44)-°  <**> 


In  these  equations: 


(2-95) 


These  quantities  are  the  modified  cyclotron  frequency,  the  Alfven  velocity  for 
a  fully  ionized  gas  and  the  Alfve'n  velocity  operator  for  a  partially  ionized  gas. 
It  is  the  reciprocal  of  «.  ,  as  it  appears  in  (2-93),  that  has  physical  meaning 


1  mi-  me  1 


which  is  the  difference  in  the  time  required  for  an  ion  and  an  electron  to  subtend 

Uo 

one  radian  of  arc  about  the  Bc  field  "lines".  In  a  fully  ionized  medium  —  is 

Qj°  1 

Just  the  distance  a  constant  phase  point  of  the  Alfven  wave  moves  in  the  time  —  . 

Poa  . 

At  low  frequencies  <u  « - ,  and  therefore,  1/F  » 1,  indicating  strong 

Pno 

coupling.  At  these  frequencies  P*  is  essentially  the  total  gas  density  (p0  +  pno)  • 
Thus,  the  Alfven  velocity  operator  is  real  and  becomes  the  Alfven  velocity  for 


the  whole  gas,  that  is: 


(2-96) 


When  cj£— — ,  the  coupling  is  weak  and  U*  is  complex,  then  UQ  is  the  Alfven 
Pno 
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velocity  of  the  plasma,  and  the  role  of  the  neutral  gas  is  to  dissipate  the  wave 
energy  through  the  plasma -neutral  collisions.  It  can  he  shown  (Reference  2-2) 
that  for 

(2-97) 

the  dispersion  equation  (2-93)  becomes  that  of  magneto-ionic  theory  for  wave 
propagation  along  a  magnetic  field  (Reference  2-4). 

Equation  (2-93)  is  the  dispersion  relationship  which  determines  the 
existence  or  non-existence  of  an  Alfven  type  plane  wave  under  the  conditions  of 
total  gas  neutrality,  negligible  gravitational  gradients  and  turbulence,  scalar 
gas  pressure  and  small  transverse  perturbations  of  a  stationary  gas.  The  positive 
and  negative  signs  imply  that  there  are  different  dispersion  equations  for  Jeft 

(~''j 

and  right  circularly  polarized  plane  waves. 

2. 4. 3-1  The  Wave  Equation  for  Alfven  Waves  in  a  Partially  Ionized  Gas.  -  An 
examination  of  the  wave  equation  to  which  (2-93)  applies  is  useful  for  several 
reasons .  A  better  understanding  can  be  obtained  of  what  combination  of  basic 
parameters  (collision  frequencies,  cyclotron  frequencies,  densities,  etc.)  are 
effective  in  defining  the  spatial  characteristics  of  a  wave.  The  use  of  the 
classical  form  of  a  wave  equation  allows  a  somewhat  easier  analysis  in  situations 
involving  boundary  conditions .  In  addition  Alfven  wave  existence  criteria  may 
be  readily  formulated  from  the  wave  equation.  To  derive  this  equation  in  such 
a  way  that  it  explicitly  displays  the  perturbed  component  of  the  magnetic  flux 
density,  a  directly  measurable  variable,  equation  (2-68)  is  solved  for  jj  and 
(2-4p)  is  used  to  get  an  equation  similar  to  (2-68)  but  in  terms  of  b' .  Com¬ 
bining  the  restrictions  in  Section  2.4,  which  were  used  for  the  specialized 
dispersion  relationship,  with  the  form  of  (2-88)  to  eliminate  the  pressure  terms, 
equations  (2-52)  and  (2-53)  with  (2-44)  and  (2-45)  can  be  combined  to  give: 
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- 

x* 

(Vx  b ')  xTz 

xTz 

+ 

w2B  2  " 

i 

6) 

VxVx  b' 

Jr  0 

20*U2 

u2 

0 

h  ~  ' 

V  x 

(Vx"b')  > 

**Z 

- 

0 

m2B0 

0>a>o 

\  I 


(2-98) 


2 

where  terms  in  (u/c)  and  (V0/u)  have  been  neglected.  The  complex  operator 
i k  has  been  replaced  by  the  differential  operator  V  in  this  expression.  If 
(2-98)  is  multiplied  by  (-<u2)  and  the  operator  substitution  ieo  =  d/dt 
is  used,  equation  (2-98)  is  found  to  be  the  desired  general  Alfven  wave  equation. 


For  plane  waves  in  an  unbounded  medium,  equation  (2-93)  is  the  corresponding 
dispersion  relationship.  If  an  experiment  is  performed  in  a  cylindrical  chamber, 
where  the  geometry  is  such  that  the  radial  dimension  is  sufficiently  less  than 
the  length,  and  the  constant  magnetic  field  B0  is  along  the  cylinder  axis,  the 
apparatus  is  referred  to  as  a  hydromagnetic  (HM)  wave  guide  (Reference  2-5). 

Under  such  circumstances,  cylindrical  coordinates  are  used  in  (2-98),  and  a 
dispersion  relationship  which  exhibits  discrete  frequency  modes  is  acquired  from 
the  usual  conditions  on  a  solution  to  the  differential  equation.  Since  nearly 
fully  ionized  stationary  gases  can  be  produced  in  an  HM  wave  guide,  the  assumption 
of  a  complete  plasma  medium  is  warranted.  In  the  next  section,  a  discussion  of 
the  wave  equation  for  a  fully  ionized  gas  and  a  development  of  the  wave  equation 
for  wave  guides  are  presented. 

2.4.3. 2  The  Wave  Equation  for  Alfven  Waves  in  a  Fully  Ionized  Gas ♦  -  Appendix 
A  is  a  derivation  of  the  wave  equation  for  a  fully  ionized  gas  using  the  original 
approach  of  Alfven;  however,  it  is  instructive  to  examine  equation  (2-98)  in 
this  case  and  compare  the  result  with  that  found  in  the  Appendix.  By  dividing 
(2-98)  by  the  coefficient  of  its  second  term  and  defining: 
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(r 


U? 


C= 


\*CJ 


ft 

^  =  2  —  i 

B0 


U2<u 

o 


(2-99) 


equation  (2-98)  can  be  written  as: 

£b'+ V  x  j  £z  x[iz  x  V  x  b'] !  +  £tVx  (Vx  b')]  -  i/V  x  [  tz  x  (Vx  b')l  =  0  (2-100) 

After  repeated  use  of  certain  vector  identities  and  rearranging  terms,  equation 
(2-100)  further  reduces  to: 

V2(£b'-b'i  )  -yfrfax  ~  )  +  czx(vx  ~ )  -fb'>=0  (2-101) 

dz  dz 

For  complete  ionization  of  a  stationary  gas,  the  following  conditions  exist: 

The  neutral  density  pn  is  zero,  and  consequently. 

a  =  /3=0 


and: 


Since  F  remains'  finite: 


mi  me*'ej 


(2-102) 


/3*=  0 


Wo 


Therefore : 


U*=U„ 


'V 

0 


(2-103) 


<T= 


GJ 

<P  =  ~  i  - 


(2-10U) 
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where  the  approximation  m  «  mi  has  been  used.  After  substituting  (2-104), 
equation  (2-101)  becomes: 

i<uvei  o  _  _  db"  co  db'  co2  *  .acn 

-  V2b'+y2b't , -T ,x(Vx - )  -  i  - Vx —  +  —  b'=0  (2-105) 

“ce  ^ci  ^ z  “ci  ^ z  \J2 

o 

Alfven's  result,  equation  (A-29),  can  be  written  in  the  form: 

*  ,2?%  (2-106) 

— -  V2b  +  — —  +  — b'=0 

,.2  -,,2  ..2 

U  u  a  oz  U 

o  r°  o 

The  assumption  that  <u«<acj  is  at  least  reasonable,  if  not  necessary,  when 
discussing  Alfven  waves;  hence,  the  fourth  term  in  (2-105)  may  be  neglected. 

Upon  substituting  for  U0  and  recalling  from  (2-60)  the  expression  for  the 
conductivity,  the  coefficient  of  V2  in  (2-106)  is  seen  to  be  just  that  of  the 
corresponding  term  in  (2-105).  By  expanding  the  third  term  in  (2-105)  it  is 
found  that:  ,  —  .  ;>T' 


(2-106) 


_  /  <9b'\  d2b' 


(2-10J) 


9  ■».  Cbl 

Hence  the  terms  Vz b'  cz  andy — -that  are  present  in  (2-105),  are  missing 

dz 

from  (2-106).  For  transverse  wave  disturbances,  however,  these  two  terms 
vanish,  and  Either  equation  is  valid.  Expressing  y2  in  cylindrical  coordinates 
with  jz  along  the  wave  guide  axis,  equation  (2-106)  in  matric  form  becomes 


A  1 


0  k; 


A  1 

n  —  ob; 


(2-108) 


where  the  operator  n  is  given  by: 


n  b; 


2  2 

i  co  2  d  co 

tZ*Vr'e'Z+d'2  +Uo 


(2-109) 
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The  electric  field  E  can  "be  expressed  by  (2-32)  with  j  =  a  E.  For  the 
longitudinal  components,  Gould  (Reference  2-6)  writes: 


*  CsnJn(rsnr)e 


i  (cot  -  K  z  -  rtO) 


(2-110) 


E:-  *  D,„Jn'A.n'>' 
s,  n 


i(d)t  -  K  z  -  nd) 


(2-111) 


where  J„  is  the  n-th  order  Bessel  function,  K  is  the  longitudinal  wave 
n  sn 


number  and  fsn  is  related  to  A_n  by: 


•(U  /c)2  ~  ]/£ 


with: 


A2  =r2  0 

sn  sn  (U./c)2  +  1/0-0 


r9  <3-  1  k2 

Sn=  c2  U2(l-0  " 


(2-112) 


(2-113) 


These  are  determined  by  the  boundary  condition  imposed  on  either  b'z  or  E'z. 

The  coefficients  CSn  and  D__  are  calculated  by  utilizing  the  orthogonality 
property  of  the  Bessel  functions  and  the  physical  specifications  of  the  wave 
exciter.  The  transverse  b*  and  E'  components,  which  cause  the  torsional  mode, 
may  be  found  by  first  substituting  (2-110)  and  (2-111)  into  (2-108),  and  making 
use  of  (2-112).  Rather  than  using  this  method,  the  equations  for  b'z  and  E'z 
are  used  for  the  discussion  of  transverse  modes  for  n  =  0,  which  is  of  primary 
importance  to  investigations  in  the  EM  wave  guide.  For  Alfven  waves  the  dis¬ 
placement  current  is  negligible  compared  to  the  conduction  current.  Using 
(2-113)  in  (2-112)  and  discarding  the  terms  U0^/c^  and  the  terms  w2/c2,  which  are 
absent  in  (2-98),  the  dispersion  equation  applicable  to  the  propagation  of  Alfven 
waves  in  a  wave  guide  is  found  to  be: 

a2, + (i-i/z)  k20=o  (2-ni0 

for  the  circularly  symmetric  modes  (n  =  0). 

The  solutions  of  transverse  wave  dispersion  equations  for  an  unbounded 
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nfidium  and  a  hydromagnetic  wave  guide  were  considered  very  important  to  the 
achievement  of  the  over-all  program  objectives.  It  was,  therefore,  desirable  to 
investigate  these  equations  for  a  wide  variety  of  plasma  conditions  in  order  to 
obtain  qualitative  and  quantitative  information  as  to  the  existence  and  propagation 
characteristics  of  Alfv/n  waves  in  both  the  laboratory  apparatus  and  the  ionosphere. 


2.4.4  Parametric  Analysis  of  Transverse  Wave  Dispersion  Equations.  -  A 
parametric  study  is  valuable  in  analysis  when  the  situation  under  investigation 
is  very  complicated,  or  when  the  number  and  range  of  parameter  values  to  be  con¬ 
sidered  is  great.  Both  situations  were  found  to  exist  when  physical  inter¬ 
pretations  of  the  dispersion  relationships  were  attempted  for  the  many  possible 
experimental  and  ionospheric  conditions.  For  these  reasons  most  investigators 
were  forced  to  consider  only  limiting  cases  which,  while  important  from  a  theo¬ 
retical  standpoint,  were  not  always  applicable  to  experimental  conditions.  To 
fill  in  the  missing  knowledge  for  non-limiting  cases,  and  to  compile  quantitative 
information  on  transverse  wave  propagation  characteristics  for  a  large  range  of 
parameters  and  variable  values,  a  program  for  a  high  speed  digital  computer  was 
written.  A  summary  of  the  capabilities  of  this  program  is  presented  in  Table  2-1.' 
The  computed  curves,  which  influenced  the  choice  of  wave  exciters  and  detectors 
used  in  the  experimental  program,  are  discussed  in  section  2.5  for  the  HIT 
experiments  and  section  2.7  for  an  ionospheric  model.  The  remainder  of  this 
section  is  devoted  to  the  reformulation  of  equation  (2-94)  that  was  found  suitable 
for  the  computer  investigations. 

Since  the  dispersion  relations  may  be  considered  a6  polynomials  in  either 
co  or  k  ,  they  must  be  formulated  in  two  different  ways  for  computation. 

A 

Cage  1 

If  the  value  of  k  is  specified,  the  solutions  of  the  resulting  polynomial 
in  oo  represent  a  set  of  characteristics  frequencies  which  will  propagate  with  the 
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TABLE  2-1 


COMPUTER  PROGRAM  SUMMARY 


INPUT 

EQUA- 

(IN  ORDER 

DESCRIPTION 

OUTPUT 

TION 

VALIDITY 

OF  USE) 

NO. 

MASS  RATIO 

mi/me 

SIGN 

SIGN  OF  WAVE  POLARIZATION 

S  =  1 !  LCP  WAVE  DATA 

2-93 

Ks  =  0 

S  =  -Is  RCP  WAVE  DATA 

2-93 

K$  =  0 

KS 

LONGITUDINAL  WAVE  NUMBER  (REAL) 

Kj  =  0:  FOUR  COMPLEX 
cos  ROOTS 

2-112 

INFINITE  PLANE  OR 

(2-120) 

TORSIONAL  WAVES  IN  A 
STATIONARY  GAS  (V  =  0) 
FOR  WHICH  fi,«  B  °AND 

«, «  "0- 

K-  =  0:  co  ROOTS 

2-115 

INFINITE  PLANE  WAVES, 

TEST  =  0 

K-  =  0:  k  ROOTS 

2-125 

INFINITE  PLANE  WAVESi 

TEST  =  0 

V 

GAS  FLOW  VELOCITY 

V 

— 

INFINITE  PLANE  WAVES, 

(Kj  =  0) 

B 

AXIAL  MAGNETIC  FLUX  DENSITY 

B 

— 

0 

o 

Pn 

NEUTRAL  GAS  DENSITY 

P  n 

2-7 

,  P 

PLASMA  DENSITY 

P 

2-7 

V.i 

ELECTRON-ION  COLLISION  FREQUENCY 

B-37 

V 

ELECTRON-NEUTRAL  COLLISION 

V _ 

B-37 

FREQUENCY 

en 

Vin 

ION-NEUTRAL  COLLISION  FREQUENCY 

♦'in 

<,'in/v.n> 

B-37 

"o 

2-95 

a 

2-15 

P 

2-24 

y 

2-27 

R 

2-7 

Uo 

2-95 

OPTIONS: 

ia  /v 

in  «n 

v,/p 

3 

-S' 

3 

R 

IONIZATION  RATIO 

✓  J 

Cl) 

WAVE  ANGULAR  FREQUENCY  (REAL) 

TWO  COMPLEX  k  ROOTS 

2-125 

INFINITE  PLANE  WAVES 
(Ks  -  0), 

IN  A  STATIONARY  GAS 

(vo  =  °) 

K 

WAVE  NUMBER  (REAL) 

FOUR  COMPLEX  a)  ROOTS 

2-115 

INFINITE  PLANE  WAVES 

(K,  -  0) 

TEST 

DISPERSION  EQUATION  FORMULATION 

TEST  =  0:  co  ROOTS 

2-115 

TEST  =  0:  k  ROOTS 

2-125 
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specified,  wave  number .  In  general  the  characteristic  frequencies  are  complex, 
the  real  part  giving  the  real  frequency  in  radians  per  second  and  the  imaginary 
part  giving  the  rate  of  decay  of  the  wave  in  nepers  per  second. 

For  plane  waves  the  general  expression  given  in  u xuation  (2-9*0  may  be 
used.  This  takes  the  following  form,  after  substituting  (2-95),  (2-51)  and 
(2-16),  as  restricted  in  section  2.1,  and  taking  Pn  and  P  to  have  essentially 
their  undisturbed  values  Pno  and  P0  : 


where : 


(Aj  +  i  A2)  +  (A3  +  i  A4)  x3  +  (As  +  i  A^)  x2 

+  (A7  +  i  Ag)  X  +  (A9  +  i  A10)  =  0 


x  =  1  CO 


O 

Al=Po"or  P 


a2  =  o 


r\ 

A3  =  co0  r  [p0ap{2  +  r)  +  rx^y] 


r 2  P*Vo  "  Sl<2  Bo 


A5  =  [(1  +  0  (p0a^p  -  r  x2  p2)  +  (2  +  r)  (r  x2  ay) 

-^2r2(3p0PV02-  B^)]  +  3sB02x3r2  VQ 

A6  =  r  [(2  +  r)  a  (2pQ  coQ  p  VQ  x  -  sx2  B2) 

+  oi0  x2  r  (2s  p  B0  +  3k  Vq  y)] 

A7  =  k2  (<w0[( T  +  r)  a  (ay-  /32)  -  (2+  r) p0  rap  V2  +  2B2  r a 
-  r2  V0  (4s  x/3  B0  +  3  K2  Vo  yM  +  2s  B2  rxa  VQ  (2  +  r)i 
Ag  -  (1  +  0(popa2  x  V0  <u0  -  2  co0  r  p2  VG  k3  -  sa2  x2  B2) 

+  2(2  +  r)  a0  ax3  r  V  y  +  rx2  (2r  x  V  &>„  B2  +  3*r  x2  V2  B2 


-  Po  P  r  K  Vo  +  2s  P  Bo 

A9  =  x2  0,0  [( 1  +  r)  r  «2  p2  V2  -  (2  +  r)  x2  V2  ayr  +  B2  (a2  -  r2  V2  x2) 

-  2srx  V0  p  B0  a]  +  (1  +  r)  a2  X  s  B2  V0  -  s  r2  x3  V3  B2j 


(2-115) 


(2-116) 
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A10  =  * 2  V0  t<u0  [(1  +  r)  m  (ay-  £2)  -  2s/3  BQ  r2x2  V0 
-  y  r2  x3  V2+  2a  r  x  B2  +  (2+  r)  sar/t2  B2V0I 
r  =  Po^P 

S  =  t  1 

for  a  stationary  medium  (Vo  =  0)  the  coefficients  (2-116)  reduce 
A1  =  Po  “o  1,2  P 
A2  =*  0 

A3  *  <i)0  r  [^0ap  (2  +  r)  4-  r  y] 

A 4  -  -  sr  k2  B2 

Aj  «  0)0  [(1  +  r)  (fia  a2  p  -  r  x2  02)  +  (2  +  r)  (r  x2  ay)  +  k2  r2  B2] 

A6  =  sr  «2B0[26j0  fir  -  a( 2  +  r)  BQ] 

Ay  =  a  x2  <uQ  1(1  +  r)(ay-  /S2)  +  2r  B2] 

Ag  =  sax2  BQ  [2w0  r  £  -  a  ( 1  +  r)  B„] 

A 9  =  a2  x2  <u0  B2 

Aio  =  ° 

while  for  frequencies  well  below  the  cyclotron  frequency  <dc  ,  the 

become:  .  9 

Ai  ■k.O 

a2  =  o 

A3  =  r  [p0  ap  (2  +  r)  +  r  x^  y] 

A4  '  %0P  '2  *  V0 

A5  »  (1  +  r)(pQ  a2  p  -  r  x2  £2)  +  (2  +  r)  r  x2  ay 
-x2  r2(3p0pV2-B2) 

A6  ”  r  t2PoP“  Vo  x  (2  +  r)  +  X2  r  (2s/3  B„  +  3x  V0  y)] 

Ay  -  X2  [(1  +  r)  a  (ay-  /S2)  -  p0  pr  (2+  r)  a  V2  +  2B2  r  a 
-  r2  V0  (4sx/3  B0  +  3x2  V0  y)] 

As  ■  x  V0(l  +  r)(pQ  pa2  -  2r  /S2x2)  +  2ar  (2+  r)  x3  V0  y 
+  rx2  <2  r,<  Vo  Bo  -  P0P  «  V2  +  2$  $  B0  a) 
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coefficients 


(2-118) 
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A9-K2[r(l+r)  x2/S2V2-  r(2+r)K2V2ay 
+  B2  (a2  -  r2  V2  k2)  —  2  six  VQ  /9  B0  a] 

A10  “  V0  [,<  (1  +  r)  a.  (ay-l S2)  -  2  $  /3  B0  r2  k2  Ve 
-  y  r2  x3  V2  +  2  a  rx  B2] 


(2-118) 

continued 


If  both  conditions  are  satisfied,  that  is,  To  >  0  tod  <o«tu0,  ,  the  coefficients 
are  further  simplified  to: 


A1  -  f*0  r  P 


a2-o 


A3  -  r  [p0ap  (2  +  r)  +  r  xz  y) 


A4-0 


A5  -  (1  +  r)(p0a2p  .  rx2  8s)  +  r  (2  +  r)  x2  ay  +  x2  r2  B2 


(2-119) 


A6.2sr2K2/3B0 


A7-aK2[(T+r)(ay-82)+2rB2) 
Ag-  2sax2r0  B0 
A 9-a2x2  B2 


A,o-0 


For  torsional  plane  waves  in  a  stationary  gas,  equation  (2-112),  with  the 
definitions  (2-104),  (2-95)  and  (2-60),  yields  a  polynomial  in  <0  of  the  form: 

<a4  —  2i  p0  <7  c2  <u3  —  c2  (p2  <r2  c2+  A2  +  K2  )<u2  (2-120) 

+  i  p0  <7  c4  (A^  +  K2  )<a  +  p2  cr2  c4  U2  K2  -0 
where  the  subscript  zero  has  been  dropped  and  U0  yc  neglected. 


By  a  change  of  variable  from  uto  x>  i<u  ,  this  equation  may  be  put  into  the 


more  convenient  form: 


x4  +  2  p0  <7  c2  x3  +  c2  (p2 a2  c2+A2  +K2)x2 
+  to  <7  c4  (A2  +  K 2  )  x  +  p2  a2  e4  U2  K2  -  0 


(2-121) 
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having  all  real  coefficients. 

Case  2 

If  the  value  of  <o  is  specified,  the  resulting  polynomials  in  *  may  he  solved 
for  a  set  of  characteristic  wave  numbers.  These  also  are  generally  complex  valued, 
the  real  part  being  the  wave  number,  in  radians  per  meter,  and  the  imaginary  part 
giving  the  spatial  rate  of  decay  in  nepers  per  meter. 


Rearranging  equation  (2-120)  into  a  polynomial  in  K,_ft  yields  the  form: 

SO 


[c2  (a,2  -\ii0a  c2o)-  K2  -  [J  -  2  i  p0  a  c2  a> 3 


(2-122) 


-  C2  (jl2  CT2  C2  +  A2  )  CD2  +  i  Hoa  c4  Aj  to]  «  0 
which  is  valid  for  torsional  waves,  but  restricted  to  the  special  case  of 
a  stationary  medium  and  a  frequency  >&>  much  less  than  the  cyclotron  frequency 

For  plane  waves  in  an  unbounded  medium,  a  rearrangement  of  equation  (2-94) 
with  the  appropriate  substitutions  gives: 

( Am  ±  i  Am)  ts  ^  t  A _ L.  i  A  .)  i  ( A _ I  i  A  \ 


(Aj  +  i  A2)  K3  +  (A3  +  i  A4)  k*  +  (Aj  +  i  Ag) 

+  (Ay  +  i  Ag)  k2  +  (A^  +  i  A^q)  k  +  (Ai  1  +  i  A]2>  "  0 


(2-123) 


where: 


A]  -  sr2  V2  B2 
A2-r2V2o)0y 

A3  -  3  sr2  V2  o>  B2  -  r  (1  +  r)  V2  <u0  /S2+  r  (2  +  r)  V2  co0  a  y  +  r2  V2  a>0  Bj' 

A4  -  3r2  V2  «0  6>y  +  2 sr2  V2  coQ  fi  B,  -  sr  (2+  r)  V2  a  B2 

Aj  -  3  sr2  VQ  <u2  B2  -  2r  (1  +  r)  V0  coco0  /32  +  2r  (2  +  r)  VQ  oxo0  ay  (2-124) 

+  2r2  VQ  axo 0  B2  -  r2  V2  qm0  n0p+2sr  V0  ooQ  a/9  B0 
-»0  +  r)Vo  a2  B2 

-  3  r 2  Vc  co2  co0  y  +  4  sr2  V2  <u0  /3  BQ  -  2  sr  (2  +  r)  V0  a  a  B2 
-(Ur)Ve<Uoa(ay-/32)-2rVe<i,0aB2 
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Ay-  sr2  &>3  B2  -  r  (1  +  r)  to2  coQ  ) S2  +  r  (2+  r)  to2  <u0  ay 

-  3  r2  V2  a)2  a0  pQ  p  +  r2  o>2  o)0  B2  -  s  (1  +  r)  <u  a2  B2 
2  2 

+  2  sr  a><u0  «/3  30  —  a>0  aL  B£ 

Ag  -  r2  0?  a>0  y  -  sr  (2  +  r)  a>2  a  B2  +  2sr2  <u2  a>0  /3  BQ 


—  (1  +  r)  coa>0a  (ay  -  jS2)  -  2r  oa)Q  a  B2  (2-124) 

continued 

+  r  (2+  r)  VQ  (oa)0  p0  pa 
A9«  (1  +  t)  VQ  iaa>Q  p0  pa2  -  3  r2V0<u3a>0p0p 
2 

A iq  «  2r  (2+  r)  V0  or  a>0  pQpa 

0  a  04 

An  =  0  +  r)  coZ0)o  p0paz  -  rL  aT  b>0  p0p 
A12  =  r  (2+  r)  coJ  co0  p0  pa 

Again  this  may  be  simplified  for  the  various  special  cases. 

(a)  For  a  stationary  medium,  A3,  A2,  A^,  A^,  A5,  Ag,  A^  and  A10  reduce 


to  zero,  leaving: 

(Ay  +  i  Ag)  k2+  (All  +  '  A12^  “  ® 


(2-125) 


with: 


Ay*  sr2o>3  B2  -  r  (1  +  r)  a >2<u0  /S2+  r  (2+  r)  <u20)o  ay 
+  r2  w2  <o0  B2  -  s  (1  +  r)  <oa2  B2  +  2  sr  <o0  <0  a  /3  B0 

-*>o  “2  Bo 

Ag  *  r2  w3  a>0  y  —  sr  (2  +  r)  tu2  a  B2  +  2  sr*-  a)2  <u0  ^  BQ 
-  ( 1  +  r)  axo0  a  (ay  -  £S2)  -  2  r a)0>0 a  B2 
A^  -  ( 1  +  r)  w2  cu0  p0  pa2  -  r2  c 0*  co0  p0  p 
A12“  r  (2+  r)  or  <o0  pQ  pa 


(2-126) 


(b)  For  a)  «  tu0 , 


A|  *  0 

A2-r2  V3  y 


Ag  =  r  V2  [(2  +  r)  a  y  -  (1  +  r) /S2  +  r  B2] 
A4-r2V2(3Wy+2sjSB0) 


(2-127) 
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Aj  -  r  V0  [2(2  +  r)  coay  -  2  (1  +  r)  «  /32 
+  2ro)  B2  -  r  V2  cup0p+  2s  a£  BQ  ] 

A6  »  VQ  [3r2  <u2y  +  4  sr2  VQ  /3  BQ  -  (1  +  r)  a  (ay  -  /92)  -  2r  a  B2] 
Ay  =  r  (2  +  r)  <y2a  y  -  r  (1  +  r)  ca2  /32  -  3r^  V2  <u2  p0  p 
+  r2  a1  tu0  B2  +  2  sr  coa  (3  BQ  -  a2  B2 
Ag  -  co  [r2  <u2y  +  2  sr2  a>  /3  BQ  -  ( 1  +  r)  a  (a  y-  /32) 

-  2 r  a  B2  +  r  (2  +  r)  V2  pQ  pa] 

A9  =  (1  +  r)  VQ  <up0pa2  -  3r2  VQ  cu3  pQp 

A]0  ”  2r  (2+  r)  v0  "2^0Pa 

An  =  (1  +  r)  <u^pQpa/  -  >  a)*V0 p 

O 

A^  =  r(2  +  r)  pQ  pa 


(2-127) 

continued 


(c)  For  a  stationary  medium  and  «<<«0  all  coefficients  vanish  except 
Ay,  Ag,  A^  and  A]_2  leaving: 


A7  =  r  (2  +  r)  eu2  a  y  -  r  ( 1  +  r)  o>2  j82  +  r2aj2  B2 
+  2  sr  co2  af3  BQ  -  a2B2 

Ag  ■  r2  cu3  y  +  2  sr2  w2  / 3  BQ  -  (1  +  r)  ma  (ay  -  |9  ) 
-  2r  axz  B2 

A]]  =  (1  +  r)  or  p0  pa^  -  r*  cu*  p 
A]2  =  f  (2  +  r)  co0  n0pa 


(2-128) 


Three  of  the  above  expanded  equations,  (2-112),  (2-115)  and  (2-125),  were  put 
into  one  computer  program  with  the  selection  of  the  one  appropriate  equation 
being  made  part  of  the  input  data.  Only  the  collision  frequencies  are 
required  to  be  computed  outside  of  the  program. 


Up  to  this  point  considerable  effort  has  been  spent  establishing  a 
reasonably  firm  background  in  the  theory  of  MHD  waves,  with  the  bulk  of  the 
analysis  concerned  with  Alfven  waves .  The  sections  that  follow  apply  part  of  the 
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theory  directly  and  part  of  it  indirectly  through  the  results  of  the  computer 
studies.  Application  is  made  to  the  three  pertinent  physical  situations: 
experiments  in  McDonnell's  hypervelocity  impulse  tunnel,  experiments  in 
McDonnell's  arc  discharge  tube,  and  disturbances  in  the  F-region  of  the  iono¬ 
sphere.  In  each  case  criteria  for  the  existence  of  Alfven  waves  and  their 
propagation  characteristics  are  discussed. 

2.5  ALFVEN  WAVES  IN  THE  HYFERVELOCITY  IMPULSE  TUNNEL.  -  Section  3. 4. 2. 6 
establishes  the  plasma  environment  for  the  experiments.  The  various  means  of 
wave  excitation  and  the  actual  experiments  performed  are  delineated  in  Sections 
4.1  and  7.1  of  Vol.  III.  To  utilize  the  computer  results  for  an  unbounded  medium, 
the  experiment  instrumentation  and  tunnel  conditions  must  be  such  that  the 
density,  temperature,  ionization  ratio  and  magnetic  field  are  essentially  constant 
over  a  characteristic  distance  X0.  This  distance  is  established  by  the 
requirement  that : 


X0>Uph  *0 

(2-129) 

O 

II 

*C|o- 

+ 

3 

(2-130) 

with  ZQ  being  the  distance  from  the  point  of  excitation  to  the  detector  and  tjj, 
the  time  during  which  measurements  are  taken.  Assuming  this  to  be  the  experimental 
situation  a  set  of  "most  probable"  conditions  was  taken  as  a  reference.  Quantita¬ 
tive  investigations  in  the  six  dimensional  space  of  parameters  were  made  about 
this  set  by  varying  one  parameter  at  a  time.  These  conditions  and  the  associated 
primary  natural  frequencies  are  given  in  Table  2-2. 

2.5.1  Existence  of  Alfven  Waves  in  the  Tunnel.  -  Two  approaches  are  taken 
here;  one  involves  examination  of  the  computer  curves  for  wave  existence,  and 
the  other  utilizes  a  fairly  obvious  general  criterion  suggested  by  Watanabe 
(Reference  2-2).  For  the  computer  analysis  approach  only  a  limited  study  has 

been  performed  using  equation  (2-115)  for  <0  as  a  function  of  *,  since  the  case 
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HIT  REFERENCE  CONDITIONS 


PARAMETER 

DESCRIPTION 

VALUE 

mi/ma 

MASS  OF  N+  TO  MASS  OF  ELECTRON 

2.55  (104) 

Pn 

NEUTRAL  PARTICLE  DENSITY  CORRESPONDING 

TO  10-3  ATMOSPHERES  PRESSURE 

1.72  (10-4)kg/m3 

R 

IONIZATION  RATIO,  R  -  P 

P  +  Pn 

10-3 

T 

EQUILIBRIUM  TEMPERATURE 

103  °K 

V 

GAS  FLOW  VELOCITY 

V  «U  , 

o  ph 

10_1  WEBER/m2 

B 

o 

AXIAL  CONSTANT  MAGNETIC  FLUX  DENSITY 

V.n 

ELECTRON-NEUTRAL  COLLISION  FR EQUENC Y  (A-37) 

1.14  (107)/SEC. 

"in 

ION-NEUTRAL  COLLISION  FREQUENCY  (A-37) 

4.04(105)/SEC. 

".i 

ELECTRON-ION  COLLISION  FREQUENCY  (A-37) 

1.03(106)/SEC. 

“«• 

ELECTRON  CYCLOTRON  FREQUENCY 

1.67(109)RAD./SEC. 

wo 

«B 

a,  =  - —  =o >  .,  ION  CYCLOTRON  FREQUENCY 

0  _  m  Cl' 

mi  -  m« 

6.90  ( 1 05)  RAD. /SEC. 

for  which  equation  (2-125)  is  valid  more  accurately  represents  the  conditions 
which  existed  during  the  HIT  program.  The  case  of  interest  is  that  for  which 
a  disturbance  frequency  is  maintained  and  the  spatial  progress  is  to  be  measured. 
The  computer  program  for  determining  the  k  roots  is  not  capable  of  handling 
moving  plasmas,  which  involve  a  fifth  order  polynomial  in  k  (2-123) j  therefore, 
the  results  portrayed  in  Figures  2-2  through  2-lk  are  for  stationary  plasmas. 
Examination  of  these  figures  indicates  that  existence  is  determined  by  the  amount 
of  information  necessary  to  at  least  ascertain  the  periodicity  of  the  disturbance. 
Although  a  fairly  reliable  prediction  might  be  obtained  by  extrapolation  of  a 
detected  partial  cycle,  it  is  intuitively  more  desirable  to  detect  a  full  cycle 
of  the  disturbance  before  claiming  existence.  This  somewhat  arbitrary  criterion 
can  be  expressed  as: 

5d  > A  (2-13l) 
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EFFECT  OF  MOLECULAR  WEIGHT  ON 
ATTENUATION  DISTANCE 
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Another  method  of  deriving  existence  criteria  is  that  used  by  Watanabe  in 
Reference  2-7*  He  uses  criteria  which  essentially  reduce  equation  (2-98)  to 
the  form  of  (2-106),  without  the  first  term,  thereby  transforming  it  to  the 
original  equation  discovered  by  Alfven  for  totally  ionized  gases.  Weakly  ionized 


gases  act  the  same  in  most  respects  as  fully  ionized  gases,  if  the  coupling 
parameter,  l/F,  is  small.  The  criteria  so  developed  are  somewhat  stringent, 
since  they' are  a  result  of  requiring  the  waves  to  be  completely  free  of  frictional 
damping.  However,  if  these  conditions  are  satisfied,  the  probability  of 
generating  a  wave  is  greatly  enhanced.  Since  in  the  HIT  B0  =0.1  weber/m^  and 
B  =  6.3  (10~5)  weber/m^,  that  is,  the  vector  field  is  much  stronger  than  the 
average  of  the  random  scalar  field,  the  existence  requirements  in  the  weakly 
ionized  gas  of  the  HIT  become  equivalent  to  those  for  a  strong  field  as  defined 
by  Watanabe.  If  in  addition  the  plasma -neutral  coupling  is  small,  l/F  is  much 
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less  than  unity,  requiring: 

*>»  J~a  (2-135) 

r  n 

This  restriction  is  augmented  by  the  fact  that  in  order  to  appreciably  reduce  the 
drag  of  the  neutrals  on  the  plasma: 

co»  a  (2-136) 

As  was  pointed  out  at  the  end  of  Section  2. 4. 3-1,  equation  (2-98)  can  be  shown 
to  agree  with  (2-106),  for  plane  waves  in  a  fully  ionized  gas.  But  in  going 
from  (2-98)  to  (2-105)  and  then  to  (2-106)  a  damping  term: 

i  to  d  b , 

—  Vx  — 

ct>cj  az 

was  neglected.  This  is  justified  provided: 

GXC^i  (2-137) 

The  remaining  damping  term  in  (2-105)  can  be  reduced  to  frequency  factors  by 
making  use  of  the  relations: 

1 


and: 


Then: 


where : 


F  »  1 


I  0) 


I  Ct>  V 


V  =  V  •  +  V'  +  “ — ~  Vi¬ 
ol  on  2m-  n 


(2-138) 

(2-139) 


Hence  the  additional  requirement: 

“ce  ^ci 

co  « - 

v 

must  be  imposed.  Thus,  the  result  of  taking  (2-136),  (2-13?)  and  (2-139) 
together  implies  that  Alfven  waves  can  be  generated  if  co  is  restricted  as 
follows : 

...  T  "cij 

a  « co  .«  Mill  — - - 1  (2-l40) 
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Inserting  the  reference  conditions  for  the  HIT  requires: 

4.05  ( 105)  <<  co  .«  5. 90  ( 105)  rodi  ans/sec.  (2-141) 

The  frequency  range  is  seen  to  occur  in  the  region  of  relatively  small 
attenuation  distances  for  the  curves  of  Figure  2-2.  But  the  inequality  (2-l4l) 
cannot  he  met  under  these  conditions,  and  therefore  the  attenuation  in  this 
range  need  not  he  small. 

For  strong  coupling,  l/F«l,  implying  that: 

«  —  (2-142) 

pn  2 

CO 

Thus,  the  results  (2-55)  and  (2-5 6)  are  valid.  Upon  multiplying  (2-98)  by  — — 


the  reduction  of  damping  to  a  negligible  amountrequires  that: 

o 

2/3  F  Uo  ^ 

<<J 


(2-143) 


(2-144) 


where  (2-55)  has  been  used  for  /3*  .  Since  for  strong  coupling  and  strong  field 
conditions /9«  B0  and  F«l,  relation  (2-143)  becomes: 

ujp  .  (2-145) 


By  making  the  substitutions  (2-95)/  this  reduces  to: 


Using  (2-56)  and  (2-95)>  condition  (2-144)  becomes: 

0-  Cl) 

— ~(y-p2/a)«l 

P*l 

From  the  definitions  of  a  ,  /S  and  y  it  can  be  shown  that: 

/32 

—  «Y 

a 

and  therefore,  (2-14?)  yields  the  requirement: 

p  "ce  wci 

CO  «  - 

Pn  v 

where  v  is  given  by  (2-138).  Thus  the  condition  for  wave  existence  is: 


P 

0  <  co  «  —  Min  a,  ct>cj , 

Pn  L 


(2-146) 


(2-147) 


(2-148) 


(2-149) 
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The  HIT  reference  parameters  then  give  the  following  range  for  a>  s 

0  <cu  «0.2  radian/ sec.  (2-150) 

Since  the  total  test  time  was  only  about  5  milliseconds,  this  criterion 
could  not  be  used  to  establish  the  exciter  frequency. 

2.5.2  Alfven  Wave  Characteristics  for  Tunnel  Conditions.  -  The 

characteristics  of  interest  are  the  frequency  range,  wave  lengths,  phase 

velocities  and  attenuation  distances  for  Alfven  waves  that  can  be  generated  in 

the  HIT.  Figure  2-2  shows  the  attenuation  distance  as  a  function  of  the  real 

part  of  the  wave  number  for  the  reference  gas,  illustrating  the  difference  in 

propagation  characteristics  of  right  (plus  sign)  and  right  (minus  sign)  circularly 

polarized  waves.  The  curve  for  the  right  circularly  polarized  wave  (rcf)  exhibits 

a  relative  minimum  attenuation  distance  at  Re  [*]«  6. 31  x  10-2  radians/meter. 

Both  curves  change  direction  at  high  frequencies,  and  eventually,  each  achieves 

a  slope  of  minus  one, this  having  been  verified  analytically.  The  second  change 

always  occurs  when  a>  is  greater  than  either  a  or  ,  and  consequently,  the 

coupling  factor  is  much  less  than  —  ,  implying  that  no  appreciable  coupling 

Pn 

exists  between  the  plasma  and  neutral  particles.  Furthermore,  such  frequencies 
are  well  above  the  electron-ion  and  ion-electron  collision  frequencies,  and 
the  electrons,  therefore,  tend  to  move  independent  of  the  much  more  massive  ions. 
Because  of  the  high  frequencies  involved,  this  occurs  beyond  the  region  where 
Alfven  waves  should  exist.  The  wave  lengths  for  the  RCP  curve  soon  become  so 
short,  that  the  curve  has  no  physical  realization.  On  the  other  hand,  the  ICP 
waves,  while  exhibiting  the  long  wave  lengths  associated  with  magneto-ionic 
disturbances,  are  so  rapidly  attenuated  that  they  also  loose  their  physical 
significance.  Figure  2-3  relates  phase  velocity  to  Alfven  velocity  for  the  RCP 
waves .  There  is  an  inflection  point  near  &>0  where  the  wave  begins  to  change  its 
identity.  This  occurs  far  below  the  ion  plasma  frequency  <upi  which  is  defined 
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in  the  usual  way: 

M  -V2 

pl  \fo  ml  / 

For  the  LCP,  Up^  exceeds  c  at  a  point  near  <o  =  2(10^)  radians  per  second. 

The  result  of  varying  the  parameters  about  the  reference  conditions  is 
shown  in  Figures  2-4  through  2-13-  They  portray  the  influence  of  the  parameters 
R  (ionization  ratio),  pn  ,.  B0,  T  and  mj_.  Only  the  curves  for  the  RCP  wave  are 
given  since  the  LCP  wave  reaches  a  critical  point  at  a>0  ,  and  thereafter,  its 
behavior  becomes  extremely  complex. 

Figure  2-4,  for  a  partially  ionized  gas,  may  bo  compared  with  Figure  A-l, 
for  a  fully  ionized  gas.  It  is  observed  that  the  region  of  increasing  attenua¬ 
tion  distance,  d,  does  not  occur  for  a  fully  ionized  gas.  The  relationship 

between  r0,  and  R  in  the  figures  is: 

106  R 

I*  <*  - 

°  1  -  R 

Hence  R  -  10" 2  corresponds,  approximately,  to  rD  =  10^.  The  partially  ionized 
and  fully  ionized  propagation  characteristics  are  thus  seen  to  merge  quite 
smoothly.  The  region  of  increasing  d  occurs  when  the  neutral  particles  pre¬ 
dominate.  For  the  conditions  considered  in  Figure  2-4,  the  neutral  gas  is  by 
far  the  major  component;  it  thereby  becomes  a  mechanism  for  energy  exchange  with 
the  plasma.  A  much  smaller  attenuation  distance  results  for  a  given  o  in 
the  partially  ionized  case .  In  the  region  of  increasing  attenuation  distance 
less  attenuation  may  be  exhibited  than  for  the  fully  ionized  medium  since  not 
only  is  the  frictional  force  reduced  in  this  region,  but  also  the  Lorentz  force 
which  gives  rise  to  the  transverse  oscillations.  The  increase  in  d  beyond  the 
minimum  of  a  or  <o0  then  represents  a  transition  from  Alfven  waves  to  the  waves 
of  magneto-ionic  theory.  Here  the  phase  velocity  is  always  greater  than  the 
Alfven  velocity  Uq,  while  in  other  portions  of  the  curves,  the  phase  velocity 
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may  have  almost  any  value,  depending  on  the  degree  of  ionization.  For  large 
enough  a  ,  Up^  will  always  exceed  c. 

Figure  2-5  summarizes  the  effect  of  co  on  attenuation  distance  for  various 
ionization  ratios.  In  this  figure  the  curve  for  a  given  ionization  (R)  crosses 
only  curves  of  higher  charge  concentration.  Hence,  for  all  values  of  co  below  an 
intersection,  the  wave  corresponding  to  higher  ionization  has  the  least  attenua¬ 
tion.  The  tendency  for  the  neutral  particles  to  support  and  then  separate  from 
the  gas  disturbance  is  clearly  exhibited  in  the  range  co  <  a . 

The  resulting  of  holding  R  constant  at  10  ^  and  varying  pn  is  portrayed  in 
Figure  2-6.  For  the  rarefied  gases  of  the  upper  curves  the  attenuation  distance 

is  nearly  constant  over  several  orders  of  magnitude  of  co  .  The  inequality  (2-148) 

-6 

is  approximately  satisfied  for  the  Pn  =  10  curve .  As  Pn  increases ,  the  neutrals 

take  more  and  more  energy  from  the  plasma  motion  at  low  co  ,  and  the  flattened  regions 

indicating  unified  gas  motion  diminish  until  at  the  curve  of  greatest  Pn  the 

Aifven  wave  is  practically  non-existent.  Figures  2-7  and  2-8  show  that  the  effect 

of  increasing  the  stationary  magnetic  field  B0  is  similar  to  decreasing  Pn  - 

Hence,  combinations  of  Pn  and  Bo  in  an  experiment  should  be  able  to  satisfy 

most  requirements  for  long  and  constant  attenuation  distances  over  a  broad 

range  of  vave  lengths  or  frequencies. 

The  position  of  the  relative  attenuation  minimum  depends  on  the  value  of 

a  =  a (T)  .  Thus  in  the  curves  of  Figure  2-9  the  minima  shift  to  the  right  as 

the  temperature  is  increased.  For  fixed  co  greater  than  a  and  fixed  wave  number, 

Re(*),  there  is  an  attenuation  distance  associated  with  each  value  of  T..  In 

this  region  the  phase  velocity  of  the  wave  is  independent  of  temperature,  since 

✓ 

multiple  particle,  and  therefore,  Aifven  wave  field,  oscillations  occur  before 
collisions  with  neutrals  can  introduce  particle  drag. 

In  Figure  2-10,  showing  the  effect  of  molecular  weight,  curves  are  shown 
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for  singly  ionized  helium,  monatomic  nitrogen,  diatomic  nitrogen,  and  argon.  <> 

A  gas  pressure  of  10" 3  atmospheres  has  been  chosen,  so  that  the  gas  density 

(pn  +  p)  is  proportional  to  the  molecular  weight. 

9 

Over  the  frequency  range  0  <  a  <  10  radian/second,  the  phase  velocity 
for  a  fully  ionized  gas  is  nearly  constant  and  equal  to  the  Alfven  velocity  U0. 

Figure  2-11  shows  the  phase  velocity  as  a  function  of  angular  frequency  for  the 
standard  condition  R  *  10“^  and  the  fully  ionized  condition  R  =  1,  with 
p  =«  I.72  x(l0-10) kg/m^.  The  curve  for  R  •  10“^ can  be  represented  by  a  simple 
exponential  equation. 

The  ratio  of  phase  velocity  to  Alfven  velocity  for  changing  u  is 
presented  in  Figure  2-12  with  R  as  a  parameter.  Here  it  is  seen  that  for  inter¬ 
mediate  frequencies  <0  <  10^  radians/second,  the  phase  velocity  for  waves  in  a 
partially  ionized  gas  is  much  less  than  the  phase  velocity  for  a  fully  ionized  gas, 
indicating  the  effect  of  the  neutral  gas  inertia.  At  the  attenuation  minima 
displayed  in  Figure  2-4  the  phase  velocities  approach  the  velocity  computed  for 
a  fully  ionized  gas  having  the  reference  neutral  density.  The  phase  velocities  have 
a  tendency  to  attain  the  same  value  in  the  frequency  range  10^<  <u<10^  radians  per 

second.  All  curves  for  R  <  10  ^  fall  nearly  on  top  of  the  R  =  10“6  curve. 

Figure  2-13  depicts  the  relationship  of  to  ^  for  the  HIT  reference 

conditions .  The  quantity  is  referred  to  as  the  "damping  factor"  by 

U  '  . 

Tanenbaum  and  Mintzer  in  Reference  2-2.  Since  2ir-^  is  the  wave  length  for  an 
Alfven  wave  in  a  fully  ionized  gas  and  d  is  the  attenuation  distance  in  a 
partially  ionized  gas,  the  damping  factor  provides  a  conservative  lower  limit 
for  8  in  (2-131).  Thus,  the  detector  sensitivity  requirement  varies  appreciably 
with  o)  as  portrayed  in  the  figure . 

2.6  ALFVEN  WAVES  IN  THE  ARC  DISCHARGE  TUBE.  -  A  description  of  the 
conditions  prevalent  in  the  arc  tube  is  given  in  Section  3*4. 3.2.  This  facility 
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is  an  example  of  an  HM  wave  guide  and  is  similar  to  the  one  used  by  Wilcox  et  al 

(Reference  2-8).  A  discussion  of  the  existence  and  propagation  of  waves  in 

A  this  device  must  proceed  from  the  solutions  (2-110)  and  (2-111)  of  the  longitudinal 

components  of  "the  wave  equations  for  b  and  E.  Since  the  allowed  values  of  wave 

number  Ks  are  discrete,  each  mode  s  (corresponding  to  the  s-th  zero  of  the  Bessel 
function  J0^has  its  own  propagation  characteristics.  Therefore  both  existence 
conditions  and  propagation  features  can  be  determined  from  a  mode  analysis. 

2.6.1  Existence  of  Alfven  Waves  in  the  Arc  Tube.  -  Three  modes,  or  types, 
of  propagation  can  be  considered  -  the  transverse  magnetic  mode  (TM),  the  trans¬ 
verse  electric  mode  (TE),  and  the  transverse  electromagnetic  mode  (TEM).  The 
TM  modes  are  derivable  from  E^and  are  called  E  modes,  and  TE  modes  are  derivable 
from  b/  and  are  called  B  modes.  The  TEM  modes  are  characterized  by  both  E/ and 
tz  vanishing. 


REPORT  NO.  A219 
30  NOVEMBER  1963 


M  H  D  Wave  Investigation  3ornpovetmSer 

II  -  theory 


ally,  that  j_(R,z)  =  0  for  a  similar  facility.  Using  equation  (2-151),  this 


leads  to  the  boundary  condition: 


-<l(AsoR)»0 


(2-158) 


On  the  other  hand,  Gould  (Reference  2-6)  puts  e£  (R,z)  =  0  and  finds  that  the 


condition: 


Jo(AsoR)-0  (2-159) 
is  appropriate.  For  the  TM  mode  of  excitation  used  in  the  arc  tube,  E/  is  given 
by  the  expression  for  a  concentric  electrode  capacitor,  namely: 


r  In  - 
R' 


,  R'  <  r  <  R 


(2-I60) 


,0<r<R' 


where  VQ  is  the  initial  voltage  and  R  and  R  are  the  radial  distances  of 
the  inner  and  outer  electrodes  from  the  tube  axis.  The  coefficients  DgQ  in 
equations  (2-151)  and  (2-152)  can  be  determined  from  the  orthogonality  property 
of  the  Jn  and  the  above  relation  for  E/.  Substituting  (2-l60)  into  (2-151)  and 
integrating  over  r  from  the  tube  axis  to  outer  electrode: 


C  ~~R  J1  <A*o'>  rdr  -  f  DS0  M  -  -7 -)/R  J1  (Aso  '>  J1  (Ato'>  * 

r  In  —  ^  so 

R' 

Then  using  (2-158)  and  the  recurrence  relation  for  Bessel  Functions: 


(2-161) 


j2(AtoR>-Jo(At0R) 


the  integration  yields: 


H)^ 


R  lnfe) 


Jo  (AgR)  -  J0  (A^RQ 

j o  (AS0R) 


(2-162) 


Then  substituting  into  (2-152)  the  value  of  Dso  obtained  from  (2-162)  gives  the 
magnetic  flux  density  for  the  torsional  mode  which  was  experimentally  investi¬ 
gated  in  the  McDonnell  arc  tube.  The  derived  result  is: 
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J0(ASOR)-J0(A$0R')  ...  .  ifet-M 

- - -  .  J,  (AMr)e 

kma$0j0  (AsoR) 


(2-163) 


2.6.2  Alfven  Wave  Characteristics  in  the  Arc  Tube.  -  For  the  case  of  a 

strong  axial  magnetic  field  B0  and  infinite  conductivity,  the  field  lines  are 

fixed  in  the  plasma,  and  the  Alfven  waves  travel  along  them  unattenuated.  In  the 

case  of  TM  modes  there  is,  in  addition,  no  dispersion.  The  magnetic  field  for 

the  TM  modes  is  tangential,  and  hence,  a  torsional  wave  propagates.  For  TE  modes 

there  exists  a  radial  component  of  t>,/  causing  E^to  encompass  the  b"0  lines.  The 

attenuation  distance  -jr — j- — rt  is  obtained  from  the  roots  of  dispersion  equation 

I lm  (Ksom 

( 2-114) .  In  the  case  of  small  damping  of  the  wave  and  an  Alfven  velocity  UQ  much 
less  than  the  velocity  of  light  c,  equation  (2-114)  can  be  manipulated  to  give: 


2^o  aUo 


(2-164) 


The  value  of  E&  along  the  tube  axis  is  just  the  z  -  component  of  the  gradient  \ 
of  the  potential  V(r,0,z,t),  hence: 

E'«  — 
z  <9  z 

For  given  r,  9,  and  t,  this  electric  field  is  proportional  to  the  voltage 
associated  with  the  attenuated  wave,  therefore: 

(2-165) 

Utilizing  the  roots  of  (2-114),  the  phase  velocity  along  the  tube  axis 
can  be  found  from: 


z“  dz  "  ~  d 

u*o 


u 


Re  (cos) 

ph‘  (Kto) 


(2-166) 


/  _ 

2.7  ALFVEN  WAVES  IN  THE  IONOSPHERE.  -  The  ionospheric  region  of  interest 
and  its  pertinent  characteristics  are  defined  in  Section  3*2.  An  altitude  of 
300  km  has  been  chosen  for  reference  conditions.  Two  models  can  be  examined. 
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an  unbounded  medium  and  an  array  of  wave  channels.  A  wave  channel  model  seems 
most  reasonable,  particularily  since  there  are  two  wave  restraining  phenomena, 
the  magnetic  pressure  and  the  ionized  gas  pressure.  The  existence  criteria  are 
the  same  as  those  hypothesized  in  paragraph  2.5*1  for  both  models,  but  the 
propagation  characteristics  must,  be  considered  separately.  The  conditions  for 
applicability  of  the  infinite  medium  model  are  outlined  at  the  beginning  of 
Section  2.5*  While  a  discussion  of  propagation  in  an  infinite  medium  is 
straightforward,  the  discussion  for  a  wave  channel  is  considerably  more  complex, 
and  only  a  cursory  examination  is  presented  here. 

2.7*1  Existence  of  Alfven  Waves  inlhe  Ionosphere.  -  Applying  the  criterion 
expressed  by  (2-133)  to  the  curves  of  Figure  2-15,  it  is  found  that  8  need 
not  be  greater  than  unity  to  detect  LCP  waves  in  the  range  <u  =  10  to  100  radians 
per  second.  For  RCP  waves  the  results  are  even  more  favorable.  It  can  be 
demonstrated  that  for  very  low  frequencies,  in  the  VLF  band  of  co  <  10~2  radians 
per  second,  the  curves  remain  together  and  achieve  a  slope  of  minus  two  as  <u 
tends  to  zero.  This  implies  that  if  any  portion  of  the  perturbed  field  can 
be  detected  at  these  frequencies,  the  wave  will  be  found  to  exist  according  to 
the  requirement  (2-131). 

Since  B0  «  2-75  (10"5)  weber/m2  and  /3  **1.1*3  (10-10)  weber/m2  at  300  km, 
the  strong  field  criteria  of  Watanabe  for  negligible  damping  can  be  used.  For 

O 

<w»l.l*5  (10-J)  radians  per  second,  the  plasma-neutral  coupling  is  weak.  There¬ 
fore,  Alfven  waves  with  <u  in  the  range  a  =  O.85  to  o>0  =  1.23  (lO2)  radians  per 
second  can  be  excited  by  a  suitable  technique.  Considering  the  case  where  the 
coupling  is  strong,  extremely  low  frequency  (ELF)  distrubances  occur  for  any 
<u«l.l*5  (10'3)  radians  per  second.  The  computed  curves  of  Figures  2-15  and  2-16 
begin  in  the  VLF  range.  In  the  ELF  range  the  traces  are  qualitatively  the 
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IONOSPHERIC  PLASMA  WAVE  ATTENUATION  vs  FREQUENCY 

(300  KM  ALTITUDE) 


DISTURBANCE  FREQUENCY,  to  -  RADIANS  PER  SECOND 


REPORT  NO.  A219 
30  NOVEMBER  1963 


IVIMD  Wave  Investigation 

1 1  -  theory 

as  those  for  the  HIT  reference  condition,  in  the  VLF  range,  Figure  2-2. 

2.7*2  Alfven  Wave  Characteristics  in  the  Ionosphere.  -  An  infinite  gas  model 
has  been  assumed  and  the  conditions  at  the  reference  altitude  have  been  used  in 
equation  (2-125)  to  obtain  the  curves  of  Figures  2-15  and  2-1 6.  The  attenuation 
distance  d,  while  reasonably  long  even  for  large  a>  ,  varies  appreciable  with 
wave  number  as  indicated  in  the  figures.  The  curves  begin  in  a  region  of 
constant  d,  which  is  the  Alfven  wave  region  for  VLF  disturbances .  This  region 
has  ended  before  =  <u0  =  123  radians  per  second  for  both  curves ..  When  o 
approaches  <d0  the  LCP  curve  loses  its  monotonic  character.  Above  =  10^  radians 
per  second,  the  LCP  curve  for  the  ionosphere  has  the  same  characteristics  as  the 
LCP  curve  for  the  HIT  reference  conditions.  Their  RCP  curves  are  also  similar 
when  the  slope  reversal  for  <u  >  107  is  considered.  Again,  for  Alfven  waves,  o>0 
is  near  the  transition  point. 


The  ionosphere  may  be  represented  to  seme  degree  by  an  arrangement  of 


channels,  or  ducts,  whose  axes  are  defined  by  the  relation: 


{■8 


1(b'o-V)'bo 

Mo 


(2-167) 


where  p  is  the  partial  pressure  of  the  charged  particles.  Equation  (2-167) 


expresses  the  requirement  that  no  net  force  be  exerted  on  a  fluid  element, 


that  is: 


A-*- 

Dv  =0 


(2-168) 


(c.f.  Appendix  A,  equation  (A-22),  without  the  displacement  current).  If  the 


spatial  rate  of  change  of  B0  along  an  arbitrary  flux  line  is  negligible, 
(2-167)  can  be  integrated  to  yield: 

r2 


(2-169) 


The  usual  concept  of  the  magnetic  field  being  frozen  in  the  fluid  is  associated 
with  (2-169).  The  right  side  of  (2-I67)  can  be  evaluated  by  means  of  a  magnetic 
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field  nodel  such  as  a  dipole,  or  it  can  be  determined  semi-empirically  by  using 
satellite  and  rocket  probe  data  to  evaluate  the  coefficients  in  a  spherical 


harmonic  expansion  of  B0 .  It  is  found  in  the  ionosphere  that: 


(2-170) 


hence,  using  (2-169),  the  waves  will  propagate  along  lines  of  constant  flux  density. 

If  the  phase  velocity  for  the  ionospheric  reference  conditions  is  examined  over 
the  ELF  to  intermediate  frequency  range,  the  curve  of  Figure  2-17  results.  Up  to 
w  „  v  . ,  corresponding  to  strong  coupling,  the  phase  velocity  is  essentially  con¬ 
stant  and  is  given  hy  U#  in  the  limit  of  F  =  0.  Between  v  ni  and  v  in  the 
plasma -neutral  coupling  diminishes  until,  at  <0  »  =  UQ  for  the  plasma. 

As  the  transition  point  or  is  approached,  the  Alfven  wave  ceases  to  exist  since 
it  does  not  propagate  or  is  highly  attenuated.  For  v>a>  Q>  the  LCP  wave,  whose 
velocity  curve  essentially  matches  that  of  the  RCP  wave  when  co  <  co  Q,  continues  to 
propagate.  However,  it  soon  looses  its  physical  identity  in  that  exceeds  c, 
and  the  corresponding  wavelengths  become  excessively  large.  The  characteristics 
of  this  curve  are  discussed  qualitatively  in  Reference  2-2. 

In  the  light  of  certain  calculations  and  satellite  observations,  the  assumption 
of  thermal  equilibrium  in  the  ionosphere  may  not  be  entirely  justified  (Reference  2-9 )■ 
If  it  is  not,  the  restrictions  of  paragraph  2.4  would  have  to  be  removed.  Not 
only  would  the  development  of  the  wave  equation  and  associated  dispersion  equation 
be  much  more  laborious,  but  an  additional  relationship  to  account  for  the  mechanism 
sustaining  the  deviation  from  equilibrium  would  have  to  be  incorporated.  From 
Reference  2-9,  a  reasonable  estimate  of  the  deviation  puts  the  electron  temperature 
at  five  times  the  ion  temperature.  Since  the  mean  plasma-neutral  collision 
frequency  and  effective  scalar  magnetic  field  are  roughly  inversely  proportional 
to  the  square  root  of  temperature,  the  effect  of  the  non -equilibrium  is  probably 
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small.  The  other  quantity  affected  by  this  condition  is  the  resistivity  n0A#  > 
which  is  inversely  proportional  to  the  three -halves  power  of  the  temperature 
(Reference  2-4).  Thus,  the  resistivity  will  be  somewhat  smaller  in  the  direction 
of  B0.  U 
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3*  THE  IABORATORY  PLASMA  ENVIRONMENT 

3*1  INTRODUCTION »  -  Prior  to  and  concurrent  with  the  theoretical  study  of 
wave s  in  plasmas,  theoretical  and  experimental  investigations  of  actual  plasma 
environments  were  conducted  in  McDonnell's  laboratory  facilities.  Consideration 
of  the  HIT  flow  field  properties  originally  suggested  the  possibility  of  simulating 
ionospheric  conditions  and  conducting  plasma  wave  experiments.  A  knowledge  of  the 
gas  conditions  in  terms  of  constituents,  total  gas  density,  temperature,  ioni¬ 
zation  and  flow  velocity  and  an  estimate  of  the  gas  turbulence  was  essential  to 
the  selection  of  the  wave  theory.  Since  numerical  estimates  for  the  important 
parameters  involved  in  the  wave  theory  were  not  available,  a  considerable  amount 
of  preliminary  experimentation  was  necessary.  The  result  of  these  experiments 
was  the  development  of  both  a  theoretical  model  for  the  HIT  and  the  necessary 
instrumentation  for  plasma  wave  investigations.  Parameter  values  which  could 
not  be  measured  directly  were  obtained  by  utilizing  the  HIT  model. 

Before  Alfven  wave  experiments  could  be  sensibly  performed  in  the  shock  and 
arc  discharge  tubes,  quantitative  results  from  the  HM  wave  guide  theory  were 
required.  Again  the  problems  of  parameter  estimation  and  appropriate  instru¬ 
mentation  arose.  Preliminary  experimentation  was  performed  in  each  apparatus, 
and  a  literature  survey  was  conducted  in  order  to  establish  suitable  models. 
Application  of  the  models  to  determine  the  parameter  values  of  the  theory  resulted 
in  certain  changes  in  the  two  devices. 

The  remainder  of  this  volume  is  devoted  to  the  description  of  the  plasma 
environment  in  the  HIT,  the  electromagnetic  shock  tube  and  the  sure  discharge 
tube,  with  the  particular  goal  of  determining  the  extent  of  their  similarity  to 
the  ionospheric  environment. 
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The  generation  of  a  plasma  environment,  the  control  of  that  environment,  and 
the  accurate  measurement  of  the  plasma  conditions  are  active  areas  of  present  day- 
study.  These  investigations  have  largely  neglected  the  study  of  plasmas  for 
simulating  ionospheric  conditions  because  of  the  particular  difficulty  of  scaling 
the  size,  ionization  percentage,  and  low  gas  density.  As  far  as  Alfven  wave 
experimentation  is  concerned,,  the  only  published  data  has  been  obtained  in  liquid 
metals  and  nearly  fully  ionized  gases,  none  of  which  can  be  easily  related  to 
ionospheric  conditions.  This  section  therefore  details  the  theoretical  considerations 
which  were  necessary  in  designing  the  plasma  environments,  and  in  determining  their 
plasma  characteristics. 

Section  3*2  outlines  the  general  characteristics  of  the  ionosphere  which  were 
to  be  simulated,  including  the  important  gas  and  magneto-ionic  conditions  and 
the  method  by  which  each  of  these  parameters  was  derived.  Section  3-3  discusses 
the  importance  of  various  parameters  for  the  experimental  study,  and  Section  3.4 
treats  the  methods  by  which  these  parameters  may  be  generated. 

The  major  theoretical  activity  in  the  development  of  a  laboratory  plasma 
environment  was  the  study  of  the  ionized  flow  in  the  hypervelocity  impulse  tunnel. 

Such  a  tunnel  had  not  previously  been  used  specifically  for  its  ionized  flow,  thus 
a  major  effort  was  required  to  predict  the  ionized  flow  conditions.  The  results  of 
this  study  indicated  two  highly  ionized  regions  of  flow,  one  region  behind  the 
starting  shock  and  a  second  region  in  the  early  portion  of  the  blow-down.  However, 
wide  variation  had  to  be  allowed  in  the  parameter  values,  because  the  complicated 
physical  and  chemical  processes  were  not  fully  explained  by  theory.  The  detailed 
analyses  of  the  conditions  in  the  arc  chamber,  in  the  shock  excited  flow  and  in  the 
blow-down  flow  are  discussed  in  Section  3*4.2. 5  and  3. 4. 2. 6. 

During  the  course  of  the  program  the  uncertainties  in  the  tunnel  flow 
conditions  and  in  the  experimental  results  made  a  redirection  of  effort  necessary. 
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The  economics  of  operating  a  large  and  highly  complex  HIT  facility  demanded  that 
only  a  few,  carefully  prepared  developmental  experiments  be  programmed.  To  prepare 
for  such  shots  the  electromagnetic  shock  tube  and,  later,  the  arc  discharge  tube 
were  constructed  and  employed.  Section  3-^-3  discusses  the  characteristics  of 
these  two  plasma  generators.  The  arc  discharge  tube  experiments  proved  to  be  of 
most  use  to  the  program. 

3.2  GENERAL  CHARACTERISTICS  OF  THE  IONOSPHERE.  -  The  upper  portion  of  the 
atmosphere  has  a  number  of  ionized  regions  which  make  up  the  ionosphere.  These 
regions  are  created  through  the  interaction  of  solar  radiation  and  the  atmospheric 
constituents  which  is  called  photo-ionization.  Due  to  the  nature  of  the  process 
and  the  variation  in  atmospheric  density  and  composition,  there  exist  several 
regions  of  relative  maximum  ionization  at  any  instant,  this  having  been  confirmed 
by  sounding  experiments.  The  lack  of  a  complete  understanding  of  the  processes 
involved  and  a  firm  knowledge  of  the  spectral  distribution  of  the  impinging  solar 
radiation  makes  It  impossible  to  predict,  theoretically,  the  observed  electron 
density  distribution. 

A  typical  distribution  of  ion  densities  as  obtained  from  radio  measurements 
indicates  three  distinct  layers  of  ionization,  the  E  layer  beginning  at  about 
100  kilometers,  the  FI  layer  at  200  kilometers,  and  the  F2  layer  at  4-50  kilometers. 
A  fourth  layer,  the  D  layer,  exists  only  in  the  daytime  and  can  be  detected  by  its 
large  absorption  of  radio  waves.  Rocket  measurements  during  periods  of  high 
absorption  have  determined  the  position  of  the  D  layer  to  be  in  the  altitude  range 
from  60  to  60  kilometers.  These  layers  may  act  as  dielectric  wave  guides  in 
channeling  MHD  waves,  changing  however,  as  a  function  of  season  and  time  of  day. 

In  many  parts  of  the  ionosphere  additional  ionized  formations,  which  are 
stable  for  several  hours,  frequently  occur,  and  extend  from  ten  to  hundreds  of 

f  i 

kilometers.  These  additional  layers,  which  do  not  appreciably  affect  the  behavior 
of  the  rest  of  the  ionosphere,  are  called  sporadic  layers.  The  most  regularly 
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occurring  of  these  is  the  ESpQr  layer,  located  slightly  above  the  E  layer.  It 
exists  frequently  enough,  that  from  25  to  90  percent  of  the  time,  long  range  propa¬ 
gation  in  the  15  mcps  band  is  made  possible  in  the  mid  latitudes.  A  sporadic  E2 
layer  can  also  be  observed,  but  less  frequently.  During  solar  flare  disturbances, 
radio  sounding  of  the  ionosphere  is  made  useless  by  abnormally  high  absorption  in 
the  D-layer,  because  of  increased  electron  concentrations  in  this  region.  Also 
magnetic  storms  greatly  affect  the  conditions  in  the  ionosphere. 

The  electron  concentration  profile  above  the  maximum  of  the  F2  region  is 
known  only  from  measurements  made  by  artificial  earth  satellites.  The  slope  of 
the  curve  of  electron  concentration  versus  altitude  has  been  found  to  be  smaller 
above  the  position  of  the  maximum  than  below  it. 

The  ionosphere  should  not  be  thought  of  as  a  quiet  medium,  but  rather  a 
turbulent  one.  This  is  particularly  evident  from  the  rapid  fluctuations  of  radio 
returns  from  the  ionosphere.  In  addition  to  the  turbulence,  the  ionosphere 
exhibits  a  regular  drift,  the  various  ionized  formations  moving  between  10  and 
1200  meters  per  second. 

The  standard  ionospheric  model  chosen  for  the  MHD  experiment,  Tables  3-1 
and  3-2,  give  an  average  value  of  the  important  parameters  required  for  computing 
the  wave  propagation  and  wave  generation  characteristics.  These  important  parameters 
include  both  non-ionized  gas  conditions,  Table  3-1,  which  are  taken  from  the  ARDC 
Model  Atmosphere  (Reference  3-1 ),  and  the  magneto-ionic  condition,  Table  3-2,  which 
are  computed.  The  following  paragraphs  indicate  the  meaning  and  give  a  derivation 
of  each  of  the  pertinent  parameters. 

Atmosphere  Parameters  -  Altitude,  as  used  in  this  report,  is  the  geometric 
altitude  measured  from  sea  level.  The  altitude  used  in  the  ARDC  calculations  of 
gas  conditions  is  the  geopotential  altitude,  which  is  related  to  the  potential 
energy  of  a  unit  mass  at  that  altitude.  Because  of  the  change  in  gravity  as  a 
function  of  altitude,  this  measurement  differs  from  the  geometric  altitude. 
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TABLE  3-2 

MAGNETO-IONIC  GAS  CONDITIONS 


GEOMAGNETIC  ALTITUDE 
(KILOMETERS) 

100 

150 

200 

250 

300 

MAGNETIC  FIELD  (GAUSS) 

.300 

.293 

.287 

.281 

.275 

ION  DENSITY  (PART/CM3) 

1.3  x  105 

2.8  x  105 

4  x  105 

1.1  x  106 

1.8  xlO6 

PLASMA  FREQ.  (RAD/SEC.) 

2.05  x  107 

3  x  107 

3.6  x  107 

5.9  x  107 

7.5  x  107 

IONIZED  TO  NEUTRAL 

PARTICLE  RATIO 

2.4  x  10-8 

9  x  10-6 

5.4  x  10“5 

4.2  x  10-4 

1.7  xlO-3 

ELECTRON  CYCLOTRON  FREQ. 
(RAD./SEC.) 

5.26  x  1046 

5.15  x  10+6 

5.04  x  1046 

4.94  x  10+6 

4.84  x  104* 

ION  CYCLOTRON  FREQ. 

(RAD./SEC.) 

1.79  x  10+2 

1.75  x  10+2 

• 

1.71  x  10+2 

1.68  x  10+2 

1.64  x  102 

AVERAGE  ELECTRON  SPEED 
(CM/SEC.) 

8.95  x  106 

2.07  x  107 

2.34  x  107 

2.34  x  107 

2.35  x  107 

AVERAGE  ION  SPEED 
(CM/SEC.) 

5.23  x  104 

1.21  x  105 

1.37  x  105 

1.37  105 

1.38  x  105 

ION  NEUTRAL 

MEAN  FREE  PATH  (CM) 

3.12  x  101 

5.44  x  103 

2.29  x  1044 

6.45  x  1044 

1.63  x  10+5 

ELECTRON  NEUTRAL 

MEAN  FREE  PATH  (CM) 

1.75  x  102 

3.06  x  104 

1.29  x  105 

3.64  x  105 

9.1xl05 

ELECTRON-ELECTRON  PARTICLE 

MEAN  FREE  PATH  (CM) 

8.17  x  103 

1.085  x  TO5 

1.25  x  10S 

4.5  x  104 

2.9  x  104 

ION-ELECTRON  MEAN  FREE 

PATH  (CM) 

1.14  x  104 

1.51  x  105 

1.75  x  10+5 

6.4  x  1044 

4.0  x  10+4 

ION-ION  MEAN  FREE  PATH 

(CM) 

8.17  x  103 

1.085  x  105 

1.25  x  10S 

4.5  x  104 

2.9  x  104 

COLLISION  FREQ.  ION 

NEUTRAL  (ENCOUNTER S/SEC.) 

1.68  x  103 

22.3 

6.0 

2.1 

8.5  x  10_1 

COLLISION  FREQ.  ELECTRON 

NEUTRAL  (ENCOUNTERS/SEC.) 

5.11  x  104 

6.77  x  102 

1.8  x  102 

6.4  x  101 

2.6  x  101 

COLLISION  FREQ.  ELECTRON 

ION  (ENCOUNTERS/SEC.) 

7.85  x  102 

1.37  x  102 

1.3  x  102 

3.7  x  102 

5.9  x  102 

COLLISION  FREQ.  ELECTRON- 

ELECTRON  (ENCOUNTERS/SEC.) 

1.1  x  103 

1.9  x  102 

1.87  x  102 

5.20  x  102 

8.10  x  102 

COLLISION  FREQ.  ION-ION 
ENCOUNTERS/SEC.) 

6.4 

1.1 

1.1 

3.0 

4.8 

EFFECTIVE  ELECTRON 

COLLISION  FREQ. 
ENCOUNTERS/SEC. 

5.3  x  104 

103 

5.0  x  102 

9.6  x  102 

1.4  x  103 

EFFECTIVE  ELECTRON  MEAN 

FREE  PATH  (CM) 

1.69  x  102 

2.07  x  104 

4.7  x  104 

2.4  x  104 

1.6  x  104 

EFFECTIVE  ION  COLl  ISION 

FREQUENCY  (ENCOUNTER/SEC.) 

1.68  x  103 

23.4 

7.1 

5.1 

5.6 

EFFECTIVE  ION  MEAN  FREE  PATH 

31.1 

5.18  x  103 

1.9  x  1044 

2.7  x  104 

2.5  x  104 

DEBYE  SHIELDING  RADIUS  (CM) 

2.74  x  10-1 

4.33  x  10_1 

4.1  x  10”’ 

2.5  x  10_1 

1.9  x  10-1 
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TABLE  3-2 

3-IONIC  GAS  CONDITIONS 


100 

150 

200 

250 

300 

350 

400 

450 

500 

.300 

.293 

.287 

.281 

.275 

.269 

.263 

.257 

.251 

1.3  x  105 

2.8  x  105 

4  x  105 

1.1  x  106 

1.8  x  106 

1.8  x  104 

1.2  x  106 

8.3  x  105 

6.9  x  10s 

2.05  x  107 

3  x  107 

3.6  x  107 

5.9  x  107 

7.5  x  107 

7.5  x  107 

6.2  x  107 

5.1  x  107 

4.7  x  107 

2.4  x  10-8 

9  x  10~6 

5.4  x  10-5 

4.2  x  10-4 

1.7  x  10“3 

4.2  x  10-3 

6.4  x  10~3 

9.3  x  10-3 

1.6  x  10“2 

5.26  x  1044 

5.15  xlO44 

5.04  x  1044 

4.94  x  1044 

4.84  x  1044 

4.72  xlO6 

4.62  x  106 

4.51  x  106 

4.40  x  104 

1.79  x  10+2 

1.75  x  10+2 

1.71  x  10+2 

1.68  x  10+2 

1.64  x  102 

1.60  x  102 

1.57  x  102 

1.53  x  102 

1.50  x  102 

8.95  x  106 

2.07  x  107 

2.34  x  107 

2.34  x  107 

2.35  x  107 

2.38  x  107 

2.41  x  107 

2.46  x  107 

2.51  x  107 

5.23  x  104 

1.21  x  105 

1.37  x  105 

1.37  105 

1.38  x  105 

1.39  x  105 

1.41  x  105 

1.44  x  105 

1.47  x  105 

3.12  x  101 

5.44  x  103 

2.29  x  1044 

6.45  x  1044 

1.63  x  10+5 

3.93  x  10+S 

9.04  x  10+5 

1.90  x  1044 

3.85  x  1044 

1.75  x  102 

3.06  x  104 

1.29  x  105 

3.64  x  105 

9.1  x  105 

2.22  x  106 

5.1  x  106 

1.07  x  107 

2.15  xlO7 

8.17  x  103 

1.085  x  105 

1.25  x  105 

4.5  x  104 

2.9  x  104 

3.0  x  104 

4.8  x  104 

7.2  x  104 

105 

S’) 

1.14  x  104 

1.51  x  105 

1.75  x  10+5 

6.4  x  1044 

4.0  x  1044 

4.2  x  1044 

6.6  x  1044 

1.03  x  10+5 

1.4  x  10+5 

8.17  x  103 

1.085  x  105 

1.25  x  105 

4.5  x  104 

2.9  x  104 

3.0  x  104 

4.8  x  104 

7.2  x  104 

105 

1.68  x  103 

22.3 

6.0 

2.1 

8.5  x  10-1 

3.6  x  10_1 

1.6  x  10_1 

7.6  x  10“2 

3.8  x  10~2 

5.11  x  104 

6.77  x  102 

1.8  xlO2 

6.4  x  101 

2.6  x  101 

1.1  x  101 

4.7 

2.3 

1.1 

7.85  x  102 

1.37  x  102 

1.3  x  102 

3.7  x  102 

5.9  x  102 

5.7  x  102 

3.6  x  102 

2.4  x  102 

1.8  x  102 

1.1  x  103 

1.9  x  102 

1.87  x  102 

5.20  x  102 

8.10  x  102 

7.95  x  102 

5.02  x  102 

3.42  x  102 

2.51  x  102 

6.4 

1.1 

1.1 

3.0 

4.8 

4.6 

2.9 

2.0 

1.5 

5.3  x  104 

103 

5.0  x  102 

9.6  x  102 

1.4  x  103 

1.4  x  103 

8.6  x  102 

5.8  x  102 

4.3  x  102 

1.69  x  102 

2.07  x  10^ 

4.7  x  104 

2.4  x  104 

1.6  x  104 

1.7  x  104 

2.8  x  IQ4 

4.2  x  104 

5.8  x  104 

1.68  x  103 

23.4 

7.1 

5.1 

5.6 

4.9 

3.1 

2.0 

1.5 

31.1 

5.18  x  103 

1.9  x  1044 

2.7  x  104 

2.5  x  104 

2.8  x  104 

4.6  x  104 

7.2  x  104 

1.0  x  105 

2.74  x  VO-1 

4.33  x  10-1 

4.1  x  10_1 

2.5  x  10_1 

1.9  x  10_1 

2.0  x  10-1 

2.4  x  10-1 

3.0  x  10-1 

3.3  x  10-1 
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The  temperatures  at  altitudes  inaccessible  to  balloons  are  derived  from 
measured  pressures  or  densities.  The  ARDC  calculations  utilize  a  molecular-scale 
temperature  that  is  independent  of  molecular  weight.  Since  the  molecular  weight 
is  not  well  known  at  altitudes  above  90  kilometers,  the  molecular-scale 
temperature,  Tj^,  gives  a  more  accurate  indication  of  temperature.  It  is  defined 
by 


where 

T  »■  temperature  in  °K 

M  =  molecular  weight 

Mq  =  sea-level  value  of  molecular  weight 

Since  (T/M)  is  directly  proportional  to  the  square  of  the  particle  kinetic 

energy,  a  measurement  of  the  energy  enables  calculation  of  (t/m)  and  therefore,  Tm 

Molecular  weight  is  a  dimensionless  quantity  in  which  the  naturally  occurring 
mixture  of  oxygen  isotopes  has,  by  definition,  a  value  of  l6.  The  atmosphere 
defined  by  this  model  is  assumed  dry,  with  a  sea  level  molecular  weight  Mq  equal 
to  28.966,  The  mean  molecular  weight  curve  from  the  1959  model  results  from  an 
attempt  to  make  use  of  the  best  available  experimental  data,  and  is  about  10  to 
20  percent  heavier  than  that  of  the  1956  model. 

In  Table  3-1  the  number  density  is  computed  from  the  density  and  mean 
molecular  weight.  The  particle  speed  is  the  mean  of  a  Maxwellian  distribution 
of  speeds  for  all  air  particles  within  a  given  element  any  volume,  assuming  that 
all  air  molecules  have  the  mass  associated  with  the  mean  molecular  weight. 

The  mean  free  path  is  the  mean  value  of  the  distance  traveled  by  each  of  the 
molecules  between  collisions  with  other  molecules.  The  expression  for  mean  free 
path,  L,  is  obtained  from  kinetic  gas  theory,  assuming  elastic  collisions  between 
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spherical  molecules  and  is  given  by: 


L- 


where : 

a  =  average  effective  collision  diameter 
n  =  particle  density 

The  mean  collision  frequency  of  the  molecules  is  the  ^verage  velocity  of  the 
molecules  divided  by  their  mean  free  path. 

Magneto-ionic  Parameters  -  The  magneto -ionic  properties  of  the  ionosphere, 

Table  3-2,  were  not  given  for  the  ARDC  atmosphere  and,  therefore,  had  to  be 
calculated.  The  magnetic  field  strength  has  been  extrapolated  from  calculated 
field  strengths  at  the  magnetic  equator  assuming  a  dipole  source.  Original 
calculations  gave  a  field  strength  of  0.3116  gauss  at  the  earth's  surface,  and  a 
value  of  0.2341  gauss  at  an  altitude  of  0.1  earth  radii.  The  field  gradient  was 
assumed  constant  between  these  values. 

Ion  densities  were  available  from  many  sources.  Those  densities  used  in 
the  model  atmosphere  of  Figure  3-1  were  obtained  from  a  composite  of  rocket  and 
satellite  measurements  in  mid  latitudes  during  the  afternoon.  The  plasma  frequency, 
*»  ,  is  a  direct  measure  of  the  electron  concentration  and  is  defined  by: 


where : 

ng  =  electron  concentration  in  particles/cm^ 
e  =  the  charge  on  the  electron  (4.8  x  10~10  e.s.u. ) 

=  the  mass  of  the  electron  in  grams. 

The  cyclotron  frequency  of  charged  particles,  u  c,  is  the  angular  oscillation 
rate  of  the  particles  around  a  magnetic  line  of  force  given  as: 

<Bo 
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LOO  ELECTRON  DENSITY  -  ELECTRONS/CM3 


FIGURE  3-1 
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where: 

B 

q. 


.10, 


•  the  magnetic  field  strength  in  gauss 
>■  the  charge  on  the  particle 
m  -  the  mass  of  the  particle  in  grams 
c  *  the  speed  of  light  in  cm/sec  (3  x  10xw ) • 

The  ion  cyclotron  frequency  has  been  computed  for  the  atomic  oxygen  ion  which 
was  assumed  to  be  the  predominant  positive  ion. 

The  ion  suid  electron  average  speeds  are  found  from  the  kinetic  theory  of 
gases  to  be  equal  to: 


2 

V“ 


Yt 


where: 

V  *  the  average  speed  in  centimeters  per  second 

k  ■  Boltzmann's  constant  (1*38  x  10-1^  ergs/deg  K) 

T  =»  temperature  in  °K 

m  *  mass  of  the  particle. 

For  molecules  of  the  same  size,  one  molecule  will  encounter  another  molecule, 
if  its  center  is  at  a  distance  less  than  one  molecular  diameter,  d,  from  the  line 

generated  by  the  motion  of  the  center  of  the  moving  particle.  Thus,  the  area  for 

2  _  2 
collision  will  be  n  d  ,  and  the  volume  of  collision  will  be  V  n  d  .  If  the 

average  concentration  of  molecules  is  n  per  cubic  centimeter,  the  number  of 

—  ,  2 

collisions  per  second  will  be  nV  n d  . 

The  mean  free  path,  L,  is  the  average  distance  traveled  in  one  second, 

divided  by  the  number  of  collisions  per  second.  For  molecules  of  equal  size: 

L  -X_  1 

n  V  it  n  nr 

For  the  mean  free  path  of  two  different  types  of  particles  having  masses  m-^  and 
m2,  radii  r^  and  Tq,  and  mean  speeds  and  V2,  the  expression  for  mean  free  path 
of  a  particle  of  type  1  colliding  with  particles  of  type  2,  is: 
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l 


where  =  rl  +  r2  ant^  n2  the  number  of  particles  of  type  2  per  cubic  centimeter. 
When  Vq  equals  Vg  this  equation  reduces  to: 


For  a  collision  between  two  neutral  particles  or  between  a  neutral  particle  and 
a  charged  particle,  the  effective  diameter  is  of  the  order  of  10  centimeters. 

The  AKDC  model  atmosphere  uses  an  average  effective  collision  diameter  of 

-8 

3- 65  x  10 ”  centimeters.  For  these  calculations  this  number  is  a  constant  for 
ions  and  neutral  particles.  The  assumption  of  equal  effective  velocities,  while 
valid  for  ion  and  neutral  particle  encounters,  is  not  true  for  neutral  and 


electron  encounters.  If  the  electron  is  considered  as  a  point  particle  moving 
at  a  relatively  high  speed  through  the  gas,  the  cross-sectional  area  of  the  cylinder 
of  interception  is  just  the  cross-sectional  area  of  one  molecule,  which  results  in 
a  mean  free  path  four  times  as  large  as  that  of  the  molecules.  Assuming  the 
neutral  particle  velocity  to  be  negligible  compared  with  the  electron  velocity, 
the  mean  free  path  of  electrons  with  neutrals  is  found  to  be: 

1 

L  « -  or  5.64  L. 

e  o  n 

n2  n  rj 

where  1^  is  the  mean  free  path  of  neutrals  with  neutrals.  The  mean  free  paths  for 
electron-ion  and  ion-ion  pairs  can  be  calculated  by  considering  the  Coulomb 
interaction.  For  a  collision  between  an  electron  and  an  ion,  the  effective 
diameter  is  of  the  order  of  10"^  Z^  (300/T)  centimeters,  where  Z^  is  the  number 
of  unit  charges  on  the  ion.  For  a  collision  between  two  ions,  it  is  about 
10  ^  Z i  Zj  (300/T),  where  Z^  and  Zj  are  the  number  of  unit  charges  on  the  two 
different  ions  (Reference  3-2).  The  collision  frequencies  for  each  pair  of 
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particles  is  simply  the  ratio  of  the  average  particle  velocity  to  the  mean  free 
path. 

A  The  effective  collision  frequency  is  the  total  number  of  collisions  with  all 

particles  per  unit  time.  The  effective  mean  free  path  is  the  average  length  of 
path  traveled  between  collisions  with  all  particles,  and  it  is  computed  from  the 
particle  velocity  and  ratio  of  collision  frequency  to  mean  free  path  relationship. 

The  Debye  shielding  distance,  h,  is  a  measure  of  the  distance  over  which 
the  electron  concentration  can  deviate  appreciably  from  the  positive  charge 
concentration,  and  is  defined  as: 


where  T  is  the  temperature  in  °K  and  ne  is  the  electron  concentration  in 
particles  per  cubic  centimeter. 

The  Alfven  velocity  is  defined  in  the  cgs  system  of  units  by: 

U0«B„ 

o  o 

where : 

2 

P±  ■  density  of  the  ions  in  gm/cm 

Bo  *  magnetic  flux  density  in  gauss. 

3-3  SIMUIATION  PARAMETERS.  -  The  experimental  environment  can  only  partially 
simulate  and,  thus,  limits  the  number  of  Alfven  wave  parameters  that  may  be  studied. 
This  is  because  the  low  ionization  concentrations,  the  small  magnetic  fields  and 
the  low  densities  of  the  ionosphere  would  normally  result  in  very  large  charac¬ 
teristic  dimensions  for  these  waves.  Laboratory  experiments  are  limited  to  small 
characteristic  dimensions;  typical  dimensions  prior  to  the  initiation  of  the 
McDonnell  experiments,  were  a  small  fraction  of  a  meter.  The  characteristic 
dimension  of  the  McDonnell  HIT  experiments  was  approximately  one  meter.  The 
following  paragraphs  indicate  the  importance  of  the  major  parameters  and  how  each 
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influences  the  experimental  design. 

In  the  ionosphere  two  types  of  Alfven  wave  propagation  may  exist:  a  plane 
wave  traveling  along  a  magnetic  field  line,  or  a  guided  wave  traveling  between  the 
dielectric  boundaries  of  various  ionospheric  layers.  In  the  laboratory  only  guided 
waves  exist,  because  of  the  boundaries  of  the  confined  plasma.  To  obtain  plane 
wave  propagation,  the  attenuation  across  the  magnetic  field  lines  must  be  great 
enough  and  the  wave  generation  must  be  spatially  confined  to  a  region  sufficiently 
far  from  the  walls  to  eliminate  the  boundary  effects.  Therefore,  in  practice,  a 
plasma  of  large  cross  section  is  required  to  propagate  plane  waves.  In  wave  guide 
propagation,  transverse  magnetic  (TM)  and  transverse  electrical  (TE)  modes  may  exist. 
In  general,  both  of  these  modes  have  cutoff  frequencies  which  decrease  with 
increasing  experimental  cross  section;  and  for  a  fixed  frequency,  the  attenuation 
decreases  for  increasing  cross  section.  The  physical  size,  therefore,  determines 
the  type  of  propagation  that  may  exist,  its  minimum  frequency,  and  its  rate  of 
attenuation. 

6 

The  ion  concentration  in  the  ionosphere  is  a  maximum  of  about  2  x  10°  ions 

per  cubic  centimeter  in  the  F2  region.  The  corresponding  neutral  particle  density 

q 

is  about  10  particles  per  cubic  centimeter,  or  the  ionosphere  is  about  0.1  percent 
ionized.  The  ionization  concentration  determines  the  velocity  of  propagation  of 
the  wave,  that  is,  the  Alfven  velocity.  The  wave  number,  as  determined  by  the. 
propagation  velocity,  is  influential  in  determining  attenuation.  Wave  attenuation 
is  inversely  proportional  to  the  conductivity  and  hence  the  ionization  concen¬ 
tration,  and  is  also  affected  by  ion-neutral  particle  damping,  which  in  the  first 
approximation,  is  proportional  to  the  ion-neutral  collision  frequency  and  inversely 
proportional  to  the  Alfven  velocity.  The  desired  condition  of  small  wave  damping 
in  a  partially  ionized  gas  cannot  be  met  in  the  laboratory,  due  to  the  normally 
high  neutral -ion  collision  frequency  and  low  conductivity.  However,  the  require- 
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sent  for  high  conductivity  can  be  relaxed  if  low  disturbance  frequencies  and 
large  plasma  dimensions  are  used.  The  only  possibility  of  achieving  low  damping 
in  a  partially  ionized  gas,,  is  through  the  generation  of  a  low  density  plasma  or  a 
very  cold  plasma. 

The  magnetic  field  intensity  influences  the  propagation  velocity,  the  wave 
number,  and  the  wave  length,,  and  therefore,  the  attenuation.  The  magnetic  field 
determines  the  ion  cyclotron  frequency,  which  sets  the  upper  limit  of  disturbance 
frequency.  For  Alfven  waves  to  predominate  over  acoustical  disturbances  in  a  plasma, 
tne  magnetic  pressure  \^Jm\ist  greatly  exceed  the  hydrostatic  pressure.  The 
amplitude  of  the  wave,  in  theoretical  considerations,  is  usually  assumed  to  be 
small  with  respect  to  the  ambient  parameters. 

In  none  of  the  theoretical  studies  of  this  program  were  ionization  and 
recombination  processes  included.  Because  of  the  long  recombination  times 
large  regions  of  the  ionosphere  can  not  change  appreciably  within  the  period 
of  the  Alfven  disturbance.  Thus,  frozen  conditions  for  the  characteristic 
times  associated  with  experiments  must  be  maintained.  Other  considerations 
in  the  choice  and  design  of  a  plasma  environment  are: 

1.  Diagnostic  equipment  must  be  available  to  define  plasma  and  experimental 
conditions1  easily  and  accurately. 

2.  The  plasma  generation  mechanism  should  be  completely  decoupled  from  the 
wave  excitation  mechanism. 

3.  The  propagation  and  the  excitation  of  the  wave  should  be  capable  of 
being  studied  independently. 

4.  Control  of  the  simulation  should  be  possible  for  a  wide  range  of 
conditions . 
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3.4  SIMULATION  FACILITIES 

3.4.1  General.  -  The  generation  and  propagation  of  Alfven  waves  have 
been  reported  in  liquid  metals  by  Lundquist  (Reference  3-3.)  and  Lehnert 
(Reference  3-4)  and  in  gaseous  plasmas  by  Wilcox  (References  3-5 >  3-6  and 
3-7)  and  Jephcott  (References  3-8  and  3-9).  The  liquid  metal  experiments 
certainly  cannot  be  scaled  to  the  ionospheric  conditions.  The  gaseous  plasma 
experiments  are  very  pertinent  and  provide  some  of  the  best  MHD  experiments 
conducted  to  date.  These  experiments  have  measured  for  high  frequency 
Alfven  waves  (500,000  cps),  the  phase  velocity,  attenuation,  energy  distri¬ 
bution  and  reflection  characteristics.  The  waves  generated  were  torsional 
waves,  produced  by  a  high  current  coaxial  discharge.  The  gas  used  for  most 
of  the  experiments  was  hydrogen,  at  a  pressure  of  100  microns  of  Hg.  When 
waves  were  observed,  the  gas  was  nearly  fully  ionized  and  the  magnetic  field 
strength  was  above  7000  gauss.  The  attenuation  effect  of  high  ion- neutral 
collision  frequency,  which  prohibits  the  study  of  the  waves  in  a  small 
volume,  partially  ionized  gas,  can  be  reduced  in  a  facility  such  as  a  Hyper¬ 
velocity  Impulse  Tunnel  (HIT)  since  it  is  adaptable  to  the  generation  of 
large  plasma  environments. 

The  Hypervelocity  Impulse  Tunnel  (Figure  3-2)  is  a  facility  having  a 
test  section  50  inches  in  diameter,  at  a  distance  of  twenty  feet  from  the 
throat  of  the  expansion  chamber.  The  HIT  achieves  its  particular  local  flow 
by  expanding  the  arc  chamber  high  pressure  and  high  temperature  gases  through 
a  nozzle,  resulting  in  a  reduction  of  the  local  temperature,  and  thereby 
lowering  the  speed  of  sound,  and  an  increase  in  the  local  flow  Mach  number, 
without  affecting  the  stagnation  temperature.  For  the  McDonnell  facility, 
Mach  numbers  up  to  27  can  be  generated  for  periods  up  to  100  milli seconds. 

The  end  of  useful  flow  occurs  when  dissociation  and  ionization  of  the  gas 
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commences.  The  purpose  of  this  program  was  to  extend  the  normal  operating 
range  of  the  HIT  into  that  required  for  ionospheric  simulation,  that  is,,  by 
operating  when  the  flow  was  partially  ionized.  Since  this  facility  had  not 
been  operated  under  such  conditions  prior  to  the  start  of  this  experimental 
program,  a  detailed  study  of  the  facility  was  made  in  order  to  optimize  the 
flow  for  the  generation  of  Alfven  waves,  their  propagation  and  measurement. 

During  the  course  of  the  program  two  additional  plasma  generators 
served  as  useful  tools  augmenting  the  tunnel  experiments.  These  facilities 
were  capable  of  inexpensive  tests  for  developing  instrumentation  and  consisted 
of  an  electromagnetic  shock  tube  and  an  arc  discharge  tube.  The  theory  of 
operation  of  both  of  these  facilities  is  given  in  Sections  3* ^*3.1  and 
3- 4. 3-2. 

3A.2  Plasma  Flow  in  the  Hypervelocity  Impulse  Tunnel 

3. If. 2.1  Introduction.  -  While  many  sources  describe  flow  conditions 
existing  in  a  hypersonic  nozzle  or  discuss  the  flow  phenomena,  few  present  a 
comprehensive  discussion  covering  all  phases  of  nozzle  or  tunnel  operation 
at  conditions  suitable  for  plasma  experimentation.  This  section  will 
describe  the  operation  and  flow  phenomena  of  the  arc  chamber  and  hypervelocity 
nozzle,  of  the  McDonnell  Hypervelocity  Impulse  Tunnel,  and  will  present  the 
thermodynamic  conditions  for  various  phases  of  operation,  comparing  several 
analytical  techniques  used  in  describing  the  flow. 

3 .4.2. 2  Tunnel  Construction.  -  The  McDonnell  Aircraft  Corporation 
Hypervelocity  Impulse  Tunnel  facility  consists  of  an  arc  chamber  in  which  a 
high  pressure,  high  temperature  gas  is  created  by  an  electrical  discharge 
from  a  capacitor  bank,  a  divergent  conical  nozzle  leading  to  a  constant  area 
test  section,  and  a  reservoir  which  contains  a  conical  diffuser  to  eliminate 
shock  reflections  through  the  test  section.  Initially,  the  arc  chamber  is 
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separated,  from  the  test  section  by  a  diaphragm  near  the  nozzle  throat  which 
ruptures  during  pressure  build-up  following  initiation  of  the  electrical  dis¬ 
charge.  A  schematic  drawing  of  the  McDonnell  tunnel  is  shown  in  Figure  3-2. 

The  arc  chamber  consists  of  an  outer  case,  designed  to  withstand 
100,000  psia,  into  which  a  cartridge  assembly  is  inserted  prior  to  each 
firing.  The  components  of  the  cartridge  assembly  are  an  internal  liner,  the 
throat  and  diaphragm  assembly,  the  electrode  assembly,  the  dump  port  assembly, 
and  fittings  for  transducers  and  gas  filling  ports.  The  electrical  discharge 
is  initiated  by  exploding  a  trigger  wire  which  runs  from  a  trigger  assembly 
within  the  center  electrode  to  the  outer  electrode.  The  exploding  wire  ionizes 
the  gas  in  the  gap  between  electrodes,  increasing  the  electrical  conductivity 
in  this  region  and  creating  a  current  path  for  initiating  the  capacitor  bank 
discharge. 

Temperatures  and  pressures  achieved  in  the  arc  chamber  are  determined  by 
the  quantity  of  energy  available  from  the  arc  discharge  and  by  the  volume  of 
gas  into  which  the  energy  is  absorbed.  The  volume  of  the  internal  liner  can 
be  varied  from  100  to  1200  cubic  inches.  The  maximum  storage  capacity  of  the 
capacitor  bank  is  7> 000, 000  joules  at  12,000  volts;  maximum  discharge  current 
is  3,600,000  amperes.  The  energy  level  of  the  discharge  is  varied  by  control¬ 
ling  either  the  voltage  or  the  number  of  capacitors  in  the  bank. 

The  diaphragm,  which  is  commonly  of  Jfylar,  is  placed  slightly  forward  of 
the  nozzle  throat.  The  nozzle  throat  diameter  is  adjustable  and  may  be  varied 
from  less  than  0.1  inch  to  greater  than  1.5  inches,  although  the  larger  sizes 
would  not  normally  be  used  for  aerodynamic  testing.  Nozzle  throat  sections 
are  commonly  machined  from  copper  or  tungsten. 

A  conical  nozzle  with  a  total  included  angle  of  10  degrees  leads  from 
the  throat  to  a  50- inch  diameter  test  section.  The  distance  from  the  nozzle 
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throat  to  this  test  section  is  approximately  23  feet,  with  ports  for  ion¬ 
ization  gages  or  other  probes  provided  at  intervals  along  the  nozzle  wall. 

The  test  section  has  a  37*5  inch  diameter  viewing  window,  which  permits  high 
speed  photographs  of  the  flow  and  microwave  interferometer  measurements  of 
electron  co,  centration. 

The  test  section  is  followed  by  a  fixed  conical  diffuser  with  an  area 
ratio  (diffuser  throat  area  to  50- inch  diameter  test  section  area)  of  0.5. 

From  the  diffuser  the  flow  passes  into  the  reservoir.  A  vacuum  system,  con¬ 
sisting  of  both  mechanical  and  oil  diffusion  pumps,  is  connected  to  the 
reservoir  and  is  capable  of  evacuating  the  total  volume  of  the  tunnel  to  1 
micron  Hg  in  15  minutes. 

3* ^.2. 3  Operating  Regimes.  -  The  analysis  of  the  HIT  flow  is  divided 
into  the  characteristic  regimes  of  the  arc  chamber,  the  starting  shock,  the 
blow-down  flow,  and  the  flow  breakdown. 

The  arc  chamber  is  first  charged  with  the  desired  gas.  As  the  capacitor 
bank  is  discharged  through  the  electrodes,  the  gas  is  heated  and  the  pressure 
increases  rupturing  the  diaphragm  near  the  nozzle  throat  and  initiating  a 
shock  wave  which  moves  downstream  into  the  nozzle.  The  arc  chamber  pressure 
then  continues  to  rise  -until  the  discharge  is  complete.  The  analysis  of  the 
arc  chamber,  discussed  in  paragraph  -j.4.2.4,  will  be  concerned  with  the  equi¬ 
librium  concentrations  of  electrons  and  ionized  species  at  arc  chamber  con¬ 
ditions,  and  with  the  changes  which  occur  as  the  arc  chamber  gas  discharges 
through  the  throat. 

Rupture  of  the  nozzle  diaphragm  initiates  a  shock  wave  which  travels 
downstream  into  the  low  pressure  gas  in  the  tunnel.  During  this  "starting 
cycle"  the  tunnel  resembles  a  shock  tube  with  a  divergent  nozzle.  The  initial 
shock  is  followed  by  the  contact  surface  and  then  by  a  second  shock  or  rare- 
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faction  wave,  depending  upon  conditions  existing  in  the  tunnel.  The  region 
between  the  initial  shock  and  the  contact  surface,  referred  to  as  the  "shock- 
excited"  region,  is  brought  to  high  temperature  and  pressure  by  the  shock. 

This  region  is  discussed  in  detail  in  paragraph  3-^-2.5  because  it  was  deemed 
useful  for  the  plasma  experiments. 

After  the  shock  system,  created  when,  the  diaphragm  bursts,  has  decayed, 
the  arc  chamber  discharges  into  the  tunnel  in  a  quasi- steady  flow  referred 
to  as  the  blow-down  flow.  This  flow  reaches  the  test  section  at  a  relatively 
high  velocity,  but  low  temperature.  Because  of  the  rapid  expansion,  the  nozzle 
flow  may  occur  under  nonequilibrium  conditions,  where  chemical  and  electron 
recombination  reactions  are  "frozen"  at  some  point  in  the  nozzle.  If  freezing 
occurs,  the  electron  concentration  may  remain  at  a  relatively  high  level 
despite  the  rapid  cooling  of  the  gas.  Flow  conditions  in  the  blow-down  region 
are  examined  in  detail  in  paragraph  3. 4. 2. 6. 

The  final  stages  of  discharge  from  the  arc  chamber  are  characterized  by 
a  region  of  unsteady  flow  referred  to  as  flow  breakdown  and  the  flow  direction 
may  reverse  as  the  result  of  an  increase  in  pressure  within  the  reservoir. 

'3.4.2. 4  Arc  Chamber  Performance.  -  The  performance  characteristics  of  the 
arc  chamber  influence  the  nature  of  the  flow  during  the  starting  cycle  and  blow¬ 
down  periods.  In  the  initial  phase,  during  discharge  of  the  capacitor  bank, 
thermodynamic  changes  occur  as  a  result  of  the  finite  time  required  for  arc 
discharge,  and  in  later  phases  these  changes  occur  as  a  result  of  mass  loss 
accompanying  expansion  of  the  gases  into  the  tunnel. 

In  order  to  obtain  the  high  ion  and  electron  concentrations  required  for 
the  propagation  of  Alfven  waves,  it  was  necessary  to  work  with  unusually  Mgh 
temperatures  in  the  arc  chamber.  As  a  consequence,  the  initial  charge 
pressures  in  the  arc  chamber,  being  of  the  order  of  25  psig,  were  much  less 
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then  those  normally  used  for  aerodynamic  testing.  For  these  tests,  the  arc 
voltage  was  lower  than  usual,  and  the  damping  time  was  correspondingly  longer. 

In  Figure  3-3  is  presented  a  typical  capacitor  bank  voltage  record  for  one  of 
the  higher  temperature  tests.  For  this  test  the  arc  voltage  during  the  first 
half  cycle  is  estimated  to  be  2.5  KV,  and  tlie  maximum  current,  determined 
graphically,  is  2.24  x  10^  amperes.  The  energy  input  during  the  first 
quarter  cycle  is  found  to  be  approximately  2.5  x  10 ^  joules.  The  stored  energy 
in  the  condenser  is  4.63  x  10^  joules.  It  is  evident  that  over  half  of  the 
available  energy  is  added  within  the  first  half  millisecond.  The  energy 
stored  in  the  capacitor  bank  is  equal  to  l/2  CE2,  bus  the  energy  delivered 
to  the  gas  is  less,  and  is  determined  by  the  peak  arc  chamber  pressure  and  the 
initial  gas  density. 

The  temperature  rise  achieved  in  the  arc  chamber  as  a  result  of  the 
electrical  discharge  is  determined  from  the  measured  arc  chamber  pressure. 

A  typical  pressure  trace  from  an  oscillograph  record  is  shown  in  Figure  3-4. 

Although  the  Norwood  strain  gage  transducers  have  a  rapid  response  time,  an 
appreciable  response  lag  is  introduced  by  a  length  of  small  tubing  containing 
several  orifices  Inserted  between  the  transducer  and  the  chamber  to  protect 
the  transducer  from  the  hot  gases.  This  lag,  of  the  order  of  15  milliseconds, 
is  clearly  evident  in  the  figure .  After  peak  pressure  is  reached,  the  pressure 
decays  at  a  steady  rate  due  to  discharge  of  gas  through  the  nozzle.  An  over¬ 
shoot  at  the  peak  value  is  often  observed.  Assuming  that  the  pressure  decay 
rate  during  the  first  10  to  15  milliseconds  is  equal  to  that  measured  during 
the  latter  phase  of  the  gas  discharge,  the  peak  arc  chamber  pressure  is  esti¬ 
mated  by  extrapolating  the  measured  decay  rate  to  its  value  at  1  millisecond 
after  initiation  of  arc  discharge.  The  extrapolation  is  illustrated  in 
Figure  3-4.  The  temperature  rise  in  the  arc  chamber  is  then  obtained  from 
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tables  of  high  temperature  gas  properties,  with  the  further  assumption  that 
energy  addition  takes  place  along  a  constant  density  path.  The  pressure  decay 
is  the  result  of  a  mass  loss  through  the  nozzle,  as  well  as  a  temperature  drop. 

The  degree  of  ionization  represented  by  the  electron  concentration  is  an 
important  thermodynamic  parameter  for  plasma  studies.  Equilibrium  electron 
concentrations  in  helium,  nitrogen,  and  air  are  presented  in  Figure  3-5  as 
functions  of  arc  chamber  temperature  and  pressure.  The  assumption  of  equili¬ 
brium  conditions  in  the  arc  chamber  appears  to  be  justified  when  the -relatively 
long  holding  time  for  the  gas  in  the  arc  chamber  is  compared  with  the  ionization 
and  relaxation  times.  Petcheck  and  Byron,  Reference  3-10,  for  instance,  report 
that  at  10,000°K  and  5  cm  of  Hg  pressure  the  time  to  reach  equilibrium  ionization 
behind  a  shock  wave  in  argon  is  of  the  order  of  100  microseconds.  At  the  higher 
densities  of  interest,  the  time  required  for  ionization  should  be  proportionately 
shorter. 

The  electron  concentration  at  the  nozzle  throat  may  be  significantly  less 
than  that  computed  for  equilibrium  conditions  in  the  arc  chamber,  as  a  result 
of  the  lowered  temperature  and  density  at  the  throat. 

The  extent  of  this  reduction  for  the  case  of  nitrogen  flow  can  be  estimated  from 
Figures  3-6,  3-7,  and  3-8  where  temperature,  pressure  and  density  ratios  at  the 
throat  are  compared  with  the  arc  chamber  conditions.  For  these  figures  the  usual 
isentropic  compressible  flow  relations  have  been  used.  Equilibrium  concentrations 
at  the  throat  can  be  found  from  Figure  3-5  by  substituting  throat  properties  for 
the  corresponding  arc  chamber  values. 

3.4.2. 5  Shock  -  Excited  Flow 

General  Description.  -  If  a  duct  is  divided  by  a  partition  into  two  sections 
in  which  a  pressure  differential  exists,  a  combination  of  wave  phenomena  will 
follow  a  sudden  removal  of  this  partition.  A  shock  wave  propagates  into  the  low 
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DENSITY  RATIO  AT  ARC  CHAMBER  THROAT 


pressure  region  and  a  rarefaction  fan  propagates  into  the  high  pressure  gas.  The 
gas  particles  that  were  initially  in  contact  with  the  partition  form  a  contact 
surface,  a  temperature  discontinuity  separating  the  two  adjacent  gases.  The 
resultant  wave  diagram  is  shown  in  Figure  3-9*  The  gas  that  is  swept  out  by 
the  shock  wave  is  compressed  and  appears  as  region  2  in  the  wave  diagram. 

If  the  duct  is  of  constant  area,  as  in  the  simple  shock  tube,  the  down¬ 
stream  propagating  shock  wave  theoretically  does  not  attenuate  as  it  moves  down 
the  tube.  On  the  other  hand,  if  the  duct  has  a  divergent  cross-section,  the 
downstream  propagating  shock  wave  decreases  in  strength,  and  a  secondary  wave 
appears,  which  propagates  in  the  upstream  direction.  One  possible  wave  system 
for  the  divergent  duct  considers  that  the  rarefaction  fan,  which  propagates 
upstream,  separating  regions  3  and  4,  coalesces  into  the  secondary  wave  of 
the  divergent  channel.  Depending  upon  the  area  ratio,  pressure  ratio,  and 
initial  shock  strength,  the  secondary  wave  may  be  either  an  upstream  propa- 
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gating  expansion  wave  or  an  upstream  propagating  shock  wave.  Although  the 
wave  propagates  upstream  in  either  case,  it  travels  downstream  relative  to 
the  tube  and  increases  in  strength  as  it  does  so.  This  is  the  wave  system 
assumed  for  this  study. 

The  basic  equations,  describing  the  wave  phenomena  used  in  analyzing  the 
shock  excited  flow  in  both  the  constant  area  shock  tube  and  the  diverging 
hypersonic  nozzle,  are  the  same.  The  analysis  presented  for  the  shock  excited 
region  is  based  mainly  on  ideal  gas  relationships,  although  some  r.eal  gas 
effects  have  been  taken  into  account  in  the  discussion  of  the  "integral" 
approach  to  the  expanding  nozzle.  Upon  diaphragm  rupture,  three  disturbances 
are  created.  First,  a  shock  wave  forms  and  travels  downstream.  Secondly,  a 
contact  surface,  which  separates  the  high  pressure  driver  gas  from  the  low 
pressure  driven  gas,  propagates  downstream  behind  the  shock  wave.  Finally, 
an  expansion  wave  is  formed,  the  head  of  which  propagates  upstream  into  the 
driver  section  and  is  reflected  downstream  from  the  end  of  the  driver  chamber 
(arc  chamber).  The  tail  of  the  expansion  wave  attempts  to  travel  upstream  into 
the  driver  section,  but  can  be  swept  downstream  by  the  gas  flowing  from  the 
driver  chamber.  Testing  normally  takes  place  in  the  shock  processed  gas 
between  the  shock  wave  and  the  contact  surface.  The  test  time  is,  therefore, 
the  duration  of  time  between  passage  of  the  shock  wave  and  passage  of  the 
contact  surface  past  a  given  test  section.  In  a  closed  tube  the  primary  shock 
is  reflected  and  different  properties  are  found  in  the  region  behind  it.  Test 
time  in  the  main  shock  processed  gas  is  limited  by  either  the  shock  wave 
reflecting  off  the  downstream  end  of  the  tube  and  interacting  with  the  contact 
surface,  or  by  the  expansion  wave,  reflected  downstream  from  the  back  of  the 
driver  chamber  interacting  with  the  contact  surface.  Because  the  length  of 
the  driven  section  (nozzle-reservoir)  is  long  compared  to  that  of  the  driver 
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section  in  most  shock  tubes,  the  expansion  wave-contact  surface  interaction 
does  not  occur. 

The  theoretical  speed  of  the  shock  generated  at  diaphragm  rupture  may  be 
determined  from  the  stagnant  conditions  existing  on  both  sides  of  the  diaphragm 
before  rupture.  The  pressures  and  velocities  in  regions  2  and  3  are  equal. 

The  velocities  of  regions  1  and  h  are  equal  to  zero.  The  resulting  expression 
for  shock  strength  in  the  shock  tube  is: 


p  =  pressure 

7  =  ratio  of  specific  heats  (cp/cy) 

Ms  =  shock  Mach  number 
a  =  speed  of  sound 

The  shock  Mach  number  is  a  function  only  of  stagnant  conditions  existing 
before  diaphragm  rupture.  With  increasing  diaphragm  pressure  ratio  (ph/pi) 
shock  strength  increases  asymptotically  to  an  upper  limit,  determined  primarily 
by  the  speed  of  sound  ratio  (a4/a;]_).  A  low  value  of  7^,  as  well  as  a  high 
value  for  a^/a^,  is  desirable  to  obtain  maximum  shock  strengths.  If  the  initial 
temperature  is  the  same  on  both  sides  of  the  diaphragm,  as  is  usually  the  case 
in  simple  shock  tubes,  the  maximum  attainable  shock  strengths  are  given  in 
Table  3-3.  However,  by  going  to  a  driver  gas  with  a  much  higher  temperature  than 
the  gas  in  the  low-pressure  side,  stronger  shock  waves  can  be  generated,  which  is  sun 
advantage  of  the  arc  driven  shock  tube.  The  shock  Mach  numbers  are  of  course 
only  theoretical  values,  which  are  not  met  in  practice,  but  only  for  Mach 
numbers  of  20  or  greater  can  a  high  degree  of  ionization  be  obtained. 

The  initial  operation  in  a  divergent  channel  with  the  diaphraga  downstream 
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TAB  IE  3-3 

MAXIMUM  SHOCK  STRENGTHS  OBTAINABLE  IN  A  SHOCK  TUBE 


Gas  Combination 

<Vai 

(Pp/Pt )  max 

Ms 

Bmax 

Air -air 

1.00 

44 

6.16 

Helium -air 

2.93 

132 

10.64 

Hydrogen-nitrogen 

3*73 

574 

22.42 

Hydrogen -air 

3-65 

590 

22.5 

of  the  throat,  as  in  the  McDonnell  HIT,  is  similar  to  the  operation  of  a  shock 
tube  of  constant  area,  except  that  the  gas  velocity  in  region  4,  Figure  3-9, 
is  no  longer  zero,  but  equal  to  the  local  speed  of  sound.  With  this  one  change, 
equation  3-2  now  becomes: 


where  is  the  gas  velocity  in  region  4.  If  the  diaphragm  is  located  exactly 
at  the  throat,  or  upstream  of  the  throat,  the  expansion  wave  propagates  into 
essentially  a  stagnant  region  and  initial  operation  is  identical  to  that  of 
the  shock  tube. 

At  diaphragm  rupture,  a  high  velocity  normal,  shock  wave  travels  downstream 
into  the  undisturbed  low  pressure  driven  gas.  Sketches  of  the  pressure  distri¬ 
bution  through  the  shock  excited  region,  as  it  progresses  down  the  nozzle,  are 
shown  in  Figure  3-10.  The  downstream  propagating  shock  is  followed  by  a  region 
of  shock  processed  gas  flow  which  is  terminated  by  the  contact  surface.  Behind 
the  contact  surface  is  a  region  of  steady  flow  which  is  bounded  by  an  upstream 
propagating  disturbance.  Behind  this  disturbance  is  the  quasi-steady  blow¬ 
down  flow.  The  upstream  propagating  disturbance  starts  out  at  the  throat  as  an 
unsteady  expansion  wave,  propagating  at  sonic  velocity,  but  is  swept  downstream 
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into  the  nozzle.  It  decays  in  strength  as  it  progresses  down  the  nozzle. 

Figure  3'10a  and  10b.  At  some  nozzle  station  the  expansion  wave  vanishes 
(Figure  3-10c),  and  a  weak  normal  shock  wave  forms,  propagating  upstream 
(Figure  3-10d).  This  shock  increases  in  strength  as  it  proceeds  downstream. 

When  the  unsteady  expansion  wave  vanishes,  the  pressure  of  the  quasi-steady 
blow  down  flow  equals  the  pressure  of  the  initially  shock  processed  gas.  As  the 
entire  process  travels  further  down  the  nozzle,  the  pressure  at  the  head  of  the 
quasi- steady  blow-down  flow  decreases.  Since  the  pressure  generated  by  the 
initial,  or  downstream,  propagating  shock  is  greater  than  the  blow-down  pressure, 
an  upstream  propagating  shock  wave  is  necessary  to  provide  a  pressure  adjust¬ 
ment.  This  shock  wave  adjusts  the  flow  conditions  between  the  HIT  starting 
process  and  its  quasi-steady  blow-down  flow  out  of  the  driver  section. 

The  initial  shock  strengths  at  diaphragm  rupture  in  the  divergent  channel, 
discussed  above,  have  been  calculated  and  are  shown  in  Figure  3-11  (a-c).  The 
shock  Mach  number  in  a  divergent  channel  is  greater  than  the  shock  Mach  number 
of  a  simple  shock  tube  for  the  same  conditions,  due  to  the  addition  of  the 
velocity  term,  Ul*.  Figure  3-12  shows  this  comparison. 

The  normal  shock  decreases  in  strength  as  it  travels  down  the  divergent 
tube.  In  this  program  three  different  analyses  have  been  used  to  predict  the 
shock  decay,  because  none  seem  to  be  clearly  superior  to  the  other  two.  One 
method  is  based  on  an  approximate  analysis  presented  in  References  3-11,  3-12 
and  3-13-  A  second  approach  uses  the  method  of  characteristics  as  defined  in 
References  3-14-  and  3-15*  The  third  method  is  an  "integral"  approach, 
utilizing  the  simultaneous  solution  of  the  basic  equations  describing  the  wave 
and  flow  phenomena  existing  at  discrete  intervals  along  the  axis  of  a  diver¬ 
gent  channel.  Several  typical  sets  of  conditions  for  a  helium-air  system  will 
be  compared  in  the  following  paragraphs,  utilizing  the  three  approaches  to 
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A 


describe  the  shock  and  the  shock  excited  region. 

Analytical  Approximation  -  Chisnell  has  developed  an  approximate 
expression  for  the  primary  shock  wave,  as  it  passes  through  a  diverging  duct. 
The  equation  is: 

K 


M 

M 


si 

s2 


OO 

2 


The  exponent  Koo  is  a  function  of  the  specific  heat  ratio, 
the  expression: 


('\E)  ('•«=') 


-1 


y 


and  is  given  by 


Shock  wave  attenuation  curves  for  three  typical  values  of  y  are  presented  in 
Figure  3-13.  A  value  for  V  of  approximately  1.22  is  representative  of  con¬ 
ditions  in  shock  heated  air. 

Method  of  Characteristics  -  The  method  of  characteristics,  for  the 
solution  of  the  differential  equations  describing  simple  one-dimensional, 
unsteady  duct  flow,  has  been  programmed  on  an  IBM  7090  digital  computer  to 
arithmetically  generate  the  system  of  characteristic  lines  for  ideal  gas  flow. 
By  plotting  the  generated  points  and  drawing  in  the  network  lines,  variations 
in  properties  throughout  the  field  may  also  be  determined. 

Figure  3-14  presents  the  shock  attenuation  for  several  helium-air 
systems.  The  related  test  times  are  given  in  Table  3-4.  A  portion  of  the 
characteristics  network  at  approximately  23  feet  from  the  diaphragm  for  the 
HIT  (10  degree  conical  nozzle  with  a  1.5  inch  diameter  throat)  is  presented  in 
Figure  3-15 >  fo r  a  8l00°K  and  6200  psia  helium  driver  into  room  temperature 
air  at  5  microns  of  Hg.  A  cross  plot  of  the  conditions  at  the  test  station, 
presented  in  Figure  3-16,  gives  the  variation  of  gas  temperature  with  test 
time  in  the  shock  excited  region.  The  ordinate  length  between  the  shock  and 
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TABLE  3-4 
TEST  TIMES 


the  contact  surface  on  the  characteristics  plot  represents  the  available  time 
for  testing  in  the  shock  excited  region.  Figure  3-17  gives  the  spatial 
position  of  the  shock  and  contact  surface  at  one  initial  condition,  and  presents 
the  temperature  distribution  through  the  shock  excited  region. 

Integral  Approach  -  Rather  than  use  the  differential  equations  themselves 
to  determine  the  flow  properties  of  the  shock  excited  region,  as  was  done  in 
the  previous  section  for  the  characteristics  solution,  the  integral  equations 
describing  the  various  wave  and  flow  phenomena  are  used  here.  Briefly,  at 
selected  points  along  the  divergent  duct,  properties  of  the  blow-down  flow 
region  are  determined  from  the  isentropic  flow  relationships.  The  blow-down 
region  represents  the  region  into  which  either  the  expansion  wave  or  secondary 
shock  wave  propagates.  The  remaining  equations  describing  the  relationships 

across  shocks  or  expansion  waves  are  solved  by  adjusting  the  value  of  the 
initial  shock  strength.  Here  too,  the  solution  has  been  obtained  from  a 
digital  computer  program,  utilizing  the  ideal  gas  relationships.  Unlike 
the  characteristics  solution,  each  region  is  considered  to  have  constant 
properties,  and  the  test  times  must  be  determined  by  first  numerically  inte- 
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grating  the  relationship  between  velocity  and  distance  according  to  the 
expression: 


to  obtain  the  time  for  a  particular  wave  or  the  contact  surface  to  travel 
down  the  tunnel  to  a  specified  station.  Subtracting  the, propagation  time 
of  the  initial  shock  from  the  time  it  takes  the  contact  surface  to  travel  down 
the  tunnel,  results  in  the  available  test  time.  Curves  of  the  initial  shock 
velocity,  interface  velocity,  and  upstream  propagating  wave  velocity  are  pre¬ 
sented  in  Figure  3-18  (a-d).  In  addition,  a  hand  computation  was  performed 
on  the  helium-air  system  cases,  using  real  gas  properties  for  the  shock 
phenomenon  occurring  in  the  air,  using  ideal  gas  properties  for  helium. 

These  results  are  included  on  Figure  3-18  (a-d).  It  is  seen  that  the  solu¬ 
tion  using  real  gas  properties  predicts  lower  velocities,  than  does  the  ideal 
gas  solution.  This  is  to  be  expected,  since  the  ideal  gas  pressure  ratio 
across  a  shock  is  greater  than  the  real  gas  pressure  ratio,  for  the  same  Mach 
number  and  initial  conditions.  In  the  integral  approach  the  pressures  are 
balanced.  Thus,  assuming  the  same  conditions  exist  in  the  blow-down  region 
for  each  case,  the  pressure  ratio  across  the  initial  shock  must  be  the  same, 
regardless  of  ideal  or  real  gas  assumptions.  For  the  same  pressure  ratio  then, 
the  real  gas  assumption  predicts  a  lower  shock  velocity,  and  hence  a  shorter 
nest  time.  Where  exact  knowledge  of  the  test  time  is  not  required,  the 
assumption  of  ideal  gas  properties  is  satisfactory  for  determining  shock 
decay.  Test  times  between  the  passing  of  the  initial  shock  and  the  contact 
surface  are  compared  in  Table  3_i+  with  those  computed  by  other  methods. 

Electrical  Properties  Behind  a  Moving  Shock  -  The  electrical  properties 
existing  in  the  shock  excited  region  behind  the  initial  shock  can  be  pre- 
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dieted  from  use  of  hypersonic  gas  dynamic  charts  for  air,  such  as  are  pre¬ 
sented  hy  Feldman  (Reference  3-16)  and  by  Ziemer  (Reference  3-17)*  Equili¬ 
brium  densities  and  temperatures  behind  a  moving  normal  shock  are  presented 
in  these  references,  and  are  included  here  in  Figure  3-19  as  functions  of  the 
shock  velocity  and  the  density  ratio  in  the  evacuated  section  before  the 
shock  is  initiated.  The  equilibrium  properties  of  electron  concentration 
behind  the  shock  and  the  electron  collision  frequency  are  presented  in 
Figures  3-20  and  3-21,  and  have  been  determined  from  the  general  curves  of 
equilibrium  properties  at  the  specific  conditions  behind  the  shock,  which  are 
shown  in  Figure  3-19* 

3. 4. 2. 6  Blow-Down  Flow 

General  Description.  -  In  an  arc  driven  Hypervelocity  Impulse  Tunnel,  after 
decay  of  the  initial  shock  structure,  there  follows  a  quasi- steady  flow  from 
the  arc  chamber,  which  is  referred  to  as  the  blow-down  flow.  The  blow-down 
flow  at  the  test  station  is  characterized  by  a  high  velocity  and  a  relatively 
low  temperature ,  pressure  and  density.  Test  periods  in  this  flow  regime  are 
relatively  long,  being  of  the  order  of  100  milliseconds.  The  blow-down  flow 
is  normally  used  for  aerodynamic  testing.  During  the  period  of  blow-down  flow, 
pressure  and  temperature  at  the  test  station  varies  slowly,  in  response  to  the 
decrease  in  pressure  and  temperature  within  the  arc  chamber.  In  the  following 
discussion  flow  properties  characteristic  of  the  blow-down  regime  are  examined, 
with  particular  attention  given  to  the  prediction  of  electron  concentration 
and  other  plasma  properties  at  the  test  station. 

After  the  electrical  discharge,  the  gas  supply  in  the  arc  chamber  con¬ 
sists  of  a  complex  mixture  of  molecules,  atoms,  ions,  and  electrons  which 
are  assumed  to  be  in  chemical  equilibrium  at  the  arc  chamber  temperature 
and  pressure.  As  the  gas  discharges  from  the  arc  chamber  and  expands  through 
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the  divergent  nozzle,  the  thermodynamic  state  of  the  gas  is  altered  in  accord¬ 
ance  with  the  laws  of  compressible  flow  and  chemical  kinetics.  When  a  ga$  of 
constant  specific  heat  and  chemical  composition  flows  through  a  divergent 
nozzle,  the  response  of  the  thermodynamic  properties  to  the  changing  area 
ratio  is  rapid  with  respect  to  the  rate  of  flow  of  the  gas.  For  this  case, 
on  a  microscopic  scale,  the  rate  of  approach  to  an  equilibrium  temperature 
and  pressure  is  governed  by  the  rate  at  which  a  steady  state  Maxwellian 
velocity  distribution  can  be  achieved.  For  gases  consisting  of  species  of 
nearly  equal  mass,  the  Maxwellian  velocity  distribution  is  reached  at  a  time 
on  the  order  of  the  time  between  molecular  collisions,  or  the  reciprocal  of 
the  collision  frequency.  When  electrons  are  present,  the  masses  of  the 
particles  are  no  longer  nearly  equal,  and,  because  electrons  normally  undergo 
elastic  collisions  with  particles  of  greater  mass,  only  a  small  fraction  of 
the  kinetic  energy  of  the  electron  is  transferred  to  the  heavier  particle  upon 
colliding.  It  is  possible,  therefore,  for  electrons  to  maintain  a  kinetic 
energy,  or  electron  temperature,  which  is  not  in  equilibrium  with  that  of  the 
bulk  of  the  gas  for  a  significant  period.  If  chemical  reactions  and  variable 
specific  heats  are  allowed,  the  rate  of  approach  to  chemical  equilibrium  is 
governed  by  the  specific  reaction  rate  constants  of  the  participating  chemical 
reactions,,  and  by  the  relaxation  times  for  rearrangement  of  internal  molecular 
degrees  of  freedom. 

Exact  calculations  for  the  hypersonic  nozzle  flow  of  complex  gas  mixture s 
requires  the  simultaneous  solution  of  a  number  of  rate  equations,  one  for  each 
process,  several  of  which  may  be  highly  coupled  in  a  real  case.  The  exact 
solution  is  further  complicated  by  a  lack  of  knowledge  of  the  rate  constants 
for  many  high  temperature  processes.  For  engineering  purposes  it  is  conmon, 
therefore,  to  define  two  limiting  cases:  (l)  frozen  flow,  and  (2)  equilibrium 
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flow,  corresponding,  respectively,  to  the  assumptions  of  constant  chemical 
composition  (infinitely  slow  rate  constants)  and  equilibrium  composition 
(infinitely  fast  rate  constants).  A  comparison  between  the  frozen  flow  calcu¬ 
lations  and  the  more  appropriate  non-equilibrium  flow  analysis  is  presented 
in  the  following  paragraphs. 

Frozen  Flow  Analysis.  -  Previous  theoretical  studies  of  the  nature  of 
hypersonic  nozzle  flow  provide  the  following  description:  As  the  gases  expand 
from  the  nozzle  throat  into  the  divergent  channel  there  is  an  initial  period 
in  which  the  gases  are  in  chemical  equilibrium.  The  dissociation  and  recombin¬ 
ation  reactions  between  atomic  and  diatomic  species  requires  the  presence  of  a 
third,  inert  particle;  thus,  these  reactions  are  density  dependent.  As  the 
gases  expand  and  the  density  decreases,  there  is  a  short  transition  period, 
which  is  followed  by  a  third  region  characterized  by  constant  composition, 
i.e.,  frozen  flow.  Once  frozen,  the  flow  is  not  observed  to  return  to  an 
equilibrium  composition  at  a  more  distant  section  of  the  nozzle,  although  this 
possibility  exists  by  the  nature  of  the  chemical  rate  equations.  Freezing 
normally  occurs  a  short  distance  from  the  nozzle  throat.  As  an  approximation, 
for  engineering  calculations,  freezing  may  be  assumed  to  occur  at  the  nozzle 
throat  or  in  the  upstream  reservoir,  i.e.,  in  the  arc  chamber,  where  the 
equilibrium  composition  is  readily  obtained.  For  the  calculations  presented 
here,  freezing  is  assumed  to  occur  at  the  equilibrium  arc  chamber  composition. 

Although  the  composition  may  be  frozen,  it  is  still  necessary  to  specify 
the  specific  heats  and  the  specific  heat  ratio,  7  ,  at  every  point  in  the 
nozzle.  Specific  heats  are  functions  of  temperature  and  may  also  be  con¬ 
sidered  to  be  either  (l)  frozen,  (2)  in  equilibrium,  or  (3)  in  a  non-equilibrium 
state,  characterized  by  a  relaxation  time  for  the  rearrangement  of  the  interned 
molecular  degrees  of  freedom.  Such  relaxation  times  are  known  to  be  longer  at 
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high  temperatures,  where  they  approach  in  magnitude  the  times  characteristic 
of  the  chemical  reactions.  The  assumption  of  equilibrium  specific  heats,  or 
A  frozen  vibrational  energy,  has  been  made  for  the  calculations  presented  in  this 
section.  It  can  be  shown  that  the  predicted  density,  and,  therefore,  the 
electron  concentration,  will  not  be  appreciably  affected  by  this  assumption. 

Heims  (Reference  3-18)  and  Inger  (Reference  3-19) ,  among  others,  have 
discussed  the  calculations  for  expanding  gases  in  frozen,  dissociated  flow. 
Following  Heims,  the  local  isentropic  exponent  for  adiabatic  flow  of  a  perfect 
gas  may  be  obtained  as  follows  for  an  adiabatic  expansion: 


The  above  equations  describe  the  frozen  expansion  of  a  dissociated  gas,  the 
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contributions  from  ionic  species  and  electrons  being  negligible. 

In  Figures  3-22  and  3-23  are  presented  charts  of  temperature  and  electron 
concentration  at  a  36  inch  diameter  test  section  in  the  McDonnell  HIT.  These 
figures  describe  test  station  properties  for  the  assumed  frozen  flow  of  nitrogen,, 
initially  at  equilibrium  in  the  arc  chamber.  The  charts  apply  to  an  area  ratio 
of  8100,  corresponding  to  a  0.4  inch  diameter  nozzle  throat.  In  the  Alfven 
wave  experiments  a  solenoid  for  the  magnetic  field  B0,  inserted  into  the  diverg¬ 
ing  nozzle,  had  an  internal  diameter  of  36  inches,  which  then  became  the 
geometric  area  available  for  flow.  If  the  charts  are  to  be  used  as  labeled, 
it  is  implicitly  assumed  that  the  composition  at  which  freezing  occurs  does  not 
differ  appreciably  from  arc  chamber  conditions.  It  may  be  more  appropriate  to 
assume  that  the  flow  to  the  throat  is  in  chemical  equilibrium;  for  this  case 
chemical  composition,  temperature  and  pressure  at  the  throat  may  be  estimated 
from  the  previous  data  of  section  3* 4.2.4,  and  these  values  may  be  substituted 
for  the  corresponding  arc  chamber  values  on  the  figures. 

The  effect  of  a  divergent  nozzle  on  test  station  properties  for  nitrogen 
flow  is  illustrated  in  Figures  3-24  and  3-25;  where  temperature  and  electron 

concentrations  are  shown  as  functions  of  area  ratio  for  specific  values  of 
temperature  and  pressure  in  the  arc  chamber. 

Equilibrium  molar  compositions  at  arc  chamber  conditions  have  been  com¬ 
puted  from  data  on  equilibrium  constants,  supplemented  with  the  published  re¬ 
sults  for  nitrogen  of  Treanor  and  Logan  (Reference  3-20 )  and  the  recent  high 
temperature  calculations  of  Drellishak  (Reference  3-21). 

Non-Equilibrium  Flow  Analysis.  -  A  more  nearly  correct  computation  for 
the  flow  conditions  in  a  hypersonic  nozzle  may  be  obtained  from  a  simultaneous 
solution  of  the  flow  equations,  together  with  the  rate  equations  for  the  chem¬ 
ical  reactions  and  material,  balances  for  the  severed,  species.  Bray  (Reference 
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3-22),  Emanuel  and  Vincent!  (Reference  3-23),  Westenberg  and  Favin  (Reference 
3-24),  and  Eschenroeder  (References  3-25  and  3-26)  have  described,  in  general, 
the  approach  to  non-equilibrium  flow  calculations.  The  analysis  presented 
here  follows  more  closely  the  study  of  Westenberg  and  Favin.,  It  is  assumed 
that  the  internal  molecular  degrees  of  freedom  are  in  equilibrium  at  the  local 
temperature ,  and  that  the  electron  temperature  is  everywhere  equal  to  the 
local  flow  temperature.  Almost  no  information  is  available  concerning  the 
effect  of  electron  temperature  on  the  kinetics  of  the  recombination  reactions. 

A  stepwise  iterative  calculation  is  required.  Let  A  be  the  cross-sectional 
area  of  the  nozzle  at  any  point  z  along  the  nozzle  axis,  and  let  A*  be  the 
cross-sectional  area  at  the  nozzle  throat.  Then  the  area  ratio,  if  =  A/A*, 
is  a  function  of  z.  The  flow  through  the  divergent  nozzle  is  described  by  the 
following  equations: 

(l)  Equation  of  continuity: 

m  «.  p  v  77  A* 


(2)  Equation  of  motion: 


dp 

dz 


0 


(3)  Equation  of  energy  balance: 


dh  V  dv 

dz  J  dz 


(4)  Equation  of  state: 


P  -  P  RTJ  (2  F:) 

i 


where  the  subscript  i  refers  to  the  ith  species  and  the  other  symbols  are 
defined  as  follows: 
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m  =  mass  flow  rate  (pa/ sec) 

P  »  gas  density 

v  =  velocity  of  an  elementary  volume 
h  =  enthalpy  (cal/gm) 

J  *  mechanical  equivalent  of  heat  (erg/ cal) 

R  =  Universal  gas  constant 
T  =  gas  temperature  (°K) 

F  =  component  concentration  (gmols./gm) 

If  H  denotes  the  enthalpy  per  mole,  then: 

hj  =  Hj  Fj 


For  a  perfect  gas,  enthalpy  is  a  function  of  temperature  only,  hence: 


dHi  dHi  dT  dT 

dz  dT  d  z  dz 

The  above  set  of  equations  may  be  solved  to  yield  the  derivatives  of  tempera¬ 
ture  and  density  with  respect  to  axial  distance,  z: 


JY2 


+  TX 


dP 
d  z 


il 

dz  =  1  -  JY2  (FjCp.) 

[£_„] 


p_ 

T)  dz 


where : 


and: 


n  d,  1 

'  [r,  dz  (2  Fj)  dz  J 


~  1  T  p2  A2*  n2' 

_RJ  (2  Fj)  '  m2 


The  chemical  composition  at  each  point  z  is  obtained  from  the  appropriate 
forms  of  the  chemical  rate  equations  and  component  material  balances.  In 
general,  for  the  J-th  reaction: 
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R  K' 1  P 

where  the  summations  are  taken  over  the  components  of  the  reactants  and 
products,  respectively.  The  rate  of  reaction  may  be  expressed  by: 


The  indicated  products  are  to  be  taken  over  all  components  which  appear  as 
reactants  of  products  in  a  given  reaction.  Here  kj  is  the  reaction  rate  constant 
and  Kj  is  the  equilibrium  constant.  The  component  material  balances  are  of  the 
form: 

X  n.  :  F;  =  constant 

I  I 


where  ne,  i  denotes  the  number  of  atoms  of  element  e  in  component  i.  For  the 
nitrogen  flow  calculations  the  elemental  constituents  can  be  considered  to  be 
the  singly  charged  nitrogen  atom  sued  the  electron.  In  this  case  reactions  may 
be  selected,  such  that  the  above  set  of  equations  can  be  solved  directly  by  a 
finite  difference  technique.  In  other  cases  simultaneous  solutions  may  be  re¬ 
quired.  For  each  component  described  uniquely  by  one  of  the  chemical  reactions: 


and: 


dFi  77  A* 
dz  =  IT)  r' 


(Fi^z  +  Az  ~  (Fi^z  + 


Az 


For  the  remaining  components,  compositions  at  the  end  of  an  interval,  (Fi)z  +^z 
are  determined  from  the  material  balances.  Average  densities  and  temperatures 
over  the  interval,  obtained  by  iteration,  are  used  to  calculate  the  reaction 
rates. 

Calculations  for  the  hypersonic  flow  of  ionized  nitrogen  in  a  divergent 
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nozzle  will  be  now  described.  The  species  present  in  the  arc  chamber  are:  N2, 

t  4*  -f* 

N2,  N  ,  N  and  electrons,  e'"..  An  associate  species,  may  become  significant 
at  high  pressures,  but  will  not  be  important  here.  Because  of  the  variety  of 
possible  recombination  mechanisms,  a  large  number  of  reactions  between  the 
five  species  may  be  described. 

The  recombination  reactions  are  exothermic  and  cannot  occur  unless  some 
means  is  provided  for  the  removal  of  the  excess  energy  of  the  reaction  products. 
For  the  recombination  of  ions  and  electrons  energy  may  be  removed  (a)  by 
radiation,  (b)  by  excitation  of  at  least  two  electrons,  (c)  by  energy  of 
molecular  dissociation,  or  (d)  as  kinetic  energy  imparted  to  a  third  body. 

The  corresponding  recombination  reactions,  illustrated  for  the  species  of 
ionized  nitrogen,  are: 

(a)  Radiative  recombination: 

N++e  =  N  +  h  v 
N  j  +  e  =  N2  +  h  v 

where  represents  the  energy  quantum. 

(b)  Dielectronic  recombination: 

N+  +  e—  =  N** 

where  N**  represents  a  doubly  excited  atom. 

(c)  Dissociative  recombination: 

N  2  +  ?  =  N  +  N 

(d)  Three  body  recombination: 

N++e  +M=N  +  M 
N-j  +  e  +  M  =  N2  +  M 

where  M  represents  any  third  species.  In  many  systems  electron  loss  may  also 
occur  through  electron  attachment  to  neutral  species,  but  this  does  not  happen 
with  nitrogen  and  the  rare  gases.  Recombination  mechanisms  have  been  discussed 
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by  a  number  of  authors,  of  which  the  work  of  Brown  (Reference  3“ 27 )  suad  Massey 
(Reference  3-28)  will  be  cited,  and  by  Bialecke  and  Dougal  (Reference  3-29) , 
with  particular  reference  to  the  nitrogen  system. 

It  is  conventional  to  define  an  ionic  recombination  coefficient, , 
regardless  of  the  mechanism,  by: 

dne 

-ir“a'n»ni 

where  rig  and  n^  denote  number  densities  of  electrons  and  ions,  respectively. 

Probable  maximum  values  of  the  recombination  coefficients  for  radiative  and 
dielectronic  recombination  are  of  the  order  of  lO"-*-1  and  10"12  cm^/ sec. , 
respectively.  Dielectronic  recombination  can  usually  be  expected  to  be 
negligible  in  relation  to  the  other  mechanisms,  and  will  not  be  considered 
further  here. 

For  a  monatomic  ion  such  as  N+,  only  radiative  and  three  body  collision 
processes  need  be  considered.  The  relative  importance  of  the  two  mechanisms 
will  depend  upon  the  pressure  and,  to  a  lesser  extent,  upon  the  nature  of  the 
third  body;  third  body  recombination  will  be  dominant  at  the  higher  pressures. 

For  a  molecular  ion  such  as  dissociative  recombination  can  be  expected  to 
be  dominant,  although  three  body  recombination  may  be  a  significant  factor  at 
the  higher  pressures. 

Massey  (Reference  3-28)  has  discussed  the  theoretical  derivation  by 
J.J.  Thomson  for  the  three  body  recombination  coefficient  between  positive  and 
negative  ions.  If  appropriate  simplifications  are  made  for  electron- ion  recom¬ 
bination,  the  recombination  coefficient  for  a  three-body  process  is  given  by: 

A  .  ..3, 
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r0  -  2  e2/3  k  T 

e  =  the  charge  of  an  electron 
A  m  =  mass  of  an  electron,  gms 

Mj_  =  mass  of  an  ion,  gms 

L  =  mean  free  path  of  elastic  electron-neutral  collisions 
Thenean  free  path  L  is  given  by 


inn 


where  r2  =  molecular  radius  of  the  neutral  species 

n2  =  particle  concentration  of  a  neutral  species 
For  the  three  body  recombination  of  molecular  nitrogen  ions  at  typical  arc 
chamber  conditions  (P  =  2000  psi,  T  =  10,000  K),  it  is  estimated  from  equation 
for  ar  that  aT  is  approximately  5  x  10  cm"5/ sec.  Measured  recombination 

coefficients  for  electron  recombination  are  found  to  be  several  orders  of 
magnitude  larger,  so  that  for  this  reaction  the  dissociative  mechanism  will 
predominate.  For  the  recombination  of  N*  ions  the  situation  is  less  clear;  for 
this  case  the  magnitudes  of  the  recombination  coefficients  for  the  radiative 
and  three  body  mechanisms  are  more  nearly  equal.  Massey  (Reference  3-28)  tabu¬ 
lates  radiative  recombination  coefficients  for  atomic  oxygen  ions.  For  recom¬ 
bination  with  high  energy  electrons,  the  recombination  coefficient  is  of  the 
order  of  0.3  x  10 cm^/sec.  The  three  body  recombination  coefficient  will 
not  differ  greatly  from  the  value  computed  above  for  N^.  On  this  basis  the 
three  body  recombination  mechanism  will  predominate,  since  it  is  well-known 
that  the  recombination  of  nitrogen  atoms  to  molecular  nitrogen  takes  place 
through  a  three-body  mechanism. 

On  the  basis  of  the  preceding  discussion  the  following  reactions  have 
been  selected  to  describe  the  reacting  mixture. 
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N  +  N  +  M  m  N2  +  M 
Nj  +  ©  sa  N  +  N 
N++e  +M=N  +  M 

A  number  of  authors  have  studied  the  rate  of  dissociation  of  molecular 
nitrogen.  Wray  (Reference  3-30)  recommends  the  data  of  Byron,  which  is  used 
here.  For  the  rate  of  the  reverse  of  the  first  reaction: 

k '  =  {4.2  x  1012)  /2  (D/RT)  exp  (-D/RT) 

where  D  *  224,900  cal/gmol  is  the  dissociation  energy,  and  the  prime  refers  to 
the  reverse  reaction  above.  This  rate  expression  applies  to  reactions 
in  which  molecular  nitrogen  is  the  third  body. 

Kuhns  (Reference  3-31)  has  studied  the  dissociative  recombination  of  Hg 
ions  and  reports  a  value  of  ar  =  1*9  x  10-a  cnr/  sec  at  room  temperature  and 
10  mm  Hg.  Biale eke  and  Dougal  (Reference  3-29)  report  ar  «  1.6  x  10"^  cm^/sec 
at  300 °K  and  at  a  particle  density  of  8  x  lO1^  cm" They  have  also  found  that 
the  recombination  coefficient  decreases  with  increasing  temperature  and 
increases  with  increasing  density,  but  they  do  not  present  data  for  the  tempera¬ 
tures  and  pressures  of  interest  here.  The  pressure  dependence  in  the  observed 
data  is  thought  to  arise  from  the  presence  of  three-body  recombination.  Bates 
(Reference  3-32)  has  presented  a  theoretical  discussion  of  the  dissociative 
recombination  coefficient  and  has  concluded  that  it  will  be  temperature  depen¬ 
dent  by  a  factor  varying  from  T-1/2  to  T "3/2,  depending  upon  the  electron  energy 
as  described  by  initial  and  final  wave  functions,  which  are  in  general  unknown. 

In  view  of  these  uncertainties,  a  constant  value  of  «r  =  1,5  x  10”^  cm^/sec  has 
been  selected  for  these  preliminary  calculations. 

For  the  three  body  recombination  of  M+  ions  the  equation  of  J.J.  Thomson, 
referred  to  earlier,  has  been  used  to  compute  the  recombination  coefficient. 

This  equation  is  subject  to  a  high  pressure  limit,  but  it  can  be  shown  that 
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this  limit  is  not  exceeded  at  arc  chamber  pressures. 

Reaction  rate  constants  in  the  form  suitable  for  equation  3*^-3  have  been 
computed  from  the  recombination  coefficients  cited  above.  Equilibrium  con¬ 
stants  for  the  previously  selected  reactions  have  been  obtained  by 
rearrangement  of  the  tabulated  data  of  Treanor  and  Logan  (Reference  3-20) • 

Enthalpies  and  heat  capacities  have  been  computed  from  partition  functions,  also 
tabulated  by  Treanor  and  Logan.  Together  with  the  above  data,  the  equations 
described  in  this  section  have  been  programmed  for  solution  on  an  IBM  7090 
digital  computer. 

The  restats  of  some  preliminary  non-equilibrium  flow  calculations  for 
ionized  nitrogen  in  the  divergent  section  of  a  hypersonic  nozzle  sure  presented 
in  Figure  3-26.  These  calculations  were  initiated  at  the  nozzle  throat  with 
an  equilibrium  gas  composition  corresponding  to  a  temperature  of  8000°K,  and 
a  density  of  1.1845  x  10“2  gms/cm^,  approximately  10  times  the  atmospheric 
density.  Initial  equilibrium  mole  fractions  of  each  component  are  tabulated 
below: 


Component 

.  Mole  Fraction 

n2 

8.437  x  10-1 

N2 

5.966  x  10"5 

N+ 

3.134  x  10"5 

N 

1.561  x  10'1 

e" 

9.099  x  10-5 

Because  of  the  extremely  rapid  reaction  rates,  an  unusually  small  step  size, 

-8 

initially  of  the  order  of  10"  cm.,  was  required  for  the  iterations.  A 
nozzle  with  a  10  degree  total  divergent  angle  has  been  assumed  for  these  calcu¬ 
lations,  corresponding  to  the  construction  of  the  McDonnell  HIT.  The  total 
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nozzle  length  to  the  50  inch  diameter  test  station  is  approximately  J20  cm. 

The  preliminary  calculations  presented  here,  however,  do  not  take  into  account 
the  reduction  in  effective  area,  which  results  from  the  boundary  layer  in  the 
nozzle.  It  is  reasonable,  therefore,  to  compare  these  non- equilibrium  calcu¬ 
lations  with  the  previously  presented  frozen  flow  calculations,  at  equivalent 
effective  area  ratios.  With  the  magnetic  coil  in  place,  the  area  ratio  was 
estimated  previously  to  be  8100;  this  corresponds  to  a  distance  of  approxi¬ 
mately  360  cm.  in  the  present  calculations. 

From  Figure  3-26  it  is  seen  that  the  rate  of  recombination  of  Ng  ions  is 
very  rapid,  and  that  this  species  is  essentially  removed  from  the  gas  stream 
within  a  distance  of  6.h  cm.  from  the  throat.  On  the  other  hand,  the  loss  by 
recombination  of  IT*"  ions  is  negligible,  amounting  to  less  than  0.8  percent  of 
the  original  W*  ions  present.  This  reaction  is  essentially  frozen  beyond  the 
first  6  cm.  The  number  of  electrons  remaining  in  the  flow  is  equal  to  the 
number  of  N+  ions,  and  also  remains  constant  beyond  the  first  6  cm.  The 
reaction  leading  to  the  recombination  of  N  atoms  to  molecular  nitrogen  Np* 
becomes  frozen  at  approximately  20  cm.  from  the  throat;  approximately  85  per¬ 
cent  of  the  original  N  atoms  have  been  lost  by  recombination  at  this  point. 

The  reactions  leading  to  the  recombination  of  both  N+  ions  and  N  atoms  are 
three  body  reactions,  and  the  mass  rates  of  these  reactions  vary  with  the 
square  of  the  density.  Because  of  the  rapid  decrease  in  density  in  the  nozzle, 
these  reactions  become  "frozen"  and  do  not  again  become  significant.  Hie 
recombination  of  Ng  ions,  on  the  other  hand,  has  been  assumed  to  take  place 
through  the  dissociative  mechanism.  The  mass  rate  of  this  reaction  is 
relatively  rapid  and  varies  directly  with  the  density.  It  is  not  surprising, 
therefore,  that  the  recombination  of  ions  proceeds  to  completion. 

For  comparison,  calculations  for  the  same  initial  conditions  have  been 
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carried  out  with  the  frozen  flow  assumption,  i.e.,  all  reaction  rates  have 
been  set  equal  to  zero.  In  these  calculations,  specific  heats  have  been  com¬ 
puted  from  the  program,  and  are  at  an  equilibrium  value  at  the  local  tempera¬ 
ture.  Because  the  energy  in  vibration  is  small  compared  to  the  energy  in 
dissociation  and  ionization,  differences  arising  from  this  assumption  are 
expected  to  be  small.  A  comparison  of  the  frozen  flow  and  non-equilibrium 
calculations  for  the  same  initial  conditions  is  presented  in  Figure  3-27 • 

It  is  seen  that  density  and  velocity  at  the  test  station  do  not  differ 
significantly  for  the  two  cases.  There  is  a  significant  difference  in  the 
temperature,  however,  which  is  higher  for  the  non-equilibrium  calcinations. 

The  purpose  of  the  non-equilibrium  flow  calculations  is  twofold: 
first,  to  provide  a  basis  for  comparison  with  the  frozen  flow  calculations, 
and  second,  to  provide  an  insight  into  the  dynamic  behavior  of  the  component 
species  within  the  nozzle.  Using  the  density  and  composition  obtained  from 
the  non-equilibrium  flow  calculations  at  the  area  ratio  of  8100,  the  electron 
number  density  is  estimated  to  be  6  x  10  cm  .  For  ionized  nitrogen  in 
equilibrium  at  8000°K  and  a  density  ratio  of  10  relative  to  standard  density, 
the  pressure  is  approximately  4400  psi.  From  Figure  3-23  the  electron  number 

density  predicted  from  the  frozen  flow  calculations  is  found  to  be  approxi- 
12  -3 

mately  1.5  x  10  cm  ,  which  is  in  good  agreement  with  the  non-equilibrium 

calculations,  if  allowance  is  made  for  the  fraction  of  electrons  lost  by 

+  _ 
recombination  with  N2  ions.  At  typical  arc  chamber  conditions  of  12,000°K 

and  2000  psi  for  the  experimental  runs,  the  number  of  atomic  nitrogen  ions 
exceeds  the  number  of  molecular  nitrogen  ions  at  equilibrium;  hence,  pro¬ 
portionately  fewer  electrons  would  be  lost  by  recombination  with  Ng,  thus, 
the  electron  concentration  at  the  test  station  is  expected  to  differ  by  less 
than  a  factor  of  two  from  that  obtained  by  the  frozen  flow  calculations  of 
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CONDITIONS 

NOZZLE  TOTAL  ANGLEt  10° 
EQUILIBRIUM  SPECIFIC  HEATS. 
INITIAL  PROPERTIES  AT  THROATi 
T*  =  8000°K 

p*  *  1.1845  x  10~^  gm/cm^ 

V*=  1.84  x  105  cm/sec 
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section  3.I4-.2.6.  For  the  MHD  experiments  with  which  this  report  is  concerned, 

a  loss  of  electrons  of  this  magnitude  is  not  important.  Furthermore,  the 

+ 

elimination  of  the  molecular  nitrogen  ion,  N2>  results  in  a  significant 
reduction  in  the  complexity  of  the  plasma  wave  equations  and  simplifies  the 
experimental  analysis  of  the  observed  wave  patterns. 

The  non- equilibrium  flow  calculations  presented  here  must  be  considered 
to  be  only  approximate,  because  of  a  lack  of  knowledge  of  the  true  values  of 
recombination  rates  at  the  temperatures  and  pressures  of  importance  to  this 
study.  Calculations  have  been  carried  out  for  an  initial  temperature  of 
8000°K,  rather  than  12,000°K,  because  of  the  lack  of  complete  published  thermo¬ 
dynamic  data  at  the  higher  temperature.  This  data  is  expected  to  be  available 
soon.  A  more  nearly  exact  calculation  would  take  into  account  non-equilibrium 
specific  heats,  as  well  as  the  calculation  of  electron  temperature,  and  the 
effect  of  electron  temperature  on  recombination  rates.  Nevertheless,  it  is 
believed  that  these  calculations  present  a  reasonably  close  approximation 
to  the  true  behavior  of  ionized  nitrogen  in  a  hypersonic  nozzle. 

Notes  on  Blow  Down  Flow  Analysis.  -  The  basic  analysis  presented  for  the 
blow  down  flow  has  been  performed  for  the  geometric  area  ratio  existing  between 
the  test  section  and  the  nozzle  throat  (A/A*  =  8100).  A  few  additional  figures 
have  been  presented  to  indicate  how  certain  conditions  vary  with  area  ratio. 
Although,  it  ia  well  known,  that  a  boundary  layer  does  indeed  exist  inside  the 
tunnel,  effectively  reducing  the  core  size,  a  lack  of  data  on  core  size  made 
mandatory  certain  as sumptions  in  predicting  the  test  section  conditions. 
Certainly,  the  assumption  of  effective  area  ratio  equal  to  the  geometric  area 
ratio  predicts  lower  temperatures  and  electron  concentrations  than  might  be 
expected,  thus  providing  a  lower  limit  on  the  determination  of  these  values. 

On  the  other  hand,  a  minimum  of  data  also  means  approximate  methods  must  be 
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devised  to  provide  answers. 

As  will  be  pointed  out  in  the  sections  on  instrumentation  and  on  data 
reduction,  (Section  3*2.3  Vol.  Ill)  only  two  types  of  measurements  were 
available  upon  which  to  base  any  type  of  flow  field  analysis.  During  the 
electrical  exciter  experiments,  a  smear  photograph  captured  luminous  particle 
traces,  yielding  a  gas  velocity.  During  subsequent  runs,  a  pitot- static 
pressure  pickup  measured  the  stagnation  pressure  behind  the  normal  shock  in 
the  test  section.  The  photographically  determined  velocity  was  used  to  deter¬ 
mine  an  area  ratio  from  which  a  density  and  a  temperature  were  computed.  This 
computation  was  based  upon  the  frozen  flow  analysis  computer  analysis. 

The  usual  expressions  to  interpret  the  pressure  measurements  are  no 
longer  valid  when  the  free  stream  is  dissociated  and  ionized,  as  is  the  case 
of  the  blow  down  flow.  The  simplest  solution  has  been  to  ignore  this  pitfall, 
and  to  assume  an  effective  ratio  of  specific  heats,  Y  ,  for  use  in  the  normal 
expressions.  This  effective  Y  has  no  bearing  on  the  effective  Y  used  pre¬ 
viously  in  the  discussion  of  the  expansion  process  and  frozen  flow,  since  the 
nozzle  expansion  normally  is  based  on  state  conditions  initially  at  the 
throat,  rather  than  the  free  stream  conditions  existing  at  the  test  section. 
Referring  again  to  the  computer  results  for  different  area  ratios,  a  Mach 
number  and  pressure  was  determined  at  a  given  value  of  area  ratio.  Assuming 
an  effective  Y  of  1.3,  a  calculated  stagnation  pressure  behind  the  shock 
was  determined  and  compared  with  the  measured  value.  A  new  value  of  area  ratio 
was  then  assumed  and  the  process  repeated  until  the  calculated  and  measured 
pressure  agreed.  The  final  area  ratio  was  considered  the  effective  area  ratio 
at  which  the  other  conditions  could  be  determined.  Para.  3*2-3,  Vol.  Ill,  has 
been  prepared  in  this  manner.  Unavoidably,  the  effective  area  ratio  upon 
which  this  table  is  based  seems  much  too  great  in  that  the  predicted  tempera- 
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ture  appears  unreasonably  low.  Although  it  was  attempted,  a  different  value 
for  the  effective  y  might  have  predicted  a  more  reasonable  temperature.  At 
best  then,  this  computation  might  replace  the  assumption  of  A/A*  =  8100  as 
the  lower  limit  on  test  section  conditions,  nevertheless,  these  computations 
are  based  on  an  essentially  sound  theory  for  frozen  flow,  with  additional 
confirmation  coming  from  the  measured  electron  concentrations  at  the  test 
section. 


Discharge  Tube.  -  Two  small  plasma  facilities  acquired  during  this  program 
complemented  the  HIT.  The  addition  of  these  facilities  with  much  greater 
recycle  times  became  desirable,  as  the  number  of  necessary  development  and 
exploratory  experiments  increased.  Then  too,  the  HIT  could  not  attain  the 
high  ionization  percentages,  which  more  easily  would  support  Alfven  wave 
propagation.  The  electromagnetic  shock  tube  generates  a  moving  plug  of  plasma 
by  striking  an  arc  coaxially,  at  one  end  of  a  pyrex  glass  tube.  The  basic 
processes  involved  in  producing  the  plasma  are  largely  unknown,  however,  the 
resulting  characteristics  of  high  temperature  and  high  ionization  percentage 
can  be  easily  measured.  The  second  facility,  an  arc  discharge  tube,  produces 
ionization  by  striking  an  arc  along  the  length  of  a  pyrex  tube.  While  the 
arc  is  still  burning,  the  experiments  are  performed.  As  in  the  electromagnetic 
shock  tube,  the  basic  processes  are  not  understood,  nevertheless,  the  charac¬ 
teristics  of  the  resulting  plasma  can  be  measured. 

3. 4. 3.1  Electromagnetic  Shock  Tubes.  -  The  electromagnetic  (EM)  shock 
tube  was  originally  developed  to  surmount  the  limits  of  attainable  shock 
strength,  available  from  conventional  diaphragm  type  tubes.  Basically,  the 
operation  of  the  EM  shock  tube  involves  rapidly  discharging  a  capacitor 
through  gas  in  a  chamber,  heating  and  ionizing  the  gas,  and  accelerating  it 


JMCi 


L 


144 


IV1MD  Wave  Investigation  30RI November  tie® 

1 1  -  theory 

out  of  the  chamber  into  a  region  of  cold  gas.  The  acceleration  takes  place, 
at  least  in  part,  because  of  the  interaction  between  currents  in  the  discharge 
plasma  and  the  magnetic  field  of  the  return  currents,  through  the  external 
circuits,  supplying  the  discharge. 

If  the  "plug"  of  plasma  were  regarded  as  a  moving  piston  compressing  the 
cold  gas  ahead  of  it,  the  existence  of  a  shock  moving  ahead  of  the  driver  gas 
would  be  predicted.  Something  resembling  a  shock  has  indeed  been  observed  by 
several  investigators  under  certain  conditions  of  initial  gas  pressure.  Con¬ 
ditions  behind  the  shock  would  depend  on  the  state  of  the  gas  ahead  of  it  and 
could  be  calculated  from  the  conservation  laws,  in  combination  with  the  equations 
of  state  of  the  gas.  The  equations  would  also  have  to  take  into  account  the 
fact  that  gas  ahead  of  the  shock  may  be  excited  by  ultraviolet  or  x-radiation 
from  the  discharge.  Theories  exist  which  predict  fairly  accurately  the 
velocity  of  the  plasma  as  it  leaves  the  discharge  region. 

The  region  of  shock  heated  gas  may  not  endure  for  long,  even  in  those 
cases  where  a  well  defined  shock  is  seen  to  exist.  If  a  sample  of  hot  plasma 
which  exists  for  a  fairly  long  time  is  desired,  a  number  of  intermixing 
processes  must  be  examined.  The  initial  shock  attenuates  as  it  moves  down  the 
tube,  allowing' the  hot  driver  gas  to  catch  up;  in  the  meantime  energy  is  lost 
by  conduction  and  radiation,  and  some  of  the  electrons  recombine  with  ions. 

The  driver  gas  also  mixes  with  the  shocked  gas  by  turbulence  and  diffusion. 

Most  of  these  effects  are  rather  intractable  theoretically,  so  most  of  the 
knowledge  of  the  actual  conditions  in  the  tube  must  be  obtained  by  direct 
measurement,  using  probe,  microwave,  and  time  resolved  spectroscopic  techniques. 
However,  it  is  possible  to  construct  crude  models  of  the  interaction  of  the 
plug  of  hot  gas  with  the  cool  gas  in  the  tube.  Theory  and  experimental  measure¬ 
ments  both  suggest  that  the  gas  will  leave  the  discharge  with  a  velocity  of 
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the  order  of  10?  co/sec.  The  kinetic  energy  of  one  gram  of  this  gas  is: 

Yl  mv3  =  5  x  1013  ergs/gm  ~  10^  Col/gm 

For  comparison,  the  chart  below  gives  the  internal  energy  of  air  at  densities 
corresponding  to  the  pressure  range  over  which  tubes  of  this  type  are  operated 
(Reference  3-35) • 


T 

p/p0  =  10-6 

p/p9  =  10_1+ 

P/Po  =  10' 

3000 °K 

E  =  1,5  x  10 3  Cal/g 

1.4  x  103 

1.0  X  10 3 

5000 

8.5  x  10 3 

3  x  103 

1.5  x  103 

8000 

1.5  x  10 k 

1.0  x  kA 

.8  x  kA 

12000 

3.5  x  kA 

2.8  x  10^ 

1.5  x  lcA 

18000 

5.2  x  10^ 

3  x  kA 

2  x  hA 

24000 

9.2  x  kA 

6  x  kA 

3  x  104 

The  quantity  pc  is  the  density  of  air  at  standard  conditions,  1.29  x  10“3 
g/cm3.  For  air  at  50  microns  pressure  and  0°  C,  p/p0  ~  5  x  10"^. 

It  can  be  seen  that  the  kinetic  energy  of  a  given  quantity  of  gas  is 
considerably  higher  them  its  internal  energy  over  the  temperature  range 
shown  on  the  chart.  This  gives  some  justification  for  a  model  which  ignores 
detailed  thermodynamic  or  gas  dynamic  considerations  and  treats  the  moving 
plug  of  gas  as  a  moving  piston.  The  rate  of  change  of  momentum  of  the  plug  of 
gas  is  equal  to  the  rate  at  which  momentum  is  added  to  the  cold  gas: 

u  du  2 

dt 

where : 

M  =  mass  of  driver  gas  plug, 
m  «  mass/unit  length  of  shock  tube  gas, 
u  =  velocity  of  the  plug  of  hot  gas. 
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Upon  integrating: 


M 


m  (t  +  K ') 

where  K'  is  the  constant  of  integration.  With  the  initial  condition  u  =  Uq 
when  t  «  0  then: 

M 

U°~  mK ' 

1  muo 


and 


K  = 


K' 


thus: 


M 


(1  +  Kt) 


‘  ,LJL) 

\  muo  / 

This  can  be  integrated  immediately  to  give  distance  x  traveled  by  the  plug  of 

gas  from  the  exit  of  the  discharge  region  as  a  function  of  time,  hence: 

dx  uo 


or: 


dt  1  +  Kt 
In  (1  +  Kt) 


Streak  photos  made  by  other  investigators  of  actual  shocks  show  the  following 
features,  which  seem  to  support  this  picture.  The  velocity  of  advancement  of 
the  luminous  front  decreases  with  distance  from  the  discharge  region,  and 
the  rate  of  decrease  is  less  for  lower  pressures.  At  low  enough  pressure  the 
luminous  gas  fills  the  whole  of  the  tube.  The  distance  of  the  luminous  front 
can  be  fitted  reasonably  well  with  the  above  equation,  with  the  proper  choice 
of  the  parameter  K  as  shown  in  the  table  below: 


p(mm  Hg) 

uo  (  7  ) 

Msec 

K 

M/m 

3 

2 

.4 

5 

.3 

5 

.4 

12 

•  03 

5 

.05 

100 
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These  values  of  M/m  seem  reasonable,  if  it  is  assumed  that  at  lower  pressures 
most  of  the  gas  in  the  hot  plug  comes  from  evaporation  of  insulation  and  electrode 
materials  in  the  discharge  region. 

It  appears  that  there  is  sufficient  energy  available  to  heat  all  the  gas  in 
the  tube  to  fairly  high  temperature  and  a  correspondingly  high  degree  of  ion¬ 
ization;  and  that  there  exists  a  mechanism  for  distributing  the  energy,  namely 
the  transformation  of  ordered  velocity  into  thermal  velocity  by  collisions  at 
the  front  of  the  plug.  The  rate  at  which  thermal  energy  is  added  to  the  cold 
gas  can  be  determined  from  the  velocity  equation  above,  since  the  rate  at  which 
both  kinetic  and  internal  energy  is  added  to  the  cold  gas  is  equal  to  the  rate 
at  which  kinetic  energy  is  lost  by  the  driver  gas.  Detailed  calculations  of 
the  temperature  history  along  the  tube  require  further  assumptions  about  the 
details  of  the  interaction  at  the  front  and  the  energy  losses  by  radiation, 
diffusion,  wall  effects  and  so  on. 

Another  model,  which  is  probably  more  appropriate  for  making  estimates  of 
conditions  behind  the  shock,  is  that  proposed  by  Thornton  and  Cambel  (Refer¬ 
ence  3-36)  in  which  the  one  dimensional  blast  wave  theory  of  Harris  (Refer¬ 
ence  3-37)  is  modified  to  take  radiation  losses  into  account.  The  conclusion 
is  reached  that  radiation  losses  are  relatively  unimportant  in  shock  attenu¬ 
ation  except  near  the  discharge  region.  This  gives  some  justification  to  the 
calculation  of  the  conditions  behind  the  shock  by  classical  type  shock  relations. 
The  results  of  such  calculations  are  reproduced  below,  where  the  calculated 
temperature  and  concentrations  of  various  species  are  shown  as  a  function  of 
total  distance  traveled  by  the  shock.  The  model  used  assumes  that  the  energy 
of  the  gas  behind  the  shock  is  uniformly  distributed,  but  it  is  doubtful 
whether  completely  uniform  conditions  really  prevail.  (See  Figure  3-28) 

The  above  results  are  in  general  agreement  with  experimental  measurements. 
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TEMPERATURE  AND  EQUILIBRIUM  COMPOSITION 
OF  THE  AIR  BEHIND  THE  SHOCK  WAVE 
FOR  SHOCK  TUBE  CONDITIONS 

(VOLTAGE  =  17,200  V,  PRESSURE  =  30  p  Hg.) 


18,000  * 
o 


4  6  8  10  12 

DISTANCE  DOWNSTREAM  FROM  RING  ELECTRODE,  IN. 


FIGURE  3-28 


Thus,  it  can  concluded  that  the  electromagnetic  shock  tube  is  capable  of  gener¬ 
ating  highly  ionized  gases  near  the  discharge  region.  In  addition  it  is  a 
facility  that  is  easily  constructed  and  inexpensive  to  operate. 

3- ^.3*2  Arc  Discharge  Tube.  -  The  arc  discharge  tube  is  a  30  inch  long, 

5  inch  diameter  pyrex  tube,  provided  with  concentric  electrodes  at  each  end  - 
For  its  operation,  the  tube  is  evacuated  to  a  desired  pressure,  wnH  an  electric 
arc  is  struck  across  the  tube  between  the  center  electrodes.  The  arc  ionizes 
the  gas  and  maintains  this  ionization  as  long  as  the  current  continues.  This 
facility  at  McDonnell  was  utilized  in  the  same  manner  as  that  of  Wilcox  et  al 
(Reference  3-6),  who  showed  its  advantages  for  the  study  of  Alfven  waves  in 
hydrogen  and  deuterium.  Furthermore,  it  provided  an  environment  in  which 
instrumentation  techniques  could  be  economically  investigated.  The  following 
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paragraphs  discuss  the  characteristics  of  the  ionizing  process.  As  in  the 
case  of  the  electromagnetic  shock  tube,  the  basic  processes  remain,  for  the 
most  part,  undefined. 

Electrical  discharges  in  gases  are  divided  into  two  major  classes,  namely, 
self - sustained  discharges  and  non  self - sustained  discharges.  A  self- sustained 
discharge  is  one  which  maintains  itself  without  external  aid,  while  the  non 
self - sustained  discharge  depends  upon  on  an  external  source,  (e.g. ,  a  heated 
cathode)  that  produces  a  sufficient  number  of  electrons  and  ions  to  maintain 
the  discharge.  Electrical  discharges  can  be  further  subdivided  by  the  magnitude 
of  the  applied  voltage  for  the  discharge  and  the  amount  of  current  that  is 
observed.  Figure  3-2 9,  taken  from  Reference  3-27,  shows  roughly  the  afore¬ 
mentioned  subdivisions. 
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The  operating  characteristics  of  the  arc  discharge  tube  determine  the 
nature  of  the  discharge  process.  The  ionization  exists  for  about  .10  sec., 
with  an  estimated  maximum  current  of  41,000  amps.  Before  breakdown  occurs  the 
gas  in  the  tube  functions  as  an  insulator;  however,  once  the  ionization  process 
has  been  initiated,  it  becomes  a  good  conductor  and  the  impedance  of  the 
external,  circuit  will  greatly  exceed  that  of  the  discharge  tube.  Thus  the 
majority  of  the  energy  is  dissipated  in  the  external  circuitry.  This  is 
characteristic  of  an  arc  discharge  in  that  it  will  support  large  currents  with 
a  relatively  small  applied  potential. 

Arc  discharges  are  also  subdivided  (Reference  3-27)  according  to  the  emission 
process  taking  place  at  the  cathode  and  the  pressure  due  to  the  discharge. 

High  and  low  pressure  discharges  can  be  distinguished  on  the  basis  of  the 
temperature  of  the  positive  column  (References  3-27  and  3-38).  In  the  high 
pressure  arc,  the  positive  column  temperature  is  on  the  order  of  thousands  of 
degrees  Kelvin,  with  the  gas  components  in  approximate  thermal  equilibrium. 

On  the  other  hand,  in  a  low  pressure  arc,  the  electron  temperature  is  much 
larger  than  the  ion  or  gas  temperature,  on  the  order  of  tens  of  thousands 
of  degrees  Kelvin.  The  gas  pressure  for  the  McDonnell  experiments  was  in  the 
range  of  40  to  l'OOO  microns  of  Hg,  thus  satisfying  the  conditions  for  a  low 
pressure  arc. 

The  four  types  of  arcs  that  are  given  in  the  following  list  must  be  ex¬ 
amined  to  determine,  if  possible,  which  predominates.  The  thermionic  arc 
with  the  neutral  cathode  Is  easily  dismissed,  since  the  cathode  was  not  ex¬ 
ternally  heated  in  experiments;.,  Whether  the  arc  is  caused  by  thermionic 
emission  or  field  emission  depends  upon  the  electrode  material,  as  observed 
by  W.  Ramberg  (Reference  3-39)  and  noted  by  Cobine  and  Brown.  For  example, 

C,  Ca,  Mg,  and  W  are  thermionic  emitters;  while  Cu,  Hg,  Ag,  and  Au  are  field 
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emitters.  A  third  group  exists  that  exhibits  both  types  of  emission  processes; 
these  are  Pt ,  Sn,  Fe,  and  Cd.  Metal  arcs  are  distinguishable  by  observations 
of  the  cathode  spot,  the  cathode  temperature,  the  arc  response  to  a  magnetic 
field,  and  the  cathode  impurity  effects.  Loeb  (Reference  3-**0)  and  Thomson 
(Reference  3-^l)  hold  that  all  arcs  are  due  to  thermionic  emission.  Without 
a  very  detailed  examination  the  type  of  arc  is  difficult  to  determine. 

1.  Thermionic  Arc  a.  A  non  self - sustained  discharge  with  the 

cathode  being  heated  to  maintain  the 
discharge. 

b.  A  self - sustained  discharge,  i.e.,  the 
cathode  is  heated  by  the  discharge  only. 

2.  Field  Emission  Arc  The  electron  current  at  the  cathode  is 

produced  by  very  strong  electric  fields 
present  at  the  cathode  surface. 

3.  Metal  Arcs  This  classifies  metal  arcs  on  the  basis  of 

their  response  to  a  magnetic  field,  temper¬ 
ature  of  the  cathode,  and  the  influence 
of  impurities  on  the  cathode  surface, 
whether  or  not  the  cathode  spot  moves  on 
the  cathode  surface. 


The  burning  voltage  of  the  low  pressure  arc  discharge  is  made  up  of 
three  parts: 

1.  The  cathode  drop  region,  where  there  is  a  voltage  drop  on  the  order 
of  the  ionization  potential  across  a  thin  sheet  of  the  gas. 

2.  The  positive  column,  where  the  voltage  drop  is  a  function  of  the  gap 
length,  the  arc  current,  and  the  gas  pressure. 

3.  The  anode  diop  region,  which  closely  resembles  the  cathode  voltage 
profile. 


The  current  density  at  the  cathode  is  much  greater  than  that  in  the  glow.  The 
positive  arc  column  is  characterized  by  approximate  electrical  neutrality,  a 
current  principally  carried  by  electrons,  and  a  voltage  gradient  large  enough 
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to  produce  an  equal  amount  of  ionizing  collisions  to  balance  ions  lost  to  the 
walls  of  the  tube. 

The  analysis  of  the  low  pressure  arc  has  been  carried  out  for  a  number  of 
different  cases  by  Tonks  and  Langmuir  (Reference  3-^2)  and  Cobine  (Reference 
3-I+3) .  Following  the  analysis  of  Cobine,  there  are  four  independent  variables 
for  the  case  of  a  cylindrical  tube,  namely,  the  tube  radius  r,  gas  pressure 
Pg,  arc  current  ia,  and  the  wall  temperature  Tw,  which  must  be  related  to  the 
five  independent  variables,  i.e.,  the  axial  electric  field  Ez,  electron  con¬ 
centration  at  the  axis  neo,  electron  temperature  Te,  positive  ion  current 
density  at  the  wall  jj_,  and  the  number  of  electrons  generated  per  electron  per 
second  .  At  a  sufficiently  low  pressure,  such  that  the  mean  free  path  of  the 
gas  particles  is  comparable  to  the  tube  radius,  can  be  replaced  with  the  gas 
temperature  V 

The  ions  of  plasma  are  being  continuously  lost  to  the  walls  of  the  tube, 
and  it  is  assumed  that  the  ions  move  under  a  radial  potential  with  a  maximum 
at  the  tube  axis;  however,  since  the  reference  point  is  arbitrary  in  measuring 
potential,  the  zero  of  potential  may  be  taken  to  be  at  the  tube  axis  and,  thus, 
negative  everywhere  else.  The  origin  of  coordinates  is  taken  at  the  tube  axis, 
i.e.,  r  =  0.  The  ions  generated  at  some  point  z,  intermediate  between  the  axis 
and  the  wall,  will  have  attained  a  velocity  vz  at  a  distance  r  from  the  axis. 

Then  if  Nz  is  the  ion  production  rate  per  unit  volume  at  the  point  z,  the  ion 
concentration  at  r  is:  ,, 

Nz  z  dz 


To  account  for  all  ions  produced  at  distances  less  than  r  the  above  equation  if 
integrated  from  zero  to  r: 

1  r  Nz  z  dl 

"i  =  7  !o  ~~ 

Boltzmann's  relation  gives  the  expression  for  the  number  of  electrons/cm^  at 
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APPENDIX  A 

ALFVEN  WAVES  IN  A  FULLY  IONIZED  GAS 

/ 

The  following  derivation,  while  classical  in  the  sense  that  Alfven  first 
took  this  approach,  has  an  advantage  in  that  it  is  the  simplest  way  to  analyze 
transverse  disturbances  in  a  fully  ionized  gas.  This  approach  facilitates 
incorporation  of  wave  excitation  mechanisms,  total  plasma  motion  and  primary 
field  gradients  into  the  theory. 

Working  in  the  MKS  system  of  units,  the  Maxwell  equation  which  relates  the 
current  density  to  the  spatial  variation  of  the  magnetic  field  it  Induces  within 
the  plasma  and  vice  versa,  is: 


—  —  ao 

v«H-i.  — 


(A-l) 


As  usual,  j  is  the  current  density  due  to  the  field  E  +  V  x  B,  where  E  is  due 
both  to  the  plasma  constituents  and  existing  externally  caused  electric  fields 
and : 


B=B0+b, 


(A-2) 


where  Bo  is  due  to  a  constant  magnetic  field,  with  b  being  induced  by  moving 
charges  in  the  plasma.  V  is  the  total  macroscopic  plasma  velocity,  which  for 
an  initially  moving  plasma  of  velocity  V0  and  an  induced  plasma  velocity  v,  is 
given  by  V  =  \  +  v.  To  be  sufficiently  general,  it  is  necessary  to  introduce 
the  effects  of  B  on  the  plasma  conductivity.  Since  b  is  generally  found  to  be 
oscillatory,  the  average  of  B  over  many  cycles  of  b  is  just  B0,  neglecting 
harmonics.  A  common  but  more  restrictive  condition  is  to  consider  only  the  case 
where  b  <  <  BQ«  In  either  event  the  cyclotron  frequencies  are  due  to  B0,  and  J 
is  related  to  the  induced  field  through  the  conductivity  tensor  as: 
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where.  In  general: 


i=a.(E+VxB), 


al  1  CT12  a13 

<7  *  cr21  <722  ^23 

*31  *32  *33 


(A-3) 


(A-4) 


For  orthogonal  external  electric  and.  magnetic  fields  and  low  frequency  internal 
fields,  the  tensor  elements  of  £  have  the  familiar  forms: 


*11  =  *22" 


*21  “*12' 


<,eKe/l'ei>  *i  KAJ 

1+KA*;)2 


*33  *  *e  +  *i 


(A-5) 


m_  v  . 

°  ei 


with  the  remaining  elements  zero.  In  the  above  expressions  m,  e,  N  are  the 
particle  mass,  charge,  and  number  density,  with  the  subscripts  e  and  i  referring 
to  the  electrons  and  ions,  respectively;  wce  and  a>ci  are  the  electron  and  ion 
cyclotron  frequencies,  and  v  is  the  two-body  collision  frequency  for  an 
s-type  particle  with  the  r-type  particles.  For  thermal  equilibrium  of  the  gas, 

"ei  "  -ie. 

In  equation  (A-l)  the  convection  current  density  has  been  neglected,  and 
while  the  displacement  current  density  can  usually  be  neglected  with  respect  to 
the  conduction  current  density  for  low  frequency  plasma  disturbances,  it  will  be 
temporarily  considered  in  this  derivation  in  order  to  show  a  more  complete  form 
of  the  basic  differential  equation.  With  the  relationship: 

D-£.  E  (A-6) 


160 


REPORT  NO.  A219 


IV1MD  Wave  Investigation  30  November  ism 

!  I  -  theory 


the  displacement  current  density  becomes: 

~  (<-E) 

dt  dt  - 

where  the  permittivity  tensor  is  represented  by 


|€11  f12  £13 
k21  f22  f23 
k31  £32  e33 


(A-7) 


(A-8) 


Corresponding  to  the  situation  for  which  (A-5)  is  valid: 


0. 


0: 


<11-£22-2-'[1  +  ((i)ce/j/ej)2  +  1  +  (<uci/vje)2 


'  Ke^ei^c 


'("cAie)  ni 


2’  12  l  +  U)ce/vei)2  1  +  (<uc/vie)2 


(A-9) 


«33  =  2  —  i  (0e  +  flj) 

r\ 

fl_  =  6iZ  /W-j 

e  pe  81 

2 

n-,  =  a)pi/(UVie 

where  aJ  is  the  angular  frequency  of  the  disturbance  and  u^e  and  are  the 
electron  and  ion  plasma  frequencies.  The  condition  of  electrical  neutrality 
is  most  conveniently  expressed  by: 


V  -15=  0 

Equation  (A-l)  may  now  be  rewritten  as: 


—  VxB  =  <7.(E  +  VxB)+-  (e-1) 
M0  d' 


(A-10) 


(A-ll) 


where  the  constitutive  relationship  for  the  medium: 


B  =  Mo  H 


(A-12) 


has  been  used. 

Taking  the  curl  of  (A-U),  rearranging  terms  and  using  the  equation: 

?-B=0  (A-13) 
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gives: 


V2B  +  p0  V x  a  •  [fT+"v  x~B]  +  p0  —  7x(«.  E)  -  0 


at 


(a-14) 


A  second  equation  which  relates  V  and  B  and  their  derivatives  can  he  found 
from  the  Navier-Stokes  equation.  If  the  non -electromagnetic,  non-viscous 
forces  are  conservative,  the  net  force  on  a  plasma  element  (Reference  A-l)  is: 


A*  _  *  o—l  —  — 

p  DV  =  j  x  B  -  +  ri  (V2  V  +  -  VV-  V)  +  Nq*E 


(A-15) 


where : 


a—  av  -  - 

dv  =  —  +  (v .  v)v 

at 


(A-16) 

The  quantity  is  "the  hydrostatic  pressure  tensor,  in  general  a  funcation  of 

position  and  time,  0  is  the  non -electromagnetic  potential  function,  t)  is 

the  coefficient  of  viscosity,  N  is  the  total  electron  and  ion  concentration, 

and  q  is  the  net  charge  in  an  elementary  volume  of  the  gas.  The  last  term 

/  V 

in  (A-15)  contributes  only  about  (  q  )  of  the  first  terms,  c  being  the  vacuum 
speed  of  light,  and  it  is  hereafter  neglected.  Solving  equation  (A-l)  for  j, 
substituting  into  (A-15)  and  using  the  vector  identity: 


yields : 


Bx(?xB)»  —  V  i(B  •  B)- (B  .V)1T 
2 


—  +(V.V)V-  —  (B -V)B-  -  Bx  -  ((.  E) 

at  u0  p  p  at 

1  1  B2  p  77  —  n  — 

+  _  +  —  V(^+  , - ?  •  V)  -  —  V2  V -  o 

P  P  2  Po  3  P 


(A-17) 


(A-18) 


Equation  (A-l4)  and  (A-l8)  yield  a  system  of  six  scalar  differential  equations 
in  the  six  components  of  B  and  V.  To  facilitate  reduction  of  these  equations 
certain  relative  assumptions  must  be  made.  Since  for  tenuous  gases  viscosity 
is  practically  non-existent  (  y  «  0),  the  off-diagonal  terms  in  ^  are  negligible. 
If  the  pressure  components  along  the  principal  axes  are  nearly  equal. 
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the  approximation: 

V-i=Vp  (A-19) 

where  p  is  the  scalar  pressure,  is  valid. 

Also,  V  0  is  negligible  for  most  laboratory  situations.  Then  with  the  aid 
of  the  Maxwell  equation: 

m  ^  B 

Vx£-  (A- 20) 

Equations  (A-lk)  and  (A-l8)  become : 

V2B  +  fiL0Vx(2.  E)  +  MoVx[r.  (VxB)]--Mo—  Vx(oE)  (A-2l) 


(B  •  V)  B  -  (V  •  V)  V 


1  /B2  \  1  -  3  .  -  dV 

■/k*’]"'  *  (A_22) 


Due  to  the  form  of  the  conductivity  and  permittivity  tensors,  these  equations 
cannot  be  significantly  reduced,  unless  additional  assumptions  about  the  plasma 


characteristics,  such  as  incompressibility,  irrotational  flow,  no  pressure 
gradients  and  so  forth,  are  made.  The  problem  is  somewhat  simplified,  if  the 
particle  density  variations,  N,  the  magnetic  induction,  B,  and  the  plasma 
velocity,  V,  are  expressed  in  a  linear  form,  implying  that  all  terms  above 


the  first  order  are  negligible.  The  pressure  term  Vp  is  then  related  to  the 
particle  density  by: 

Vp-ykTVN  (A- 23) 

where  7  is  the  ratio  of  specific  heats,  k,  the  Boltzmann  constant,  T,  the 
equilibrium  temperature.  The  electric  field  E  is  related  to  the  linearized 
expression  for  B  through  equation  (A-20).  When  (A-20),  (A-2l)  and  (A-22)  are 
expanded,  the  result  is  a  system  of  linear  homogeneous  equations  in  B  and  the 
time  and  spatial  derivatives  of  B. 
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Alfven  specialized  equations  (A-21)  and  (A-22),  finding  that  they  implied 
an  oscillatory  motion  of  the  gas  transverse  to  the  direction  of  propagation  of 
the  disturbance.  Alfven' s  conditions  were: 

(1) 

(2) 

(3) 

(4) 


(5) 

(6) 

(7) 

With  these  conditions  (A-21)  and  (A- 22))  can  be  written  for  a  stationary 
plasma  (VQ  »  0)  as: 

~  -  —  V2^  (t  •  V)  V  -  (V  •  V)t  +  ("fig  •  V)V  (A-24) 

at  nao 

and: 

—  =  —  (b  .  vjb  -  (V  .  V)V  +  —  (B0  •  V)t  (A-25) 

<3*  W  Po  P 

These  equations  suggest  the  linear  relationship: 


The  gas  is  incompressible:.  V  *  V  =0 
The  gas  is  inviscid:  r/  =  0 

p  j  a  ,  and  «  are  scalar  constants. 

The  displacement  current  is  negligible,  that  is: 


<9E  -►  ^ 

c  - «a  (E  +  V  x  B) 

at 


or  equivalently: 


a2B 


aB 


e  -  «  a  [ - +  V  x  (V  x  B)] 


at" 


at 


All  externally  caused  non-electromagnetic  forces  are  negligible. 

/  b2  . 

( p  +  _£ — }  is  not  a  spatially  varying  function. 

2  Mo 

B  =  Bq  +  b  where  B0  is  constant  and  b  C  <  BQ. 


V  =  ab 


(A- 26) 


a  being  a  ccutplex  quantity  but  constant  in  space  and  time.  Upon  substi¬ 
tuting  (A-26)  and  differentiating  with  respect  to  time,  (A- 24)  becomes: 


aK 


at" 


I  o 

—  17* 


at 


Po° 


V"  “+«<B0.V>  ~ 


(A-27) 
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and  (A-25)  gives: 


—  (B0  .  V)t  +  (— ^ —  -  aj  (?•  Vi? 

*  afi0p  \ap0p  / 


(A-28) 


If  the  last  term  in  (A-28)  is  neglected,  then  combining  (A-27)  and  (A-28)  yields 
the  Alfven  wave  equation: 

bK  1  2  1  ♦ 


(B.  •  V)2? 


(A- 29) 


at2  w  dt  y-op  0 
A  relationship  between  the  mean  fluid  velocity  and  the  perturbed  field 
is  obtained  by  comparing  (A-27)  with  (A-29),  and  observing  that,  to  a  first 
approximation : 


-Yi 


and  therefore: 


0*„  P> 


(A-30) 


(A-31) 


Working  in  rectangular  Cartesian  coordinates  with  Bq  along  the  1  z  axis, 
the  plane  wave  solution  for  a  medium  of  infinite  extent  is: 


b-?0.i(<ut^Kz) 


(A-32) 


where  *  is  the  wave  number  and: 


b_  *  b  i 
o  ox 


(A-33) 


For  a  gas  with  very  high  electrical  conductivity,  equation  (A-29)  reduces  to: 


dK  Bo  *  2r 

7T  U  b 

giving  rise  to  an  undamped  plane  wave  of  phase  velocity: 


Bo"«^ 


(a-34) 


(A- 35) 
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called,  the  Alfven  velocity.  In  the  more  general  case  where  the  conductivity 
must  be  considered,  substitution  of  (A-32)  into  (A-29)  produces  the  dispersion 
equation: 

+,.,  .  ...  .-ft 

(A-36) 


_Zlh  •  1  \  V2 

K~  u°(  +'  Po'Ul) 


or,  in  rectangular  form: 


K=t 


U0V2p+  2 

A o 


y 

*2u l)  . 


nsju 


Mo2^2 


>  M.2*2  U< 


(A- 37) 


A  direct  measure  of  the  wave  attenuation  in  space  is  afforded  by  the 
distance  required  for  the  wave  amplitude  to  fall  to  l/e  of  its  initial  value 
b;0.  This  distance  is  referred  to  as  the  attenuation  distance  d  and  is  given  by 
the  reciprocal  of  the  imaginary  part  of  k  ;  thus,  in  rectangular  form: 


d  = 


•mW 


(A-30) 


from  equation  (A-37) •  The  phase  velocity  operator  for  the  damped  wave  is 
defined  by: 


4' 


+  I 


Mo  a  U 

The  correct  formula  for  the  phase  velocity  is  now: 

l) 


R.  («) 


Uph  R.  M 


(A- 39) 


(A-1»0) 


Since  the  constant  a  is  just  the  ratio  of  the  phase  velocity  to  the  constant 
component  of  the  field  (Equation  A- 35),  it  may  be  inferred  that.: 

-1/2 

V  (A-4l) 


.  B  _  B 
b 


-*-  ®0  /  CO  \ 

—  V-  —  (1  +  i  -  1 

u  Uo\  AaaU*/ 


Since  k  is  in  the  form  (* r  +  i  «  j),  b  can  be  written  as: 


MCi 
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(A- 42) 


where : 


cf>  =  tan  ^  -  ;  K|  <  0  <  Kp 

KR 


*1 


(A- 43) 

Thus  there  exists  a  phase  lead  of  the  ga3  velocity  with  respect  to  the  wave 
field. 

Figures  A-l  and  A-2  show,  for  various  angular  frequencies  and  plasma 
densities,  the  attenuation  distance  as  a  function  of  the  real  part  of  the  wave 
number.  The  parameter  rQ  is  the  ratio  of  the  plasma  density  to  a  reference 
plasma  density  of  1.72  x  10" 10  Kg/m3. 

The  result  obtained  in  (A-3l)  associates  the  magnetostrictive  pressure  of 
the  perturbation  field  with  the  disturbed  component  of  the  gas  velocity.  In 
addition,  observation  of  the  pinch  effect  makes  the  assumption,  that  the 
gradient  of  the  quantity  in  Alfv^n  condition  (6)  is  zero,  not  unreasonable.  The 
result  of  combining  condition  (6)  and  equation  (A-3l)  is  just  the  steady  state 
Bernoulli  equation  relating  the  mechanical  energy  density  to  the  magnetic 
energy  density. 

If  the  wave  disturbance  is  maintained  by  a  time  varying  magnetic  field  /3, 
perpendicular  to  the  main  field  Ed  and  fixed  at  a  point  in  space,  equations 
(A-24)  and  (A-25)  may  be  used  with: 


B  =  B0  +  b 
B0=B0rz 
f-  /8Tx+T' 
v  =  v  T  +"v' 

O  Z 


(A-44) 
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to  obtain  a  more  general  wave  equation. 

Thus  upon  substituting: 
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dX+  W  _  <2*2^  jP .  v)7>+  5 .  yjr+  (g.  .  v)7'  -  (V  •  V)P  (A-45) 

3t  3 1 


—  =  y2(iT'  •  V  )t'  +  y2  $  •  V)t'  -  (V  •  V)V  +  y2  (B0  •  V)b ' ’ 
<9t 


Using  the  linear  relation: 


(A- 46) 


v/=  a  b  =  a  (/3  +  bO  (A“47) 

and  differentiating  both  (A-4-5)  and  (A- 46)  with  respect  to  time  yields: 


at2  at2  at  dt 


(A-48) 


CUM  • 

~  (—  -  ^  [(b‘ '+' IS)  •  V]  b'+  B0  -  VQj  .  vjb"'  (A-49) 

Using  the  value  of  o(  obtained  in  (A-30)  and  the  result  from  (A- 35),  the 


relations : 


(A- 50) 


u0  ”  y  Bo 


(A-51) 


are  found.  Then  using  (A- 50)  and  (A-51),  the  final  form  of  the  wave  equation 


becomes: 


%  -  c2^2  -ro)  •  vl2r'“-  IT 


(A-52) 


as  would  be  expected  from  equation  (A-29).  No  solution  to  (A-52)  is  given 
here;  but  by  taking  the  LaPlace  transform  of  the  members  of  (A-52)  in  the  time 
domain  and  applying  suitable  initial  conditions,  a  transformed  equation  in  the 
spatial  coordinates  and  the  transformed  variable  is  acquired.  A  solution  of 
this  equation  with  appropriate  boundary  conditions  should  be  possible.  Then  the 
inverse  transform  of  this  solution  would  be  the  desired  re suit.  Since  the 
exact  form  of  0  has  not  been  specified,  any  Fourier  analyzable  function  can 
be  included,  and  a  solution  in  principle  can  be  obtained. 
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APPENDIX  B 

CALCULATION  OF  COLLISION  FREQUENCIES  IN  WEAKLY  IONIZED  GASES 


Although  particle  "collisions"  within  a  system  in  thermal  and  chemical 
equilibrium  do  not  affect  the  state  of  the  system,  they  are  largely  responsible 
for  propagating  local  disturbances  and  maintaining  equilibrium.  A  useful  model 
for  such  an  ionized  gas  is  the  Lorentz  gas  model.  The  assumption  of  a  perfect 
Lorentz  electron  gas  implies  that  the  electrons  are  considered  as  being  nearly 
free  of  fields  and  existing  in  a  "sea"  of  much  more  massive  neutrals.  The 
neutrals  consist  of  electrically  neutral  but  polarizable  gas  molecules.  For 
the  case  of  a  weakly  ionized  Lorentz  gas  it  is  proper  to  consider  the 
"collision"  term  in  the  Boltzmann  equation  as  being  due  to  the  interaction  of 
moving  electrons  and  stationary  neutrals.  This  "collision"  term  then  couples 
the  Boltzmann  equation  for  the  electrons  to  that  of  the  neutrals.  Del  Croix, 
Reference  (B-l),  designates  this  term  (Jf),  where  J  is  called  the  "collision 
operator"  and  f  is  the  phase  space  distribution  function  for  the  electrons  as 
it  appears  in  the  Boltzmann  equation,  (jf )  can  always,  in  principle,  be  calcu¬ 
lated  by  integrating  over  the  entire  ensemble  of  particle  interactions,  but  the 
effort  required  is  generally  prohibitive  due  to  the  complexity  of  the  integrand. 
With  the  assumption  of  a  perfect  Lorentz  gas,  however,  (Jf )  becomes  linear  in  f ; 
and,  therefore,  f  can  be  expanded  in  a  series  of  spherical  harmonics  such  as: 


fe  =  aM  +  X  (a'J’  cos  m <j>  +  /S™  sin  m<f>)  P'|’  (0);  m  .<  f 


(B-l) 


where  the  P  (0)  are  the  associated  Legendre  polynominals  of  order  1,0  and  0 
are  the  azimuth  and  colatitude  angles.  The  a  “  and  £  ®  are,  in  general, 
functions  of  time,  and  w  e>  is  the  magnitude  of  a  possible  electron  velocity. 
It  la  found  assuming  a  centrally  directed  electron- neutral  interaction,  that  J 


f. 
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commutes  with  the  spatial  rotation  operators.  This  coramutivity  property  pro¬ 
duces  the  eigenvalues  (-  v  |  )  corresponding  to  the  eigenfunctions  of  J  as 
exhibited  by: 

J  f7=  v|  f7  (B-2) 


The  Boltzmann  equation  for  electrons  within  a  constant  magnetic  field  Hq  and 
an  electric  field  E  is  given  by: 


where : 


3 

+  £  w  | 
dt  j-i  e| 

It  |  /.E  *  .  \ 

<?<?i  /I 

df. 

i  —  -  J„L 
aw.  •  • 

(B-3) 

n®* 

9 

1 

<,3| 

lo 

(B-4) 

b0  =  m0h0 

(B-5) 

The  terms  are  diffusion  terms  and  may  be  neglected  if  the  gas  inhomoge¬ 

neities  are  small  at  any  given  instant.  The  momentum  transfer  equation, 
obtained  from  (B-3)  in  the  usual  manner,  exhibits  a  term: 


%  s  fff  me  d  wel  dwe2  dw«3 


(B-6) 


[Jf]e-2(Jf7)-'S(-v,  #7);  m  ^  I 


(B-7) 


Del  Croix  points  out  that,  due  to  the  orthogonality  properties  of  the 
spherical  harmonics  in  (B-l),  only  the  anisotropy  characterized  by  m  <  I  ■  1 
makes  a  contribution  to  the  interaction  force,  which  flhr  the  electrons  is  Fe 
It  is  therefore,  possible  to  calculate  v  ,  by  means  of  classical  scattering 
formulas.  The  appropriate  expression  is  as  follows: 

n.2n  N„  V0“(l  -cosy)pdp  (B-8) 

O 

where  p  is  the  impact  parameter,  X  is  the  deflection  angle  and  Rq  Is  the 
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number  density  of  neutral  particles.  The  deflection  angle  is  determined  from: 


X  (p,  V  )  •  it-  2p  f r2 [1  -  (p/r)z  -  (r)/w0 j 

d 


,-1/2 


dr 


with: 


(B-9) 


(B-10) 


w0  =  P  Vo 

and  d  the  distance  of  closest  approach.  In  (B-10)  w0  is  the  relative  initial 
energy  of  the  two  "colliding"  particles,  with  M  being  the  reduced  mass  of  the 
two  particle  system,  and  we  =  V0  is  their  initial,  relative  speed.  In  general, 
0(r)  can  be  considered  as  the  interaction  potential  function  for  any  two 
particles  encountered;  and,  therefore,  this  analysis  is  correct  even  if  the 
"sea"  in  which  the  electrons  find  themselves  consists  of  both  neutrals  and 
ions.  The  three  component  mixture  still  employs  the  Lorentz  gas  model  for  the 
electrons,  but  a  small  fraction  of  the  neutrals  are  now  replaced  with  ions. 
Under  the  condition  where : 

Nn  »  N  =  N; 

n  e  I 

(B-ll) 


and: 


H  being  the  number  of  particles  the  type  indicated  by  the  subscript  per  unit 
volume  and,  T  the  temperature  for  that  specie,  the  interactions  that  dominate 
and  are  of  interest  in  plasma  work  are  the  electron-neutral,  ion- neutral  and 
electron- ion  "collisions".  Each  interaction  is  characterized  by  a  different 
0(r).  An  individual  0(r)  consists  of  both  a  classical  part  and  a  quantum 
mechanical  part,  the  latter  being  significant  only  for  extremely  small  approach 
distances,  d.  Insofar  as  randomization  of  particle  motion  is  concerned,  the 
distant  encounters,  d  large,  are  the  most  important.  For  distant  encounters 
then,  the  potential  function  takes  on  the  following  forms: 


Electron -neutral :  4>%n  •  - 


(  a  •  molecular  polarization) 


Ion -neutral: 

Electron -Ion: 


*ir 


*.  i~-r 


( singly  ionized  gas) 
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Due  to  its  dependence  on  the  potential  function,  v  1  has  the  exact,  inter¬ 
pretation  of  a  collision  frequency  only  for  impenetrable  particles  with  com¬ 
pletely  negligible  fields  ("billiard  balls").  Using  (B-7)  and  (B-l), 

Equation  (B-6),  written  in  rectangular  coordinates,  becomes: 

%  =  <"e  ///vl  Vo  hi  wx  +  ^ll  wy  +  “10  wz >  d  Vox  dVoy  dVoz 
The  macroscopic  interaction  force  on  the  electrons,  F  ,  produces  a  change  in 


the  total  macroscopic  velocity  Ve  according  to  the  relation: 


Fe  “  ne  me  Ve 


(B-15) 


where  v  e  is  the  tensor  of  averaged  interaction  quantities.  The  expression 


for  v  e  is : 


(B-16) 


Lv\  “11  Vod  Vo 


XX  OO 


(B-17) 


J«n  vo d  vo 


and  analogously  for  <C  v  e  >  and  <  v  > 

yy  zz 


The  coefficients  a  ^q,  a  and  0  ^  are  obtained  by  solving  the  approximate 


Boltzmann  equation: 


z  LE  *  _  \  d<e 

- +2  H+V0xa)J  -  =  Jf„ 

*  i=x  \me  °  /  i  Moi  • 


(B-18) 


and  retaining  only  the  first  order  terms  in  (B-l).  These  coefficients  are,  in 
general,  functions  of  vp  VQ,  B  and  E.  Obviously,  if  is  independent  of 
Vo,  —  e  tec ernes  the  scalar  collision  frequency  v  e  for  the  electrons. 
Equations  analogous  to  (B-l4)  through  (B-l8)  can  be  written  for  the  ions,  the 
essential  difference  being  the  tendency  for  the  ion- neutral  interactions  to  be 
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of  the  "billiard  ball"  type. 

An  alternative  classical  technique  for  calculating  collision  frequencies 
is  effected  by  means  of  the  differential  cross  section,  o'  ;  This  quantity  is 
defined  by: 


(B-l 9) 


where : 


(B-20) 


is  the  cross  section  of  the  incident  particle  beam,  p  being  the  impact  parameter, 
and  d  fi  is  the  element  of  solid  angle  into  which  a  number  of  the  particles  axe 


scattered.  In  terms  of  the  deflection  angle: 

d0=  2rr  sin  y  d  y 


(B-21) 


a  dj  dfi  _  P  dp 
dp  dfl  sin  y  dy 


(B-22) 


For  the  electron- ion  interaction,  with  fie±  from  (B-13),  has  the  form: 


(B-23) 


which  is  Rutherford's  scattering  formula.  The  deflection  angle  is  calculated 
from  (B-9),  provided  d  can  be  determined.  The  total  scattering  cross  section  is 


given  by: 


4r  v  /dPj  \ 

dX|j 


(B-24) 


Assuming  a  Maxwellian  distribution  of  initial  speeds,  V0,  the  kinetic  theory 
of  gases  produces  the  following  expression  for  the  toted  number  of  two  particle 
encounters  per  unit  volume  and  time : 

/.,\  3/2  /2\  1/2  *  ^  _  /,vV2kT 

«ii  - . . .  (7)  'ii,v«)dv«  (B-25) 
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Kj_ j  *  l/2  for  like  particles  (i  =  j) 

1  for  unlike  particles  (i  4  j) 


(B-26) 


ni  and  n<  are  the  number  densities  of  particle  types  i  and  j,  and  the  other 
quantities  are  as  previously  defined..  The  collision  frequency  is  directly 
related  to  (B-25)  by  means  of: 

N'ii 

(B-27) 


V'i=  "nT 


The  total  number  of  collisions  per  unit  time  per  type  i  test"  particle  is 


2Nii  +  Nij 


(B-28) 


since  for  a  two-particle  collision  involving  the  same  type  of  particles,  either 
one  is  the  "test"  particle.  In  (B-25)  the  proper  determination  of  0^  (V0) 
is  through  quantum  mechanical  considerations;  however,  a  common  technique  is 
to  assume  an  average  value  for  ^jj (v0),  commute  it  with  the  integral  operator 


and  integrate  the  remaining  expression.  The  result  of  this  procedure  is 

2K 


Nii  = 


2Kij  /  2kT  _ 

yj  _  Hj  Rj  CTjj 


Then: 


2K 


Vij  = 


2Kij.  /  2kT 

—  V  —  ni-ii 


(B-29) 


(B-30) 


Since  the  most  probable  speed  for  a  Maxwellian  distribution  of  particles  of 
mass  ^  is:  ___ 

C=_LJH!  (B-31) 

VT  ^ 

and  since  the  classical  expression  for  the  mean  free  path  of  a  particle  of 
mass  m  is: 


Lj.  -■ 


1 


n;  <7 


(B-32) 


•I 
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